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On the Smarandache additive sequence 

Fanbei Li 

School of Mathematics and Statistics, Inner Mongolia Finance and Economics College 
Hohhot 010051, Inner Mongolia, P.R. China 

Abstract In this paper, we use the elementary method to study the properties of the Smaran- 
dache additive sequence, and give several interesting identities for some infinite series involving 
the Smarandache additive sequence. 

Keywords Smarandache additive sequence, infinite series, convergent properties. 



§1. Introduction and results 

For any positive integer m > 2, let 1 < d\ < d 2 < ■ • • < d m are to positive integers, 
the Smarandache multiplicative sequence A m is defined as: If d\, d 2 , ■■■ , d m are the first 
to. terms of the sequence A m , then dk > dk- i , for k > m + 1, is the smallest number equal 
to d ^ 1 ■ d ^ 2 where cq > 1 for all i = 1, 2, ••• , in. For example, the Smarandache 

multiplicative sequence A 4 ( generated by digits 2, 3, 5, 7 ) is: 

2, 3, 5, 7, 210, 420, 630, 840, 1050, 1260, 1470, 1680, 1890, 2100, 

In the book “Sequences of Numbers Involved Unsolved Problems”, Professor F. Smarandache 
introduced many sequences, functions and unsolved problems, one of them is the Smarandache 
multiplicative sequence, and he also asked us to study the properties of this sequence. About 
this problem, Ling Li [3] had studied it, and proved the following conclusion: 

For any positive integer to. > 2, let 1 < di < d 2 < • • • < d m are m positive integers, and 
A m denotes the Smarandache multiplicative sequence generated by d\ , d 2 , ■ • • , d m . Then for 
any real number s > 0, the infinite series 

1 

I s 

is convergent, and 




E 

n£A„ 



1 

n s 



i — 1 i — 1 a * 



In this paper, we define another sequence called the Smarandache additive sequence as 
follows: Let 1 < d\ < d 2 < • • • < d rn are m positive integers, we define the Smarandache 
additive sequence D rn as: If di, d 2 , ■ ■ ■ , d rn are the first m terms of the sequence D m , then 
dk > dk- i, for k > to + 1, is the smallest number equal to ou • d\ + a 2 • d 2 + + a m • d m , 
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where cti > 1 for all i = 1, 2, • • ■ , to. For example, the Smarandache additive sequence D 2 ( 
generated by digits 3, 5 ) is: 

3, 5, 8, 11, 13, 14, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 



About the properties of this sequence, it seems that none had studied it yet, at least we 
have not seen any related papers before. The problem is interesting, because it has the close 
relationship with the positive integer solutions of the indefinite equation 

di ■ X\ + d ,2 ■ X 2 + da ■ X 3 + ■ ■ ■ + d m -i ■ x m -i + d m ■ x m = n. 



The main purpose of this paper is using the elementary method to study the convergent proper- 
ties of some infinite series involving the Smarandache additive sequence, and get some interesting 
identities. For convenience, we use the symbol a m (n) denotes the n-th term of the Smarandache 
additive sequence D m . In this paper, we shall prove the following: 

Theorem. For any positive integer m > 2, let 1 < d\ < < ■ • • < d m are to positive 

integers, and D rn denotes the Smarandache additive sequence generated by di, d 2 , d m . 
Then for any real number s < 1, the infinite series 



E 

ntzDm 



1 

n s 



E 




(i) 



is divergent; For any real number s > 1, the series (1) is convergent, and 



E 

n£Dm 



E 

n— 1 



,W 



= C(s) + 



E 

a rn {n)<M 



— ~ y — 

» n<M ^ 



where £(s) is the Riemann zeta-function, and M is a computable positive integer. 

From this Theorem we may immediately deduce the following three corollaries: 

Corollary 1. Let D 2 be the Smarandache additive sequence generated by 2 and 3, then 
we have the identity 

y, 1 _ 7T 2 157 

^ al(n) 6 144' 

n = 1 7 



E 



n * 



Corollary 2. Let D 2 be the Smarandache additive sequence generated by 4 and 6, then 
we have the identity 



E 

nG£>2 



1 



n * 



E 



1 

al(n) 



7 r 2 157 

24 " 576' 



Corollary 3. Let D3 be the Smarandache additive sequence generated by 3, 4 and 5, then 
we have the identity 



E 

ne£> 3 



E 



1 

o§(n) 



7T 2 111111 1 1 1 

~6~ ~ ~ 4 ~ 36 ~ 49 _ 64 ~ 81 _ 100 _ 121 ~ 169 _ 196' 
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§2. Proof of the theorem 

To prove our Theorem, we need a simple Lemma which stated as follows: 

Lemma. Let a and b are two positive integers with (a, b) = 1, then for any positive 
integer n > ab — a — b, the indefinite equation 

ax + by = n 



must has nonnegative integer solution. 

Proof. See Theorem 19 of Chapter one in reference [3]. 

Now we use this Lemma to prove our Theorem. First note that for any positive integer 
n > m, we have 

o-m( n ) = oq • d\ + «2 ’ d-2 + ■ ■ ■ + a m ■ d m , 
where Oj > 1 for all i = 1, 2, 3, ■ ■ • , m. 

From Lemma and its generalization we know that there exist a positive integer M such 
that for any integer n > M, the indefinite equation 

d\ • X\ + di ■ X 2 + c ?3 • X 3 + ■ ■ ■ + d m - 1 • x m —i + d m ■ x m = n 



must have positive integer solution. So from the definition of the sequence D m we know that 
for any integer n > M, there must exist one term a m (k) in sequence D m such that n = a m (k). 
Therefore, 



E 

n^Dm 



1 

n s 



E 




E 



E 



00 1 1 

E^- E E 



n° * — ' n° 

n—l n<M 



CW - E 

n<M 



a rn {n)<M 



M) 



E 



1 



a rn {n)<M 



>y 



where £(s) is the Riemann zeta- function, and the infinite series 




n—l 



is convergent if s > 1; 



It is divergent if s < 1. This completes the proof of Theorem. 

Now we prove our Corollary 1. Taking a = 2, b = 3, note that a-b — a — b = 2 • 3 — 2 — 3 = 1, 
so from Lemma we know that for all integer n > 2, the indefinite equation 2x + 3y = n has 
nonnegative integer solution. Therefore, for all integer n > 7, the indefinite equation 2x+3y = n 
has positive integer solution. So M = 7, from our Theorem we may immediately deduce that 



E 



1 

a\ (n) 



c(2)+ £ 



a m ( n )<7 






E; 

n< 7 



7T 2 1 1 1 1 

T + 4 + 9 + 25 ~ 

n<6 

7 r 2 157 

T ~ 144’ 
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7 r 



where we have used the identity £(2) = — . This proves Corollary 1. 

To prove Corollary 2, note that the greatest common divisor of 4 and 6 is (4, 6) = 2, so 
from Corollary 1 we have the identity 



^ 1 if.. 11 1 v- 1 ) 

E^) “ 4F (2)+ 4 + 9 + 25~E^} 

7 r 2 157 

24 _ 576' 

Now we prove Corollary 3. Note that for all n > 15, the indefinite equation 



3x + 4y + 5z = n 

has positive integer solution. So we can take M = 15. Therefore, 



E 



i 

a§(n) 



C(2) 



« 3(^)<15 



a§(n) 



- £ 72 



n< 15 



7T 2 1 1 1 1 1 

~6 ~ ~ 4 ~ 36 ~ 49 _ 64 _ 81 



1 

Too 



1 

121 



1 

169 



1 

196' 
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Contributions to classical differential 
geometry of the curves in E 3 

Melih Turgut and Suha Yilmaz 

Dokuz Eylul University, Buca Educational Faculty 

Department of Mathematics, 35160, Buca-Izmir, TURKEY 
E-mail: melih.turgut@gmail.com 

Abstract In this paper, a system of differential equation whose solution gives the components 
on the Frenet axis in E 3 is established by means of Frenet equations. In view of some 
special solutions of mentioned system, position vector of rectifying curves, osculating curves 
with constant first curvature, normal curves and special cases are presented. Moreover, a 
characterization for inclined spherical curves is expressed. 

Keywords Euclidean space, Frenet frame. 



§1. Introduction 

It is safe to report that the many important results in the theory of the curves in E 3 were 
initiated by G. Monge; and G. Darboux pionnered the moving frame idea. Thereafter, F. Frenet 
defined his moving frame and his special equations which play important role in mechanics and 
kinematics as well as in differential geometry (for more details see [2]). 

In the Euclidean space E 3 , it is well-known that to each unit speed curve with at least four 
continuous derivatives, one can associate three mutually orthogonal unit vector fields T, N and 
B are respectively, the tangent, the principal normal and the binormal vector fields [3]. At each 
of the curve, the planes spanned by {T, N}, {T,B} and {N, B} are known respectively as the 
osculating plane, the rectifying plane and the normal plane[5]. The curves p : I C R — > E 3 for 
which the position vector p always lie in their rectifying plane, are simplicity called rectifying 
curves [1]. Similarly, the curves for which the position vector p always lie in their osculating 
plane, are for simplicity called osculating curves and the curves for which the position vector 
always lie in their normal plane, are for simplicity called normal curves [4]. 

In this paper, position vectors of some special curves in E 3 are investigated. First, a system 
of differential equation is established. In view of some special solutions, characterizations of 
some special curves are presented. 
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§2. Preliminaries 

To meet the requirements in the next sections, here, the basic elements of the theory of 
curves in the space E 3 are briefly presented. (A more complete elementary treatment can be 
found in [5].) 

The Euclidean 3-space E 3 provided with the standard flat metric given by 



(, ) = dx\ + dx 2 + dx%, 

where (cci, X 2 , £ 3 ) is a rectangular coordinate system of E 3 . Recall that, the norm of an arbitrary 
vector a G E 3 is given by ||a|| = \J (a, a), p is called an unit speed curve if velocity vector v of 
p satisfies ||u|| = 1. For vectors v,w G E 3 it is said to be orthogonal if and only if (v,w) = 0. 
Let d = d(s) be a curve in E 3 . If the tangent vector field of this curve forms a constant angle 
with a constant vector field [/, then this curve is called an inclined curve. The sphere of radius 
r > 0 and with center in the origin in the space E 3 is defined by 



S 2 = {p= (p 1 ,p 2 ,p 3 ) G E 3 : (p,p) = r 2 }. 

Denote by {T, N, B} the moving Frenet frame along the curve p in the space E 3 . For an 
arbitrary curve p with first and second curvature, n and r in the space E 3 , the following Frenet 
formulae are given in [3], [5] 

‘ Jr/ 

N' 

B’ 

where 

(T,T) = (N,N) = (B,B) = 1, 

(T, N) = (T, B) = (T, N) = (N, B) = 0. 



1 

0 

se 

0 




T 


0 

se 

1 




N 


0 -T 0 




B 



Here, first and second curvature are defined by k = k(s) = ||T'(s)|| and r(s) = — ( N,B ’). 

It is well-known that for a unit speed curve with non vanishing curvatures following propo- 
sition holds [3]. 

Proposition 2.1. Let p = p(s) be a regular curve with curvatures k and r, then; 

i. It lies on a sphere ( S 2 ) if and only if 

r 

h 

ii. It is an inclined curve if and only if 

/i 

— = constant. 




T 



(3) 
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§3. Position vectors of some special curves in E ?> 

Let ip = <p(s) be a unit speed curve in E 3 . We can write this curve respect to Frenet frame 
{T,N,B} as 



<p = <p(s) = m\T + TO2-/V + m^B. 



(4) 



Differentiating (4) and considering Frenet equations, we have a system of differential equation 
as follow: 



m'i — to 2 « — 1 = 0 
m' 2 + rniK — to 3 t = 0 

TOg + TO2 t = 0 



(5) 



This system’s general solution have not been found. Due to this, we give some special values 
to the components. 

Case 3.1. mi = ci=constant. From (5)i we have second component 

1 



to 2 = 



1 / lV 



And (5)2 gives us third component as 



m 3 = -- - + Ci-. 

r \k J t 

m 2 and m3 satisfie (5) 3 , if we rewrite it, we have following differential equation: 



- (- 
T \ K 



a 



Cl - --=0. 



( 6 ) 



(7) 



(8) 



Considering (8) and Proposition 2.1., suffice it to say that following theorem holds: 

Theorem 3.1. Let ip = <p(s) be a unit speed curve and first component of position vector 
on Frenet axis be constant in E 3 : 

i. If first component is zero, then (p is a spherical curve (lies on S 2 ). 

ii. If <p is an inclined curve, then ip lies on S 2 . Therefore, ip is an inclined spherical curve. 

iii. All normal curves in the space E 3 are spherical curves. 

iv. Position vector all of normal curves in E 3 , can be written as follow: 

1 riY 



1 

ip = — N — 

K T \ K 



B. 



(9) 



Here, it is not difficult to see that statement ii is a characterization for inclined spherical 
curves in E 3 . 

Case 3.2. m2 = C2 = constant 7^ 0. From equation (5)i, we obtain 

m'i — 1 



From (5)2, we write 



C2 = 



TO3 = TOi — . 

T 



(10) 
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Now, let us suppose ip is an inclined curve. In this case, If (10) and (11) are used in (5)3, we 
have a differential equation as follow: 



, ,K T. T 

mi(- + -) = 0. 

T K K 



From solution of (12), we have the components 



to 1 = 



t 2 s 

n 2 + t 2 



(12) 



(13) 



and 



to 2 




(14) 



K 

TO 3 — ( 



T K Z + T Z 



l), 



(15) 



where l £ R. 

Corollary 3.2. Let ip = <p(s) be a unit speed curve and second component of position 
vector on Frenet axis be constant in E 3 . Then position vector of ip can be formed by the 
equations (13), (14) and (15). 

Case 3.2.1. m 2 = 0. From (5), we have, respectively, 



TOi = s + l 



(16) 



and 

to 3 = — (s + l) = constant, (17) 

T 

where l £ R. 

Considering (16) and (17), we give following theorem. 

Theorem 3.3. Let p = ip(s) be a unit speed rectifying curve in E 3 . Then; 

i. There is a relation among curvatures as 

K 

— (s + t) = constant. 

T 

ii. Position vector of ip can be written as follow: 

p = ip(s ) = (s + l)T + - (s + /) B, (18) 

T 

where l £ R. 

Remark 3.4. The case to 3 = constant ^ 0 is similar to case 3.2. 

Case 3.3. to 3 = 0. Thus r = 0. And here, let us suppose n = constant. Using these in 
the system of differential equation, we have second order homogenous differential equation with 
constant coefficient as follow: 



m '2 + to 2 k 2 + k = 0. (19) 

Solution of (19) is elementary. Then, we write the components 



toi = A 2 sin ks — A\ cos ns 



(20) 
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and 

TO 2 = A\ cos ks + A 2 sin ks H — , (21) 

K 

where Ai and A 2 are real numbers. 

Corollary 3.5. Let cp = ip(s) be a unit speed osculating curve with constant first curvature 
in E 3 . The position vector of ip can be composed with the equations (20) and (21). 
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Almost everywhere stability of the 
third-ordered nonlinear system 

Kaiming Wang 

School of Science, Chang’an University, Xi’an, Shaanxi, P.R. China 
E-mail: kaimingwang@yahoo.com.cn 

Abstract In order to study the stability of the third-ordered nonlinear system, density function 
is constructed, and an inequality condition is obtained. When the considered system satisfies the 
inequality for almost all initial value, its equilibrium is asymptotically stable almost everywhere. 
Keywords Amost everywhere stability, density function, nonlinear system. 



§1. Introduction 



Stability is a very important factor in the research of high-ordered nonlinear systems. The 
techniques currently available in the literature are mostly based on Lyapunov’s direct method 

[1], contraction mapping method [2], nonlinear measure method [3,4], nonlinear spectral radius 
method [5] and so on. All these methods make the research of stable property everywhere. And 
if the equilibrium of the considered system is asymptotically stable, there exists a neighborhood 
of the equilibrium which is called attraction region. In fact, many real nonlinear systems do not 
need its equilibrium acting as an attracting actor attracting any points of its neighborhood, and 
it only needs attract “enough many” points in this neighborhood. That is, the equilibrium is 
asymptotically stable almost everywhere (the equilibrium attract all points of this neighborhood 
except a set with zero Lebesgue measure). In 2001, Anders Ranzter [6] proposed the concept 
of a “dual” to Lyapunov’s stability theorem, and he called the stability derived by this way as 
stability almost everywhere. In control theory, most control systems are in high order, such as 
robotic multi-fingered grasping system [1], neural networks and so on, and up to their research 
background, it’s unnecessary or impossible for the equilibrium to attract every points in its 
neighborhood even it’s very small. So it’s valuable to study the almost everywhere stability of 
high-ordered systems. 

In this paper, the almost everywhere stability of the third-ordered system is studied. Den- 
sity functions are constructed as a tool to analyze the stability of the given system. 

For convenience, the following notations were used in this paper: 

(1) Vp = gradp=(^ JiL ... p:R n^R, 

(2) V .F = divF’=fi + ff + ••• + §£, F : R n —> R n ; 



(3) represents inner product. 
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§2. Main result and its proof 

Consider a third-ordered nonlinear system 

x + ip(x,x)x +f(x,x) = 0, /( 0,0) = 0, 

or its equivalent system 



1 x N 




( y ^ 


y 


= 


z 


V * ) 




\ -^(x,y)z- f(x,y) ) 



(1) 



(2) 



It is assumed that solutions of (1) exist and are unique. 

Theoretically, this is a very useful nonlinear system since (1) is a rather general third- 
ordered nonlinear differential equation. For example, many third-ordered differential equations 
which have been discussed in [8] are special cases of (1). However, most available research all 
focused on the study of everywhere stability and it is not necessary or impossible for many real 
nonlinear equations. 

In this section, the almost everywhere stability of third-ordered system is studied. Density 
functions are constructed as a tool to analyze the stability of the given system, and the condition 
was obtained that the given nonlinear system is asymptotically stable almost everywhere. That 
is, we will prove the following main result: 

Theorem. If there exists a positive real number a, such that the following inequalities 

(1) xf(x,y)> 0, yf{x,y)> 0; 

(2) ijj(x,y)>0; 

( 3 ) Jo (/(*> °) + y&xix, v)v + yf x {x, v ) + z^(x, y))dv+ 

+a~ 1 ip(x, y) (fj f(u, 0)dit + (vip(x, v) + f(x, v))dv)) 

< i’(x,y)z 2 { 1 - ^« _1 ) 

hold for almost all x, then the equilibrium of the nonlinear system (2) is asymptotically stable 
almost everywhere. 

First, we need a basic Lemma that is necessary in proof of Theorem. 

Lemma. Given the equation x(t) = /(x(f)), where / £ C 1 {R n , R n ) and /( 0) = 0, suppose 
there exists a non-negative scalar function p £ C 1 (R n \ {0},i?) such that p(x)f(x)/\x\ is inte- 
grate on {a; £ R n : |x| > 1} and [V • (fp)](x) > 0 for almost all x. Then, for almost all initial 
states x(0) the trajectory x(t) exists for t £ (0, oo) and tends to zero as t — > oo. Moreover, 
if the equilibrium x = 0 is stable, then the conclusion remains valid even if p takes negative 
values. 

Proof. ( See reference [6] ). 

By this Lemma, now we shall complete the proof of Theorem. 

The proof of Theorem. 

Let x = (x, y 1 z) T and 



12 



Kaiming Wang 



No. 3 



then V(x) > 0 and when i / 0 we have V(x) ^ 0. That is, V(x) is positive definite. Let 
p = V~ a , then p e C 1 (i? 3 \{0}, R) and p(x)F(x)/\x\ is integrable on {x G R 3 : \x\ > 1}. 

For di vF = —^(x, y) and 

[V • (Fp)} (5) = V"- 1 (£) (VV • F - a W • F) (x ) , 

thus the condition [V • (Fp)](i) > 0 in Lemma will be satisfied if and only if (VV • F — aVL • 
F)(x) > 0. From 

VF.F 

= v x x' + V y y' + V z z' 

ry ry 

= f(x,0)y + y ip x (x,v)vdv + y f x (x,v)dv - tp(x,y)z 2 + i>(x,y)yz, 



a _1 FV • F 



rx ry 

= — a~ 1 ip(x,y ) / f(u 1 d)du — a~ 1 %l){x,y) / vijj(x,v)dv 

Jo Jo 

ry i 

-~a~ x i>{x,y) J f(x,v)dv - -a~ 1 z 2 tp{x,y), 

then by Lemma, we know that if the inequality 

[ (f(x, 0) + yi> x (x, v)v + yf x (x, v) + zip(x, y)) du 
Jo 

/ rx ry 

+a~ 1 ip(x, v) yj f(u,0)du + J {vip(x, v) + f(x,v))dv) 



< 1 P(x,y)z 2 (l — -a x ) 



holds, the equilibrium of the considered system is asymptotically stable almost everywhere. 
This completes the proof of Theorem. 

Simplifying the condition of the inequality in Theorem, we may immediately derive the 
following two corollaries. 

Corollary 1. If there exists a positive real number a > -, such that two functions f{x,y ) 
and ip{x,y) satisfy the inequalities 

(1) xf(x,y)> 0, yf(x,y)> 0; 

(2) ip{x,y)> 0; 

(3) fg(f(x,0) + yip x (x,v)v + yf x (x,v) + ip{x, y)z)dv < 0; 

( 4 ) Jo* f(u,0)du + fg(vip(x,v) + f(x,v))dv < z 2 (a- i) 

for almost all x, then the equilibrium of the nonlinear system (2) is asymptotically stable almost 
everywhere. 

Corollary 2. If there exists a positive real number a > - ,such that two functions f(x, y) 
and ip(x,y) satisfy the inequalities 

(1) xf(x,y)> 0, yf(x,y)> 0; 

(2) i/j(x,y)> 0; 
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(3) yf(x,0) + yzil>(x,y) + y 2 (ip x (x,£)£, + f x (x,ii)) <0 ,V^,ry£ [0,y]; 

(4) xf(s,0) +y(hif>(x,h) + f(h, 0)) < z 2 (a - ^),Vs £ [0 ,x\,he [0,y] 

for almost all x, then the equilibrium of the nonlinear system (2) is asymptotically stable almost 
everywhere. 

Since the system (1) and system (2) are equivalent systems, and the equilibrium x* = 
0 of system (1) answers to the equilibrium x* = 0 of system (2), their stability properties 
are also equivalent. Thus the equilibrium of the third-ordered nonlinear system (1) is almost 
asymptotically stable when the inequality condition of Theorem is satisfied. 
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Abstract For any fixed positive integer m, if there exists a positive integer n such that 
n = m - d(n), then m is called an index of beauty, where d(n ) is the Dirichlet divisor function. 
In this paper, we shall prove that there exist infinite positive integers such that each of them 
is not an index of beauty. 
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§1. Introduction and result 

For any positive integer n, the famous Dirichlet divisor function d(n) is defined as the 
number of all distinct divisors of n. If n = p“ L • p 2 2 • • • p^ k be the prime power factorization of 
n, then from the definition and properties of d(n) we can easily get 

d(n) = («i + 1) • (a 2 + 1) • • • (a k + 1)- (1) 

About the deeply properties of d(n), many people had studied it, and obtained a series results, 
see references [2] , [3] and [4] . In reference [5] , Murthy introduced an index of beauty involving 
function d{n) as follows: For any positive integer to, if there exists a positive integer n such 
that 

n 

m= d^n)’ 

then m is called an index of beauty. At the same time, he also proposed the following conjecture: 
Conjecture. Every positive integer is an index of beauty. 

Maohua Le [6] gave a counter-example, and proved that 64 is not an index of beauty. We 
think that the conclusion in [6] can be generalization. This paper as a note of [6], we shall use 
the elementary method to prove the following general conclusion: 

Theorem. There exists a set A including infinite positive integers such that each of n € A 
is not an index of beauty. 

§2. Proof of the theorem 

In this section, we shall prove our Theorem directly. For any prime p > 5, if we taking 
to = p p_1 , then we can prove that to is not an index of beauty. In fact, if to is an index of 
beauty, then there exists a positive integer n such that n = p^ 1 • d(n). From this identity we 
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can deduce that p p 1 | n. Let n = p a ■ b and (p, b ) = 1, then from n = p p 1 • d(n) and (1) we 
have the identity 

p a ~ p+1 -b={a + l) -d(b). (2) 

From (2) we may immediately deduce that a > p — 1. If a = p — 1, then (2) become b = p • d(b), 
this contradiction with (6, p) = 1. If a = p, then (2) become p ■ b = (p + 1) • d(b), or 

6= (l + £) -d(6). (3) 

It is clear that (3) is not possible if prime p > 5. In fact if b = 1, 2, 3, 4, 5, 6, 7, then it is 
easily to check that (3) does not hold. If b > 8, then from (1) and the properties of d(b) we can 

deduce that b > ^ • d(b) > \ 1 + - ) • d(b). So (3) does not also hold. So if a = p, then (2) is 
2 \ PJ 

not possible. If a > p + 1, let a — p = k > 1, the (2) become 

p k+1 ■ b = (p + k + 1) ■ d(b). (4) 

Note that b > d(b) for all positive integer b , so formula (4) is not possible, since p k+1 > p+k+1 
for all prime p > 5 and k > 1. Since there are infinite prime p > 5, and each p p ~ x is not an 
index of beauty, so there exists a set A including infinite positive integers such that each of 
n £ A is not an index of beauty. This completes the proof of Theorem. 
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Abstract For any positive integer n, let 4>{n) denotes the Euler function. That is, <j>(n) 

denotes the number of all positive integers not exceeding n which are relatively prime to 

n. g(n) denotes the smallest positive integer m such that n \ <l>(n) ■ m. The main purpose 

of this paper is using the elementary method to prove that for any positive integer group 

/ k \ k 

(mi, m 2 , • • • , mi-), we have the inequality g I P rrii ] < O g(rm), and also for any positive 

\i=l / i= 1 

integer k, there exist infinite group positive integers (mi, m 2 , • • • , mk) satisfying the equation 

/ fc \ k 

9 E mi p(mj) . 

\i= 1 / 2=1 

Keywords Euler function </>(n), three prime theorem, elementary method. 



§1. Introduction and results 

For any positive integer n, the famous Euler function <j>(n) is defined as the number of 
all positive integers not exceeding n which are relatively prime to n. Now we define a new 
arithmetical function 5(71) as follows: g(n ) denotes the smallest positive integer m such that 
n | <j>{n) ■ m. That is, g(n) = min {m : n \ <j)(n) ■ m}. For example, the fist few values of g(n) 
are: g{ 1) = 1, g{ 2) = 2, 5 (3) = 3, g{ 4) = 2, g( 5) = 5, g{ 6) = 3, g{ 7) = 7, 5(8) = 2, g(9) = 3, 
ff(10) = 5, 5 (11) = 11, 5(12) = 3, 5(13) = 13, 5(14) = 7, • • • . 

About the elementary properties of g(n), it seems that none had studied it before. From the 
definition of g{n) we can easily deduce that for any prime p, g(p a ) = p and P(n) < g(n) < Y[p, 

p\n 

where P(n) denotes the largest prime divisor of n. If n > 1 and n 2“, then g(n) must be an 
odd number. g(n ) is also a square-free number. 

The main purpose of this paper is to study the elementary properties of g(n), and prove 
the following two conclusions: 

Theorem 1. For any positive integer k and any positive integer group (mi, m2, • • • , m*,), 
we have the inequality 

( k \ k 

]d m * ) < I \g{m.i). 

i = 1 / i = 1 

Theorem 2. For any positive integer k > 3, there exist infinite group positive integers 
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(mi, m 2 , ■ • • , rrik ) such that the equation 




Whether there exist infinite group positive integers (to 1 , to 2 ) satisfying the equation 



g(mi + to 2 ) = 5(toi) + g(m 2 ) (1) 

is an unsolved problem. If mi > 1 and ?n 2 > 1 satisfying the equation (1), then one of TOi and 
m 2 must be a power of 2. 

It is an interesting problem to study the mean value properties of g{n). For any real 
number x > 1, we conjecture that 

^g(n) ~C-z 2 , 

n<.x 

where C is a computable constant. 



§2. Proof of the theorems 

In this section, we shall use the elementary method to prove our Theorems. First we prove 
Theorem 1. For any positive integer to, since g(m) = g(m). So Theorem 1 holds if k = 1. 
Now we assume k = 2. For any positive integers to and n, if (to, n) = 1, then we have 
< t>(mn ) = <p(m)(f>(n ). This time let g(mn) = t, g(m) = u, g(n) = v. From the definition of g(n) 
we have mn \ <p(mn) ■ t, m | <p(m) ■ u, n \ <fi(n) ■ v , so we have 

mn | </>(to) ■ <p(n) ■ uv = <p(mn) ■ uv. 

By the definition of t we have t < uv. Therefore, g(mn) < g(jn) ■ g(n). 

If (to, n) > 1, without loss of generality we can assume that to = p a and n = yr , where 
p be a prime, a and /3 are two positive integers, then from the definition of g{n) we have 
g(p a pP) = p, g(p a ) = p , g(pP) = p. So g{p a -p^) < g{p a ) ■ g{p From the above we know that 
for any positive integers to and n, we have the inequality 

g{mn) < g(m) ■ g{n). (2) 



Now for any positive integer k > 2 and any positive integer group (mi, m 2 , ••• , m*,), 
from (2) we have 



9 




9 




/k - 1 \ k 

< 9 ( n rn ) • g(m k ) < ■ ■ < JJg(TOj). 

\ 2=1 / i= 1 



This proves theorem 1. 

To prove Theorem 2, we need the famous Vinogradov’s Three Prime Theorem, which was 
stated as follows: 
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Lemma 1. There exists a sufficiently large constant K > 0 such that each odd integer 
n > K can be written as a sum of three odd primes. That is, n = p± + P 2 + P 3 , where pi 
(i = 1, 2, 3) are three odd primes. 

Proof. (See reference [1]). 

Lemma 2. Let k > 3 be an odd integer, then any sufficiently large odd integer n can be 
written as a sum of k odd primes 



n = px+p2-\ 1 ~Pk- 

Proof. (See reference [2]). 

Now we use these two Lemmas to prove our Theorem 2. If k = 3, then from Lemma 1 we 
know that for any prime p large enough, it can be written as a sum of three primes: 

P = P1+P2+P3- 

Now for any positive integer k > 4 any prime p large enough, if k be an odd number, then 
by the definition of g{n) we know that g{p) = p. From Lemma 2 we may immediately get 

p = g{p) = g{p\+P 2 h \-pk) = Pi +P 2 h \-pk = g(pi) h \-g(pk)- 

If k > 4 be an even number, then k — 1 > 3 be an odd number. So if prime p large enough, 
then p — 2 can be written as a sum of k — 1 primes: 

p-2 = pi+p 2 -\ \~Pk-i- 

This implies 

p = 2 + p-L+p 2 -\ b Pk - 1 

or 

P = g{p) = g(2+p 1 +p2-\ \-pk-i) = 2 + pi +P 2 H b Pk—i = g(2)+g(p 1 ) H Vg{p k -i)- 

Since there are infinite primes p , so there exist infinite group positive integers (mi, m 2 ■ ■ ■ , mk) 
satisyfying the equation 




This completes the proof of Theorem 2. 
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Abstract It is proved that product of interval number and triangular fuzzy number is a trape- 
zoidal fuzzy number. This paper transformed a multiple attribute decision making problem 
with attribute weight of interval number and attribute value of triangular fuzzy number into 
trapezoidal fuzzy number multiple attribute decision making problem, gave TOPSIS to trape- 
zoidal number multiple attribute decision making problem, and analyzed a practical problem 
by this method. 
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§1. Introduction 

The multiple decision making (MADM) problems mainly solve to rank alternatives with 
multiple attributes among finite alternatives. It is an extremely active topic in discipline re- 
search of modern decision science, systems engineering, management science and operations 
research, etc. Its theory and method are applied to many domains: engineering design, econ- 
omy, management and military affairs, etc. Reference [1] , [2] , [3] studied MADM problems that 
attribute weight information is interval number, attribute value is crisp number from various 
angles. Reference [4]- [7] studied MADM problems that attribute weight information is interval 
number, attribute value is also interval number from various angles. At present research results 
about methods for MADM problems with attribute weight of interval number and attribute 
value of triangular fuzzy number are few. This paper focuses on this kind of MADM problems 
mainly. This paper will prove that product of interval number and the triangular fuzzy number 
is a trapezoidal fuzzy number and transform a multiple attribute decision making problem with 
attribute weight of interval number and attribute value of and triangular fuzzy number into 
trapezoidal fuzzy number multiple attribute decision making problem, then gives TOPSIS to 
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trapezoidal number multiple attribute decision making problem. Finally, an example is used to 
illustrate the proposed method. 



§2. Presentation of the problem 



Definition 1. Let a be interval number, where a = [a~ ,a + ] = {x\a~ < x < a + ,a~ < a + }. 
If 0 < a~ < a + , then a is called a positive interval number. 

Definition 2. Fuzzy number p is called a triangular fuzzy number if its membership 
function pip{x) : R — > [0, 1] is described as 



Hp(x) = 



x l 



m — l m — V 
x u 



m — u m — u' 

0 , 



x £ [l, m] 
x £ [m, u\ 

otherwise 



where x £ R, l < m < u,l,u stand for the lower and upper bounds of the support of p 
respectively, and m for the most possible value of p. The triangular fuzzy number p can be 
denoted by {l, m, u). When l = m = u, p is a real number. 

Definition 3. Fuzzy number t is called a trapezoidal fuzzy number if its membership 
function jiji'x) : R , — >• [0, 1] is described as 



K( x ) 



0, 

x — a 
b — a ’ 
1, 

d — x 



d-< 



0, 



x < a 
a < x < b 
b < x < c 
c < x < d 
d < x 



where a < b < c < d 7 a and d stand for the lower and upper bounds of the support of t 
respectively, and closed interval [ b , c] for the most possible value of t. The trapezoidal fuzzy 
number t can be denoted by (a, b , c, d) . When a = b,c = d, t becomes a ordinary interval 
number. When b = c, t turns into triangular fuzzy number. Interval number and triangular 
fuzzy number are special cases of trapezoidal fuzzy number. When a = b = c = d, it is a 
ordinary real number. 

The A-cut of a trapezoidal fuzzy number t = (a, b , c, d) is 



t\ = [a + A (b — a),d+ A (c — d)], A £ (0, 1]. 

Definition 4. [8] Let vector T — (ti, t 2 5 ■ ■ * ,tm) and V = ,v m ), where trape- 
zoidal fuzzy number ti = (ti*, <2*, = {vu, v^i, v$i, Vo), i = 1,2, - ■■ ,m. The distance 

between T and V is defined as 



D(T,V) = \\T-V\\ 



4 N 



_ V li) 2 + (*2 i ^ V 2i ) 2 + {hi 

2 — 1 



V3i) 2 + {hi 
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We will use product of interval number and triangular fuzzy number in our discussion, for 
this reason we prove the theorem as follows. 

Theorem. Let a = [ 01 , 012 ] be arbitrary positive interval number, p = ( l,m,n ) be arbi- 
trary triangular fuzzy number, then the product of a and p is 

a ■ p = [cti, 02] • ( l , to, u) = (aiZ, aim, a^rn, o^u)- 

Namely the product of positive interval number and triangular fuzzy number is trapezoidal 
fuzzy number. 

Proof. For VA G (0, 1], the A-cut of p is 

p\ = [l + A (to — l),u + A (to — u)]. 

According to extension principle we have 

a-p = |^J A{[oi, 0:2] [Z + A (to — l), u + A(m — it)]} 

0<A<1 

= A[cti/ + A(oiTO — OiZ), Oilt + A(o2TO — Ct2w)] 

0< A<1 

= (aiZ, OxTO, Ct2TO,, 02w). 



§3. Principle and step 



This paper considers the fuzzy MADM problems that the alternatives are X -\ , A' 2 , •••, 
X n and attributes are Gi,G 2 ,--- ,G m . The fuzzy MADM problems with incomplete weight 
information, Wj denotes the weight of attribute Gj, Wj G [oy,/3j], where, 0 < ay < pj < 1, 
j = 1, 2, • • • ,m, w 1 + W 2 + • • • + w m = 1, is known to us. Triangular fuzzy number (Z^, rriij, Uij) 
is the value for alternative Xi with respect to attribute Gj, i = 1,2, •• • ,n, j = 1,2, - ,m. 

Principle and step for the fuzzy MADM problems that weight of attribute is interval number 
and value of attribute is triangular fuzzy number are as follows: 

Step 1. Construct triangular fuzzy number decision-making matrix A = (xij) nxm - Where 
triangular fuzzy number Xy = (lj j , , Uy) is the value for alternative Xj with respect to 

attribute Gj, i = 1, 2, • • • , n, j = 1, 2, • • • , to. 

Step 2. Normalize decision-making matrix A. 

Generally speaking, the attribute values have different physical dimension from each other, 
the decision matrix needs to be normalized so as to transform the various attribute values into 
comparable values. Therefore, we can obtain the normalized fuzzy decision matrix denoted by 
A = (Xy) nxm according to formulas in reference [9]. 

Let Xy = (ly, my, u^), then 



Xi- 



= (; 



rriy 
7/ max ’ 

a j 



^min 7 min ^min 

;)j J G Jl,Xy = , - 7 - , ~ ), j G J2 , 

Uij ! I lij Lij 



where u“ ax = maxit^-, if j G Ji, Z“ ln = min ly, if j G J 2 , where J\ and J 2 are the set of benefit 
J i J i 

criteria and cost criteria, respectively. 
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Step 3 . Construct the weighted normalized fuzzy decision matrix. 



Y = Wig) A = 



WlXll 


W 2 X 12 • 


Wm%lm 




2/u 


2/12 ’ 


Him 


W 1 X 21 


W 2 X 22 ■ 


' ^m^2m 


= 


2/21 


2/22 ' 


U2m 


WlX n i 


W 2 X n 2 ■ 


' 'U) rn %nm 




Vnl 


Vn2 ■ 


Unm 



where 



Uij — Wj • X%j 

\P^j ? Pj\ ’ (J ij 5 TTlij 5 'U'ij ) 

— (Qj lij , OLj TTlij , f 3 j Tflij , ( 3 j U-ij ) 

— (P'ij i ^ij 5 C“ij i d*ij ) 5 

i = 1, 2, — ,n,j = 1,2, - - - ,m. 

Step 4 . Determine the fuzzy positive-ideal solution X + and fuzzy negative-ideal solution 
X - , respectively, and 

^■ + = ,y™), x ~ = (2/f>% ,y ™) 



where 



yf = (max ay, max fey, max cy, max dy) = (a+, &+, c+, d+), j = 1,2,- 






2/,- = (.min ay, min by^min cy, min dy) = (a^ ,bj ,Cj ,dj ),j = 1,2,--- ,m. 



Step 5 . Calculate the distance of each alternative from fuzzy positive-ideal solution X + 
and fuzzy negative-ideal solution X - , respectively. 

The distance of each alternative from X + and X - can be calculated as 



Df = D(X i ,X+) = 



1 



^[(ay - a+) 2 + ( bij -b +) 2 + (cy - c+) 2 + (dy - d+) 2 ] , 
j=i 



D 7 =D(Xi,X-) = 



, “ a i ) 2 + “ b j ) 2 + ( Ci i - S' ) 2 + ~ d j ) 2 ]’ 

\ i=i 



£ = 1, 2, • ■ ■ , n. 

Step 6. Calculate the relative closeness coefficient of each alternative to the ideal solution. 



a = 



i): 



D~ + D 7 



,i = 1,2, •• • , n. 



Step 7 . Determine the ranking order of all alternatives according to the closeness coefficient 
in the descending order and select the best one. 
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§4. Numerical example 

This example is taken from [10], one of the Middle East country intends to purchase a cer- 
tain number jet fighter from US, The American Pentagon’s officials have provided information 
relative to 4 kinds of jet fighter types which are permitted to sell. The Middle East country 
sends out the expert group to carry on the detailed inspection to these 4 kinds of jet fighters, 
Inspection results are shown in Table 1. Which kind of jet fighter should they purchase to make 
total utility value of decision making biggest? 



Table 1: Expert group inspection result to 4 kinds of jet fighters 



Type 


Maximum 
speed (Ma) 


Cruising 
radius (mile) 


Peak load 
(pond) 


Price 

(/$ 10 6 ) 


Reliability 


Maintainability 


Xi 


2.0 


1500 


20000 


5.5 


High 


Very high 


x 2 


2.5 


2700 


18000 


6.5 


Low 


General 


x 3 


1.8 


2000 


21000 


4.5 


High 


High 


X 4 


2.2 


1800 


20000 


5.0 


General 


General 



weight Wj of attribute Gj cannot determine completely (j = 1,2, ■ ,6), but according to the 

expert analysis we know w\ £ [0.15,0.20], W 2 £ [0.1,0.15], W 3 £ [0.1,0.15], W 4 £ [0.1,0.15], 
w 5 £ [0.15,0.20], w 6 £ [0.2,0.25], 

Step 1. Correspond to normalized fuzzy decision-making matrix as follows: 

' (0.8, 0.8, 0.8) (0.56,0.56,0.56) (0.95,0.95,0.95) (0.82,0.82,0.82) (0.6, 0.8,1) 

(1,1,1) (1-bl) (0.86,0.86,0.86) (0.69,0.69,0.69) (0,0.2, 0.4) 

(0.72,0.72,0.72) (0.74,0.74,0.74) (1,1,1) (1,1,1) (0.6, 0.8,1) 

(0.88,0.88,0.88) (0.67,0.67,0.67) (0.95,0.95,0.95) (0.9, 0.9, 0.9) (0.3, 0.5, 0.7) 

Step 2. Construct the weighted normalized fuzzy decision matrix. 

(0.12,0.12,0.16,0.16) (0.056,0.056,0.084,0.084) 

(0.15, 0.15, 0.2, 0.2) (0.1, 0.1, 0.15, 0.15) 

(0.108, 0.108, 0.144, 0.144) (0.074, 0.074, 0.111, 0.111) 

(0.132,0.132,0.176,0.176) (0.067,0.067,0.1005,0.1005) 

(0.082,0.082,0.123,0.123) (0.09,0.12,0.16,0.20) (0.14,0.18,0.225,0.25) 

(0.069, 0.069, 0.1035, 0.1035) (0, 0.03, 0.04, 0.08) (0.06, 0.1, 0.125, 0.175) 

(0.1,0.1,0.15,0.15) (0.09,0.12,0.16,0.20) (0.12,0.16,0.2,0.25) 

(0.09,0.09,0.135,0.135) (0.045,0.075,0.1,0.14) (0.06,0.1,0.125,0.175) 

Step 3. Determine the fuzzy positive-ideal solution X + and fuzzy negative-ideal solution 
X~ , respectively, 

X+ = [(0.15, 0.15, 0.20, 0.20), (0.1, 0.1, 0.15, 0.15), (0.1, 0.1, 0.15, 0.15), 

(0.1, 0.1, 0.15, 0.15), (0.09, 0.12, 0.16, 0.20), (0.14, 0.18, 0.225, 0.25)] 

A'" = [(0.108, 0.108, 0.144, 0.144), (0.056, 0.056, 0.084, 0.084), (0.086, 0.086, 0.129, 0.129), 

(0.069, 0.069, 0.1035, 0.1035), (0, 0.03, 0.04, 0.08), (0.06, 0.1, 0.125, 0.175)] 



(0.095, 0.095, 0.1425, 0.1425) 
(0.086,0.086,0.129,0.129) 
(0.1,0.1,0.15,0.15) 
(0.095, 0.095, 0.1425, 0.1425) 



(0.7, 0.9,1) 
(0.3, 0.5, 0.7) 
( 0 . 6 , 0 . 8 , 1 ) 
(0.3, 0.5, 0.7) 



Step 4. Calculate the distance of each alternative from fuzzy positive-ideal solution and 
fuzzy negative-ideal solution, respectively, as Table 2. 
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Table 2: The distance measurement 





w 


x 2 


*3 


W 


D+ 


0.03522 


0.06953 


0.03124 


0.05555 


D~ 


0.06885 


0.03740 


0.06758 


0.03412 



Step 5. Calculate the relative closeness coefficient of each alternative to the ideal solution 
as 

Ci = 0.66157, C 2 = 0.34976, C 3 = 0.68387, C 4 = 0.38051. 

Step 6. According to the closeness coefficient in the descending order, the ranking order of 
four alternatives can be determined as follows: 

X 3 >- X-t X 2 . 

Therefore, the best selection is X 3 . 

The example result by method of this paper proposal is consistent with result by method 
of fuzzy compromise decision method which is first synthesis, then weighted decision making 
pattern introduced in reference [10]. Generally speaking, except some decision-making situa- 
tions with quite obvious differences, the result which the different method obtains isn’t always 
consistent. Therefore in the actual decision-making, we should use several different methods to 
solve, synthesize and compare result from each method in order to obtain a relatively reliable 
conclusion. 

§5. Conclusion 

This paper has proven the product of interval number and the triangular fuzzy number 
is a trapezoidal fuzzy number. We have proposed new method for multiple attribute decision 
making problem with interval number attribute weight and triangular fuzzy number attribute 
value. A numerical example is used to illustrate the method. Through the example, we can 
find this method is reliable, the thinking is clear and is one practical decision making method. 
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Abstract For any positive integer n, we define a new Smarandache function G(n) as the 

771 

smallest positive integer m such that P </>(fc) is divisible by n, where is the Euler 

k = 1 

function. The main purpose of this paper is using the elementary methods to study the 
elementary properties of G(n), and give three interesting formulas for it. 

Keywords A new Smarandache function, series, inequality. 



§1. Introduction and results 

For any positive integer n, we define a new Smarandache function G(n) as the smallest 

m 

positive integer nn such that <j>(k) is divisible by n. That is, 

k = 1 

m 

G(n) = mm{m : n | 4>{k), me N}, 

h — 1 

where 4>(n) is the Euler function, N denotes the set of all positive integers. For example, the 
first few values of G(n) are: G( 3) = 7, G(4) = 4, G(5) = 11, G(6) = 7, G( 7) = 29, G(8) = 5, 
G(9) = 9, G(10) = 11, G(ll) = 23, G(12) = 7, G(13) = 43, G(14) = 29, G(15) = 11, G(16) = 5 
• • • . About the properties of this function, it seems that none had studied it yet, at last we 
have not seen any related papers before. Recently, Professor Zhang Wenpeng asked us to study 
the arithmetical properties of G(n). The main purpose of this paper is using the elementary 
methods to study this problem, and prove the following three conclusions: 

Theorem 1. For any prime p, we have the calculating formulae 

G(p) = min{p 2 , q(p , 1)}; 

G(p 2 ) = q(p, 2), if q(p , 2) < p 2 ; G{p 2 ) = p 2 , if q(p , 1) < p 2 < q{p, 2); 

G{p 2 ) = q(p, 1), if P 2 < q(p, 1) < 2 p 2 ; and G{p 2 ) = 2 p 2 , if q(jp , 1) > 2 p 2 , 
where q(p, i) is the i-th prime in the arithmetical series {np+ 1}. 
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series 



Theorem 2. G(n) is a Smarandache multiplicative function, and moreover, the Dirichlet 
G(n) . 

y — is divergent. 



n= 1 



Theorem 3. Let k > 2 be a fixed positive integer, then for any positive integer group 
(mi, m2, • • • , nik ) , we have the inequality 

G(mim2 • • • mk ) < G{m\)G{m 2 ) • • • G(mfc). 



Note. For any positive integer n, we found that n = 1, 4, 9 are three positive integer 
solutions of the equation G(n) = n. Whether there exist infinite positive integers n such that 
the equation G(n) = n is an interesting problem. We conjecture that the equation G(n) = n 
has only three positive integer solutions n = 1, 4, 9. This is an unsolved problem. 



§2. Proof of the theorem 



In this section, we shall complete the proof of our theorems directely . First we prove Theo- 

m s 

rem 1. Let G(p) = m, namely p | (f>{k), p \ ]^[ 0(fc) , 0 < s < to. If in = q^q^ 2 • • • q? r be the 

k= 1 fc=l 

factorization of m into prime powers, then (fi(rn') = l l (qi — 1 ) 92 2_1 (92 — 1) • • • qr r l (qr — 1), 
from the definition of G(n) we can deduce that p divides one of 1 < i < r. So p divides 

q or q t — 1. We discuss it in the following two cases: 

(i) . If p | qi~ 1, then we must have m = qi = lp+1 be the smallest prime in the arithmetical 
series {kp + 1}. 

(ii) . If p | then we have have on = 2, m = qf = p 2 . 

Combining (i) and (ii) we may immediately deduce that G(p) = min{p 2 , Ip + 1}, where 
Ip + 1 be the smallest prime in the arithmetical series {kp+ 1}. Similarly, we can also deduce 
the calculating formulae for G(p 2 ). This proves Theorem 1. 

Now we prove Theorem 2. For any positive integer n > 1, let n = p^p^ 2 • ■ ’P“ r be the 
factorization of n into prime powers, G(p“ ; ) = m*, 1 < i < r, m = max{mi, m- 2 , • • • , m r }, 
then from the definition of G(n) we have divides ^(l)</>(2) • • • ^(m,) for all 1 < i < r. So pf ' 
divides (f>(l)<j>(2) ■ ■ ■ <j>(m). Since ( Pi , pj) = 1, i ^ j, so we must have n = p^p^ 2 • • •p“ 1 ' divides 
Therefore, G(n) = to = maxlG^" 1 ), G(p 2 2 ), ■ ■ • , G(p“’ )}. So G(n) is a 
Smarandache multiplicative function. From Theorem 1 we may immediately get 



y- G(n) y- G(p) y^ 1 _ 

/ ^ n 2 / -J 7)2 / v 7) 



n — 1 



So the Dirichlet series 



G(n) . „ 

— is divergent. 
n- 



Finally, we prove Theorem 3. First we prove that the inequality G(toiTO 2 ) < G(toi)G(to 2 ) 
holds for any positive integer m\ and TO 2 . 

Let G(to i) = it, G(m 2 ) = u, from the definition of G(n) we can easily get 



TOi | n </>(*) , rn 2 | n# 

2=1 2=1 



28 



Jing Fu and Yu Wang 



No. 3 



Without loss of generality we can assume u < v, then 

uv u uv u 

nw)=n*(*)- n 

k = 1 k = 1 k=u-\-l k—1 

Notice that 



</>( 2) | 0(2i>), 0(3) | 0(3w), ••• , 0(u) | 0 (mu), 0(w + 1) | 0(u+l), , (t>{v) | 0(u), 

this means that 

V UV 

i n 

2—1 fc=u+l 

or 

U V uv 

n^on^o i n ^ ■ 

2 = 1 2=1 k—1 

Hence 

UV 

777-1777-2 | n 

k= 1 

From the definition of G{n ) we know that G{m\m, 2 ) < uu, or 

G(toiTO2) < G(mi)G(rri2). 

If k > 3, then applying the above conclusion we can deduce that 

G(miTO2 • • • TOfe) = G(m\[rn.2 ■ ■ ■ m^)) < G(m\)G(rn2 ■ ■ ■ mu) 

< G(TOi)G(m 2 )G(TO3 • • • TO*;) 



< G(toi)G(to 2 ) • • • G(mfc). 

This completes the proof of Theorem 3. 
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Abstract In this paper, we use the elementary methods to study some divisibility of the 
Smarandache numerical carpet, and obtain two interesting congruent theorems. 

Keywords Divisibility, Smarandache numerical carpet, elementary method. 

§1. Introduction and results 

In his book “Only problems, not solutions”, F. Smarandache introduced many interesting 
sequences and arithmetical functions, at the same time, he also proposed 105 unsolved problems 
and conjectures. One of the problems is that he asked us to study the properties of the numerical 
carpet. It has the general form 

1 

1 a 1 
1 a b a 1 
1 a b c b a 1 
labcdcbal 
labcdedcbal 
labcdcbal 
1 a b c b a 1 
1 a b a 1 
1 a 1 
1 

On the border of level 0, the elements are equal to“l”, they form rhombus. 

Next, on the border of level 1, the elements are equal to “a”, where “a” is the sum of all 
elements of the previous border, the “a” s form a rhombus too inside the previous one. 

Next again, on the border of level 2, the elements are equal to “6”, where “5” is the sum 
of all elements of the previous border, the “6”s form a rhombus too inside the previous one. 
And so on • • • 
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It can be described as 



1 



1 



1 a\ . 
1 di d2 . 

ai a 2 a 3 . 

1 di 02 . 
1 Oi . 



Cll 


1 




02 


Ol 


1 


«3 


«2 


Ol 1 


02 


«1 


1 


Ol 


1 





1 

Let n , k be integers with n > 0 and k > 0, the n — th square can be presented as a set of 
sub-squares. In [2], Mladen and T. Krassimir gave and proved the general formula: 



k—2 

a k = 4 (k + 1) JJ(4n + 1 - 4 i), 
i = o 

where integer k > 2 . 

Let C(?z, fc) denotes a* in n-th square, then 



C(n, k ) = 1, for k = 0; 

C(n, fc) = 4(n + 1), for k = 1; 

k—2 

and C(n, k ) = 4(n + 1) (4 n + 1 — Ai), for 2 < k < n + 1. 

i— 0 

The main purpose of this paper is using the elementary methods to study the divisibility 
properties of the sequence {C{n, fc)}, and give two interesting congruent theorems. That is, 
we shall prove the following conclusions: 

Theorem 1. Let n be a positive integer. Then for any prime p with 4n— Ak+9 < p < 4n+l 
and p = 1 mod 4, we have the congruence 



C(n, k) = 0 mod p. 



p < n + 



Theorem 2. Let n be a positive integer, p be any prime with n — k + 3 - 

n + 1 



n — k + 2 



< 



and p = 3 mod 4, then we have the congruence 



C(n, k) = 0 mod p. 
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§2. Proof of the theorems 

In this section, we shall use the elementary methods to complete the proof of the theorems. 
First we prove Theorem 1. For any positive integer n and prime p = 1 mod 4, if An — Ak + 9 < 
p < An + 1, then it is clear that An — Ak + 9 < 4n + 1 or k > 2, from this we have 



k—2 

C(n, k ) = 4(n + 1) (4n — Ai + 1) 

i=0 

k—2 

= 4(n + 1) ]^[(4(n - i) + 1). 

i—0 

It is obviously that C (n, 2) = 4(n + 1) (4n + 1) , so if k = 2, then from An — Ak + 9 < p < An + 1 
we know that p = An + 1, and C(n , 2) = 0 mod p. This means that Theorem 1 is true. 

Now we assume that 2 < k < n, so that 



k—2 



k—2 



C(n, k) = 4 (n + 1) J^[(4(n — «) + !)= (n + 1) • A k ■ (n — i + 4) mod 



i = o 



i= o 



k—2 



= (n + 1) • 4 fe • ( n - i 



■i = o 



3p + l 



mod p. 



(1) 



where x denotes the solution of the congruence equation x ■ x = 1 mod p, and 4 = — . 

Since An ~ Ak + 9 < p < An + 1, so there must exist an integer 0 < i < k — 2 such that 
p = n — i H — — . Therefore, from (1) we may immediately deduce that 



k—2 

C(n, k ) = 4(n + 1) (4(n — i) + 1) = 0 mod p. 

i= o 



This proves Theorem 1. 

Now we prove Theorem 2. Let n be a positive integer, p 



n — k + 3 + 



n — k + 2 



< p < n 



n + 1 



be a prime with 
and p = 3 mod 4. 



Note that 4 = 



p+ 1 

4 



, then from the definition of C(n, fc) we have 



k—2 



k—2 



C(n, k) = 4(n + 1) ]^[ (4(n — *) + !)= (n + 1)4 (n — * + 4) mod 



i=0 



i=0 



= (n + 1) 



p+ 1 



k—2 k—2 



n 

i= 0 



P+1 



mod p. 



(2) 



P -|“ 1 P 1 

From congruence (2) we know that if n — fc+2H — — < p < n-\ — , then C(n, k ) = 0 mod p. 

p + 1 P + 1 . 

The inequality n — k + 2 H < p < n H implies 



n — k n+ 1 

n — k + 3 H — <p<nl . 



3 



3 
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Since prime p must be a positive integer, so we have 

n — k + 2 



n — k + 3 + 



< p < n + 



This completes the proof of Theorem 2. 



n + 1 
3 
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Abstract For any positive integer n > 2, the Smarandache Inferior Prime Part {p P (n)} 
is defined as the largest prime less than or equal to n; The Smarandache Superior Prime 
Part (Pp(n)} is defined as the smallest prime greater than or equal to n. In this paper, we 
defined some determinants involving the Smarandache prime part sequences, and introduced 
two conjectures proposed by professor Zhang Wenpeng. 
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§1. Introduction 



For any positive integer n > 2, the Smarandache Inferior prime part {p p (n)} is defined as 
the largest prime less than or equal to n. The Smarandache Superior prime part {P p (n)} as the 
smallest prime greater than or equal to n. For example, the first few values of these sequences 
are p p ( 2) = 2, p p ( 3) = 3, p p ( 4) = 3, p p ( 5) = 5, p p ( 6 ) = 5, p p { 7) = 7, p p ( 8) = 7, p p ( 9) = 7, 

Pp(10) = 7, p p { 11) = 11, Pp { 12) = 11, Pp { 13) = 13, p p (14) = 13, Pp ( 15) = 13, Pp (16) = 13, 

Pp( 17) = 17, Pp ( 18) = 17, p„(19) = 19, 

P p (l) = 2, P p ( 2) = 2, P p (3) = 3, P p (4) = 5, P p ( 5) = 5, P p (6) = 7, P p ( 7) = 7, P p ( 8) = 11, 

P p ( 9) = 11, P p (10) = 11, Pp(ll) = 11, Pp( 12) = 13, Pp( 13) = 13, Pp(14) = 17, P p (15) = 17, 

P p (16) = 17, Pp( 17) = 17, P p ( 18) = 19, Pp(19) = 19, 

From the definitions of these sequences we know that for any prime q , we have p p {q) = 

P p (q) = Q- 

In his books “Only problem, Not solutions” W and “Sequences of Numbers Involved in 
Unsolved Problem” I 2 1, Professor F. Smarandache asked us to study the properties of these se- 
quences. About these two problems, some authors had studied them, and obtained some 
interesting results, see references [3], [4], [5] and [6]. For example, Xiaoxia Yan [6] studied the 



the asymptotic properties of 




and proved that for any positive integer n > 1, we have the 



asymptotic formula 



Sn 

In 



1 + 0(n 3 ), 



where I n = {p p ( 2) + p p { 3) H \-p p (n)}/n and S n = {P p (2) + P p (3) H b P p {n)}/n. 



1 This work is supported by the Shaanxi Provincial Education Department Foundation 08JK433. 
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The sequences {p p (n)} and {P p (n)} are very interesting and important, because there 
are some close relationships between the Smarandache prime part sequences and the prime 
distribution problem. Now we introduced two determinants formed by the Smarandache prime 
part sequences. They are defined as follows: For any positive integer n, c(n), C(n ) are n x n 
determinants, namely 



and 



Pp( 2) 


Pp( 3) 


P P {n+ 1) 


p p (n + 2) 


P p (n + 3) 


p p (2n + 1) 


p p {n(n - 1) + 2) 


p p (n(n — 1) + 3) •• 


p p (n 2 + 1) 


^p(I) 


P P ( 2) 


••• p p{n) 


P p {n+ 1) 


Pp{n + 3) 


■■■ -P p (2n) 


P p (n(n - 1) + 1] 


P p {n{n - 1) + 2) 


••• P P (n 2 ) 



For example, by definitions and calculating, we can find some values of these determinants 
as the following table: 



n 


c(n) 


C(n) 


n 


c(n) 


C(n) 


2 


1 


4 


3 


14 


4 


4 


0 


188 


5 


-96 


-1424 


6 


0 


0 


7 


1.0214 x 10 5 


37536 


8 


0 


0 


9 


0 


0 


11 


7.8299 x 10 7 


1.1478 x 10 s 


13 


9.8338 x 10 8 


-2.7958 x 10 9 


17 


8.2462 x 10 14 


1.3164 x 10 13 


19 


-1.9608 x 10 15 


1.629 x 10 15 


23 


2.4545 x 10 20 


3.156 x 10 18 


29 


8.9308 x 10 27 


-6.5008 x 10 27 


31 


1.1095 x 10 29 


-2.2835 x 10 28 


37 


9.9497 x 10 37 


-7.4692 x 10 38 


41 


-1.1502 x 10 42 


3.9244 x 10 45 


43 


1.5152 x 10 45 


-2.1671 x 10 44 


47 


1.1606 x 10 51 


4.06 x 10 5 ° 


53 


2.9359 x 10 59 


5.8735 x 10 59 


59 


3.8402 x 10 67 


-3.1043 x 10 69 


61 


-3.614 x 10 7 ° 


6.9858 x 10 72 


67 


3.2341 x 10 8 ° 


-9.5374 x 10 78 


71 


-2.219 x 10 85 


4.7688 x 10 86 


73 


-3.6656 x 10 85 


-3.6985 x 10 89 


79 


-4.2038 x 10 98 


-3.8762 x 10 97 


83 


1.6966 x 10 104 


2.1389 x 10 104 


89 


-1.0695 x 10 113 


5.7824 x 10 110 


97 


2.8.98 x 10 124 


1.9968 x 10 124 









About the elementary properties of these determinants, it seems that none had studied 
them, at least we haven’t seen any related papers before. Recently, Professor Zhang Wenpeng 
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asked us to study the properties of these determinants, at the same time, he also proposed 
following two conjectures: 

Conjecture 1. For any composite number n > 6, we have the identities c(n) = 0 and 
C(n) = 0; 

Conjecture 2. For any prime q , we have c(q) ^ 0 and C(q) ^ 0. 

From the above table, we believe that these two conjectures are correct. About Conjecture 
1, we have solved it completely, which will be published in Pure and Applied Mathematics. 
But for Conjecture 2, it’s still an open problem. We think that Conjecture 2 is very interesting 
and important, since if it is true, then we can get a new discriminant method to distinguish 
prime from integer numbers through calculating these kinds of determinants. So we introduce 
the conjectures in this place, and hope more scholars who are interested in to study it with us. 
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Abstract The main purpose of this paper is using the elementary methods to study the rela- 
tionship between the continued fractions and Dirichlet L-functions , and give some interesting 
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§1. Introduction 

Pn 

Let xq, Xi, X 2 , ■ ■ ■ be integers with Xj > 1, j > 1. We can express the original fraction — - 

Qn 

in the form xq + — L_ . n > 0. Such an expression is called a continued fraction, and 

‘i' OeS .V 

' x n 

we denote it by (xo, X\, ■ • ■ , x n ). 

About the characteristic properties of continued fractions, it can be found in the references 
[1] and [2]. In this paper, we shall study the relationship between the Dirichlet L-functions 
and the continued fractions, and obtain some interesting identities. That is, we shall prove the 
following main conclusion: 

Theorem. Let p be a prime, a be a positive integer with (a, p) = 1 and 1 < a < p. Then 
we have the identity 



x(°) 1-^(1) x)l 2 

X mod p,x(~ 1) = — 1 

' n o 

^(-1)*+^ + + - + -(-! + (-1)") , 

, n Qn P * 

.k — 0 J 

where x be the Dirichlet character mod p, L(l,x) be the Dirichlet L-function corresponding 
character % and a/p = p n /<ln = (a^O) x h ■ - • , x n ). 

p — 1 

For a = 2, 3,4 and , from the above Theorem and the definition of continued fractions 

2 

we may immediately obtain the following four Corollaries. 

1 This work is supported by the Shaanxi Provincial Education Department Foundation (08JK398) of 

P.R. China 



7T 2 (P- 1) 

12 p 
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Corollary 1. Let k > 1 be an odd prime, then we have the identity 



E 

X mod k,x(— 1)=— 1 



x(2)|£(i,x)l 2 = g^ 



1 3 

2 ~ k 



Corollary 2. Let p be an odd prime and \ be the Dirichlet character mod p. Then 

7T 2 (p-l) 2 (p-l0) 



E X(3)|£(1,X)I ^ ( p _ lj(p - 2)(p - 5) 



\ mod p,x(-l)=”l 



if p = 1 mod 3; 
if p = 2 mod 3. 



36 p 2 

Corollary 3. Let p be an odd prime and x be the Dirichlet character mod p. Then 



E X(4)|L(1, X )| 2 = 

X mod p,x(-l)=-l 



7T 2 _ (p- l) 2 (p- 17) 

48 p 2 

7 r 2 (p— l)(p 2 - 6p + 17) 



if p = 1 mod 4; 
, if p = 3 mod 4. 



48 p 2 

Corollary 4. Let p be an odd prime and x be the Dirichlet character mod p. Then 



E x 

X mod p,x(-l) = -l 



P~ 1 



\ l ( i , x )\ 2 = - 



7 r 2 (p- l)(p 2 -4p + 5) 

2V ' 



It is clear that these four Corollaries are an extension of Walum [3] . 



§2. Some Lemmas 

To complete the proof of the Theorem, we need the following three Lemmas. 

Lemma 1. Let Xq,Xi,--- ,x n be integers and (0 < i < n) be natural numbers, 

continued fraction (xo,Xi, ■ ■ • , Xj) = — . Then 

Qi 

PiQi-i -Pi-iQi = (-1)* +1 , -1 < i < n. 



and P —2 = 0,p_i = l,g -2 = 1, 9—1 = 0,pi = XiPi-i +p l _ 2 ,qi = Xiqi-i + 9i- 2,0 < i < n. 
Proof. (See reference [1] ). 

Lemma 2. Let k and h be integers with k > 3 and (h, k) = 1. Then we have 

= S x{h)\m, X )W 

d\k x mod d,x{~ 1)=— 1 



where S(h, q ) = 



a=l 



a/i\ 

q ) 



I I is the Dedekind sum and <f>(k) is the Euler’s function. 



Proof. (See reference [4]). 

Lemma 3. Let a and q be two positive integers with (a, q) = 1. Then we have the identity 



S{a,q) 



1 

12 



]T(-i) fe+1 x fc + (-i) 



a 

q 



-1 + (- 1 )" 

8 
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where a/q = p n /q n = (x 0 ,xi, ■ ■ ■ ,x n ). 

Proof. Let continued fraction (xo,Xi,--- ,x n ) = — and suppose — = — . For the 

q n q n q 

successive terms — and — — from Lemma 1 we know that 
Qn Qn— 1 



(1) 

(2) 



Pn9n-1 - 9nPn-l = (-1)" +1 , Tl > -1. 

Using the properties of Dedekind sums and ( 1) we get 

S(p n ,q n ) = S{p n q n _iq—[, q n ) 

= S (<^TT((-1)” +1 +Pn-iqn) i q n ) 

= S{(-l) n+1 q^,q n ) 

= (-1 ) n+1 S(q^,q n ) 

= (-1 ) n+1 S(q n ^,q n ), 

where q n -i denotes the solution x of the congruence equation xq n — i = 1( mod q n ). From the 
reciprocity formula and Lemma 1 we obtain 

qn + q ^.- i + 1 1 



1 ) Qn) “k 9n— l) — 



( q n - 1 , q n 

12 
l 

12 



( 3 ) 



+ 



qn qn— 1 
qn— 1 qn— 2 



1 1 



qn qn— 1 

Hence, by expression (3) and Lemma 1 we obtain 
S(qn—1 • qn) 

, . 1 / Ori-1 . 9n-2 . 

h X 



_|_ /_]\n+l ( ElL — Pn-1 

q n q n - 1 



1 

4 



( 4 ) 



— S{c[ni qn— l) ~k 12 
= —S(q n - 2 , q n -i) + ^ 
= S{q n - 3 ,q n - 2 ) - ^ 



Qn Qn— 1 
(Zn-l , <?n-2 



9n <?n— 1 
Qn—2 . Qn—S 



_|_ j _ ^"-1 

q n Qn - 1 



_!_ (_^n+l ( ElL — Pn-1 



Qn— 1 Qn—2 



+ X n — i + ( — l) r 



1 

12 



Qn—1 Qn—2 



Qn 



q n - 1 

= S{q n —3iQn—2) H - Y 2 



_|_ ^_^n+l | Pji _ Pn—1 



Qn Qn—1 



Qn Qn— 1 
Pn—1 Pn—2 
Qn—1 Qn—2 

1 

4 



1 

4 

1 

■4 



<7n-l 9n-3 



qn qn—2 

If n is an even number, then from (4) we have 

1 



x„-x n _ 1 ) + (-ir+ i (^-^ 

9n qn—2 



S(q n -i,q n ) = s(q-i,qo) + 



1 

12 

1 

12 



12 

^ ( 1 ) X n —i + 

i—0 

E(-i • 9n_1 



+ (-ir +1 f ^ ^ + vVi) 4 ^ 

V 90 / V9n <W ^ 



5 ) 



n— 1 



+ (-i) n+1 ^ + (-i) n ^ 

q n q n qo 



Li=0 



Qn 



+ (-1) 



i+l Pn 
Qn 
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If n is an odd number, then from (4) we have 



S(q n -i,q n ) = S(q 0 ,qi) + — 



+ — 



q n - 1 qo 



9n qij \q n qij 



n — 2 



(6) 



qn - 1 _ J_\ + fPn _ ffpXi + 1 

< 7 n Xi/ \ 9 ™ *1 



n — 2 






2—0 



(®i - l)(a;i - 2) 1 

+ 12 



1 

12 



12a’i 

n 

EC-O^n-i 



o iizf 

9n— 1 I / i\n+lPn A Y'V -.xi 

T ( 1) 2^0 T / ( 1) X n —i 

qn q n Xi ^ 

2 — 0 



, 2=0 



qn-1 /_Y'j n + 1 — 

q n q n 



Now combining (2), (5) and (6) we obtain 

1 



S(Pn,q n ) = 



12 



E(-i) fc+ V + (-l)™- 1 ^ + Pn 



k—0 



Note that p n /q n = a/q we may immediately get 



SO, 5) = ^ 



E(-i) fe+1 ^ + (-i) 



,fc— 0 



Qn Qn 



i—l Qn —1 a 
Qn Q 



-i + (-i) n 



-i + (-i) r 



This completes the proof of Lemma 3. 



§3. Proof of the theorem 



In this section, we shall complete the proof of the theorem. Using Lemma 2 we have 



E 

x mod p,x(-l) = -l 



x{a)\L{l,x)\ 2 = — — S(a,p). 



(7) 



Then from Lemma 3 and expression (7) we can easily get 

E x(a)\L(l,x)\ 2 

X mod p,x(-l)=-l 



\P~ 1) 



K 2 {p- 1) 
12 p 



£(-!)»+>« Hd - l )"-' 9 "- 1 ' ■ - 1 + <- 1 )' 



\k = 0 



Qn V , 



E(-i) fc+1 ^ + (-i) 



,k—0 



^qn-i + a + 3(_ 1 + (_ 1)n ) 
qn P 2 



This completes the proof of Theorem. 

The proof of the Corollaries. 

For any integer k > 1 and ( h , k) = 1, using the Mobius inversion formula and Lemma 2 we 



have 



X mod k,x{~ 1)= — 1 d\k 



(8) 
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Taking h = 2 in (8), note that 5(2, 1) = 0 and the continued fraction 

2 P2 /n k/d - 1 

777 = — = (0, 7 ,2) = 

fe/a q 2 2 

If k be a odd prime, then from (8) and Theorem we have 



E 

X mod fc,x(— 1)=— 1 



X(2)|L(1,X)| 2 = E ^ S <2. 



= ^2^0 ^ 2 ^ 



fc 2 -— ' d 

d\k,d<k 
ir 2 (f>(k ) y- /i(d) 

12fc ^ ~cT 

d|fc,d</c 



^(-l) fc+1 a; fc + (-!)-+ 2 
*=£ 92 P. 



7 T 2 (j){k) 

"12 . ~cT l 2 ^ ~~k/d 'k/d) 



E 

d|fc,d</c 



/i(d) (k/d— 1 k/d— 1 2 \ 



7T 2 0 2 (fc) 

12 k 



1 3 

2 ~~ fc 



This proves Corollary 1. 

Using Theorem and the definition of continued fraction we can also deduce that 



tt 2 . ( p - l) 2 (p- 10) 

E X(3)|i(l,x)| 2 = ^ p (p_ i)( > p_2)( J 5-5) 

X mod p,x(-l)=~l 
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E X(4)|L(1, X )| 2 = 

X mod p,x(-l)=-l 



7T 2 . (p- l) 2 (p- 17) 

48 p 2 

7T 2 (p— l)(p 2 - 6p+ 17) 

48 p 2 



and 



if p = 1 mod 3; 
if p = 2 mod 3. 

if p = 1 mod 4; 
if p = 3 mod 4. 



E x 

X mod p,x(-l)=-l 



p- ,r/i -.m2 7T 2 (p- l)(p 2 -4p + 5) 



|i(i,x)r = - 



24p 2 



This completes the proof of Corollary 2, 3 and 4. 



References 

[1] Pan Chengdong and Pan Chengbiao, Elementary Number Theory, Beijing University, 
Beijing, 1992, 305-315. 

[2] G. H. Hardy and E.M. Wright, An Introduction to the Theory of Numbers, Clarendon, 
Oxford, 1960. 

[3] H. Walum, An exact formula for an average of L-series Illinois J. Math, 26(1982), 1-3. 

[4] Wenpeng Zhang, On the mean values of Dedekind sums, Journal De Theorie des Nom- 
bres de Bordeaux, 8(1996), 429-442. 



Scientia Magna 
Vol. 4 (2008), No. 3, 41-46 



Research on the connectivity of public 
transportation network based on set method 

Chunming Liu*, Wansheng He*, Hongming Xia* and Zhongshe Gao* 

| Department of Computer Science and Technology, Tianshui Normal University 
Tianshui, Gansu, 741001, P. R. China 
| School of Mathematics and Statistics, Tianshui Normal University 
Tianshui, Gansu, 741001, P. R. China 
E-mail:cucmeLiu@gmail.com 

Abstract A connectivity judgment method for public transportation network based on set 
idea is studied in this article. Set is successfully applied to solve the problem of connectivity 
judgment on multi- vertex and multi-edge situation. Search level control strategy is adapted 
to control the search level of route, and three buffer tables are designed to improve the search 
speed and reduce the complexity of the algorithm. 

Keywords Traffic management, connectivity, transportation network, search level control. 



§1. Introduction 

Public transport is the main travel way for urban civilians and it is also a way advocated 
by many governments and organizations under the environmental protection concept. On one 
hand, along with the rapid construction of urban traffic, public travel becomes more and more 
convenience; on the other hand, problems come forth with multi-route selection. The civilians 
had to balance among fee, time cost, distance and some other areas, especially those in large 
cities or metropolis. And all this must be built on the basis of the pathway connectivity 
judgment, which are the major problem, as well as the common pathway connectivity judgment 
and optimization strategy under the circumstances of huge number of bus stations and bus lines, 
discussed in this article. 



§2. Some known results 

2.1. Problem description 

The pathway connectivity discussed in this article means the connectivity between one 
station and another, directly or through a third or more other stations, then find out the routes 
and the shortest one if they are interlinked. 

Basically, there are three types of public transport network distribution which include, all 
stations on the up- run line and down-run line are the same (also called Bidirectional Lines), 




42 



Chunming Liu, Wansheng He, Hongming Xia and Zhongshe Gao 



No. 3 



part of the stations are the same (also called Unilines) and Loop Lines. Meanwhile there are 
two types of fee system known as Flat Fare and Sectional Fare. Following are three lines’ 
information of a city’s public transport system. 

?135, Bidirectional Line; Sectional Fare. 

S 0417 — S 0272 — S3425 — 83476 — S2974 — 80234 — $0521 — «3806 — Si 682 — S3559 — $0928 _ S 3564 ~ 
S 0079 — 83175 — S2866 — S 3172 ~ 5 3269 — S 3603 ~ Sl 633 ~ S 2578 — 82579 ~ S2576 — 8 2577 

?i8i, Uniline; Flat Fare. 

Up-Run Line: 

-80238 “ 80540 — 80542 — S3024 _ S 0494 — S 0973 “ S3077 _ S 2082 — S 2213 ~ S22IO — S3332 — S3351 — 
S 1414 — S 0228 — S 0916 — S28II — Sooi 9 — S0966 

Down-Run Line: 

S 0966 “ 80019 — 828II — S0916 — So228 — S1414 — S3351 — S3332 — S2210 ~ S 2213 — 8 2082 _ S3077 — 
S 0973 — 80494 — S3024 ^ S0542 ~ 80541 — S0152 ~ S1765 — S3506 — S0238 

Z590, Loop Line; Flat Fare. 

S3295 — S3265 — 81334 — S1333 — S 2285 — S1341 — 80327 — 80317 — 80325 — S 0532 — S1734 — S 0314 — 
S1749 — 81831 — S1337 — S 0705 — 80704 — S3295 

Where l n is the n-th line’s No. of the city and s m is the m-th station of the city. 

2.2. Specialties of the problem 

Algorithm for connectivity of graph (include Breadth First Search and Depth First Search), 
adjacency matrix and genetic algorithm are all popular methods adopted to judge pathway con- 
nectivity [ 1 ]. Outwardly, public transport network can be convert to a digraph whose vertexes 
are stations and edges are lines between two stations. However there are actually more than one 
line that go through each station in the transport network, namely, there are unfixed number 
of edges between two adjoining vertexes. Therefore, the commonly used router algorithms for 
connectivity of graph are inapplicable since the time complexity and space complexity increase 
rapidly, for instance, there are more than 850 bus lines and more than 6400 bus stations in city 
Beijing [ 2 ], since the transport network is a n-th Iterated Line Digraphs [ 3 ]. 

Moreover, on one hand, because of the characteristics of intersecting and overlapping, there 
are always more than one line between two adjoining stations. On the other hand, even if there 
are route between two stations, it may need ongoing transferring, however, it is contrary to the 
public’s travel habits which are the less transferring, the better. According to the experience 
of livelihood, the public trend to reduce their transferring frequency when other factors’ (such 
as cost, time, etc.) difference is not very significant. That is, if there are more than one route 
between two stations, direct line has priority over transferring and transferring once has priority 
over transferring more times and so on, and other factors will be taken into considering only 
when transferring frequencies are equal. Therefore, the following discussing was based on the 
public’s travel habits. 



§3. Set method for connectivity 



3.1. Set description 
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Provide elements It , Sj and sets L t , Sj, Aj, which are 

If Element that indicate line No. i; 

Sj : Element that indicate station No. j; 

Lf. Set of all lines that pass through station sf 
Sj : Set of stations that line lj pass through; 

Aj : Set of stations that can be reached start from Sj ; 

B(j'. j) : Set of all stations that can be reached with line l- t start from station Sj ; 

The following expressions can be drown out: 

1). Direct Line: the following expression is obtained if there have direct line run from start 
station s start to destination station Sdest- 



S start * Sdest '■ ' ^ }s start, Sdest} — A s t ar t H Ablest • 



( 1 ) 



Where s start — ► Sdest indicate there are direct line run from s start to Sdest- 
Therefore, when {s start, Sdest} Q A s t a rt D Ad es t is establish, namely, there are direct line 
run from s s t. a rt to Sdest, the collection C of all direct lines are the intersection of set L start 
whose elements are lines that run through s s tart and set Ld es t whose elements are lines that 
run through s de st- 



C — L start bl Ldest • (2) 

2). Transfer Once: The circumstance of transfer once from s start to Sdest can be converted 
to the iterate of two direct lines, that is, the direct line from the start station s start, to the transfer 
station Stransfer and the direct line from the transfer station Stransfer to the destination station 
Sdest ■ Therefore the set D of transfer stations is: 



D — A s tart 0 Adest • ( 3 ) 

3). Transfer multi-times: Similarly, The circumstance of transfer n-times from s start to 
Sdest can be converted to the iterate of ( n + 1) direct lines. 

Direct line and transfer- needed circumstances are schematically shown in Fig.l. 




as Uiteof l.iiw bj Transfer Onao c) Transfer 'SViw 



Figure 1: The schematic view of (a) direct line, (b) transfer once and (c) transfer twice of 
two stations 



3.2. Optimization of algorithm 
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As has been analyzed that the public always trend to choose routes with less transferring 
frequencies, therefore, direct lines should be searched in top priority; deepen search level when 
direct lines are not available, and so on and so forth, until searching accomplished (succeed 
or failed) [4]. In addition, three memory tables known as line-station hash table, line-station 
searching linked-list and station-line searching linked-list are designed to improve the searching 
speed and reduce the searching complexity. 

3.2.1. Route searching strategy 

Search- level control strategy was introduce to integrate the Breadth-First Search method 
(BFS) and Depth-First Search method (DFS). The concrete steps are list in the following: 

1) . Construct set L s t ar t, Ldest i A start and A dest for s start and Sdest- 

2) . Figure out the intersection of L sta rt and Ldest , L start ALd e st- Then find out all identical 
elements in both set L start and set Ld es t via BFS, thus the result C is the set of all direct lines. 

a) . If C is not empty, record all elements in C which are all direct lines. 

b) . If C is empty, deepen the search levels and continue with step 3); 

3) . Figure out the intersection of A sta rt and Ad es t > A sta rt fl Ad es t- Then find out all identical 
elements in both set A sta rt and set Ad es t via DFS, thus the result set D is the set of all transfer 
stations. 

a) . If I? is not empty, record all elements in D which are all transfer stations. Then search 
direct lines with s s tart and D, Sdest and D , respectively. 

b) . If I? is empty, keep on deepening the search levels. 

3.2.2. Hash table structure of line-station 

Hash table structure of line-station stores the relationship of lines and stations. When the 
hash value, which are corresponded a line and a station, is “1” , the corresponding line runs 



through the corresponding station, 
is schematically shown in Fig. 2. 


and 


“0” 


, doesn’t. 


The hash table structure of line-station 




Sl 


S2 


S3 


Sn— 1 


Sn 


h 


0 


1 


0 ••• 


1 


0 


h 


1 


0 


0 ••• 


1 


0 


h 


0 


0 


1 ••• 


0 


1 


lm — 1 


0 


1 


0 ••• 


1 


1 




0 


0 


1 ••• 


0 


1 



Figure 2: Schematic view of the hash table structure of line-station 

As can be clearly seen from Fig. 2 that row direction of the hash table represents the lines 
and column direction, the stations. Moreover, the hash value “1” in row direction marks all 
stations that the corresponding line runs through, and in column direction marks all lines that 
run through the corresponding station. As is known to all that the time complexity of search 
algorithm using hash table structure is the constant of 0(1) [5]. 
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Two important logic relationships are implicated in the above hast table structure, known 
as: 

1) . There are direct line between every two stations whose hash value is “1” in the same 

row. 

2) . Transfer is available between every two lines whose hash value is “1” in the same 
column. 

Actually, the purpose of constructing the hash table of line-station is to perform high speed 
search utilizing these two relationships. 

3.2.3. Line-station searching linked-list 

Line-station searching linked-list stores the stations that every line runs through succes- 
sively in a duplex linked-list. The list is sorted by line NO. to locates elements with binary 
search. Line-station search list structure is schematically shown in Fig. 3. 




nr r V i s - 



Figure 3: Schematic view of line-station searching linked-list structure 
Since the time complexity of binary search is 0(log(n)) , the total time complexity of the 
line-station duplex searching linked-list is 

0(log(m)) + 0(log(max(spZ))), (4) 

where m is the count of all lines; ma x(spl) is the max count of stations that a line runs 
through. 



3.2.4 Station-line searching linked-list 

The station-line searching linked-list is also designed as a duplex linked-list, every of whose 
nodes stores all lines that run through the station. The list is sorted by station NO. to locates 
elements with binary search. Line-station search list structure is schematically shown in Fig. 



4. 
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Figure 4: Schematic view of station- line searching linked-list structure 
Similarly, the total time complexity of the station-line duplex searching linked-list is 
0(log(n)) + 0(log(max(7ps))) (5) 

where n is the count of all stations; max(Zps) is the max count of lines that runs through a 
station. 
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3.3. Implementation of algorithm 

The following route searching steps are acquired with integrated consider of search-level 
control strategy and optimization methods: 

With regard to the start station s s tart. and the destination station Sdest 

1) . Set L s tart and set Ldest are obtained from station-line searching linked-list by binary 
searching and select the one with less elements, that is min {L start, Ldest) ■ 

2 ) . Pick out all elements from min (L sta rt, Ldest) one by one and search for s s tart or Sdest 
in line-station searching linked-list (If min (L sta rt, Ldest) = L star t , search for s s tart , else search 
for Sdest), record matched lines and store them in set L mat ch', 

a) . If L mat ch 7 ^ 0 , output all elements in L matc h which are direct lines; 

b) . Else if L ma tch — 0, deepen the search levels and continue with step 3); 

3 ) . Figure out D A start bl A dest 

a) . If D ^ 0 and D = {d\ , then d \ , d% , ■■■d n are transfer stations. Go to step 1), 

perform direct-line searching with s s tart and elements in D, as well as elements in D and Sdest 
respectively. 

b) . Else if D = 0 deepen the search levels and continue with step 4); 

4) . Pick out s' start and s' start from A s tart and Ad es t respectively, then set the searching 
level to 1 and continue with step 1 ), namely transform the problem to search for direct-line 
between s' start and s' dest . 



§4. Summary 

Because of the characteristics of intersecting and overlapping, the commonly used router 
algorithms based on Graph Theory or Genetic Algorithm are inapplicable in the public transport 
network; therefore, a set-based connectivity judging method was proposed in this article. Search 
level control strategy is adapted to resolve the graph problem with huge number of vertexes 
and edges and a hash table and two duplex linked-list are designed to improve the searching 
speed and reduce the searching complexity. 
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Abstract In this paper, some elementary methods are used to study the property of the 
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§1. Introduction and result 



For any complex number s, let 

OO 

C(s) = 5>- S 

k — 1 

be the Riemann zeta function. For any positive integer n, let T n be a positive real number such 
that 

7 r 2n 

C(2 n) = — , (1) 

-L n 

where 7r is ratio of the circumference of a circle to its diameter. Then the sequence T = {T n }'^L 1 
is called the Smarandache-Riemann zeta sequence. About the elementary properties of the 
Smarandache-Riemann zeta sequence, some scholars have studied it, and got some useful results. 
For example, in [2], Murthy believed that T n is a sequence of integers. Meanwhile, he proposed 
the following conjecture: 

Conjecture. No two terms of T n are relatively prime. 

In [3], Le Maohua proved some interesting results. That is, if 

ord{ 2, (2n)!) < 2n — 2, 



where ord( 2, (2n)!) denotes the order of prime 2 in (2n)!, then T n is not an integer, and finally 
he defies Murthy’s conjecture. 

In reference [4], Li Jie proved that for any positive integer n > 1, we have the identity 



ord{ 2, (2n)!) = ot 2 { 2 n) 




= 2n — a(2n, 2), 



where [x] denotes the greatest integer not exceeding x. 

So if 2 n — a(2n, 2) < 2 n — 2, or a(2n, 2) > 3, then T n is not an integer. 
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In fact, there exist infinite positive integers n such that a(2n, 2) > 3, and T n is not an 
integer. From this, we know that Murthy’s conjecture is not correct, because there exist infinite 
positive integers n such that T n is not an integer. 

In this paper, we use the elementary methods to study another property of the Smarandache- 
Riemann zeta sequence, and give a general result for it. That is, we shall prove the following 
conclusion: 

Theorem. If T n are positive integers, then 3 divides T n , more generally, if n = 2k, then 
5 divides T n ; If n = 3 k, then 7 divides T n , where k 0 is an integer. 

So from this Theorem we may immediately get the following 

Corollary. For any positive integers m and n(m ^ n), if T m and T n are integers, then 



§2. Proof of the theorem 

In this section, we shall complete the proof of our theorem. First we need two simple 
Lemmas which we state as follows: 

Lemma 1. If n is a positive integer, then we have 



{Tm,T n ) > 3, (T 2m ,T 2n ) > 15, {T^mjTsn) > 21. 




(2) 



where B 2n is the Bernoulli number. 

Proof. See reference [1]. 

Lemma 2. For any positive integer n, we have 




p—l\2n 



( 3 ) 



where I n is an integer and the sum is over all primes p such that p — 1 divides 2 n. 
Proof. See reference [3]. 

Lemma 3. For any positive integer n, we have 




( 4 ) 




C(2n) = (-l) n 1 , -B 2n , n > 1, 

(2 n)\ 



( 5 ) 



where 




(6) 



Using (1), (5) and (6), we get (4). 

Now we use above Lemmas to complete the proof of our theorem. 
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For any positive integer n, from (4) we can directly obtain that if T n is an integer, then 3 
divides T n , since ( a n ,b n ) = 1. 

From (1), (2) and (3) we have the following equality 



C(2n) = — = (-ir+i. ^) 2 ^ 

U j T n 1 j 2(2n)! 



= (-l) 



n+1 



(27r) 2n 

2(2n)! 



£ \ 

p—l\2n 



Let 



P = PiP2- ■ -Ps, 

p—l\2n 

where Pi (1 < * < s) is a prime number, and pi <p% - ■ ■ < p s . 

Then from the above, we have 

(-l) n+1 • 7r 2n (-l) n+1 • (2n)! 



T n = 



(27r) 2n 

2(2n)! 



In- V 



'■W- X 



p—l\2n J 

(-l)"+ 1 -(2 n)!- JI P 

p—l\2n 



p— l|2n J 



In Up Up E ! 

p—l\2n p—l\2n p—l\2n / 
(_l)n+t . (2 n)\ -pip 2 • • -Ps 



2 2 ” 1 • (in 'PlP2 


( 1 

' "Ps -P 1 P 2 ■ ■ ■ Ps ■ 1- 


1 

b ’ ' 


- 1 )) 


V 


\p 1 


P2 


Ps) ) 




(— l) n+1 • (2n)! • Pl p 2 • • 


■Ps 




2 2n ~ 1 ■ (I n ■ P 1 P 2 ' 


"Ps- P 2 P 3 -■ Ps - P 1 P 3 • • 


■Ps - ' 


■ ■ -P 1 P 2 • • -Ps-l) 


Then we find that if Pi IP 1 P 2 ■ • 


• p s , 1 < i < s, but 







( 7 ) 



Pi t (In -PlP2 ■■■ Ps P2P3 ■■■Ps PlP3 • • -Ps P 1 P 2 ■ ■ -Ps-l) ■ 

So we can easily deduce that if T n are integers, when n = 2k, 5 can divide T n ; While 
n = 3 k, then 7 can divide T n . 

This completes the proof of Theorem. 
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Abstract In this paper, a system of generalized set-valued variational inclusions in Banach 
spaces is introduced and studied, which includes many variational inclusions studied by others 
in recent years. By using some new and innovative techniques, several existence theorems 
for the system of generalized set-valued variational inclusions in g-uniformly smooth Banach 
spaces are established, and some perturbed iterative algorithms for solving this kind of system 
of set-valued variational inclusions are suggested and analyzed. Our results improve and 
generalize many known algorithms and results. 
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§1. Introduction 

In recent years, variational inequalities have been extended and generalized in different 
directions, using novel and innovative techniques, both for their own sake and for the applica- 
tions. Useful and important generalizations of variational inequalities are set-valued variational 
inclusions, which have been studied by [1-9] . 

Recently, in [1], Chaofeng Shi introduced and studied the following class of set- valued 
variational inclusion problems in a Banach space E. 

For a given m-accretive mapping A, B, C, G : E — > CB(E), nonlinear mappings N(-, •), M(-, •) 
E x E — > E, single- valued mapping g : H — > H, any given / £ E and A > 0, find u £ E,x £ 
Au, y £ Bu, z £ Cu, v £ Gu such that 

/ £ N(x,y) - M(z,v) + Xw(g(u)), (1.1) 

where CB(E) denotes the family of all nonempty closed and bounded subsets of E. Under the 
setting of q— uniformly smooth Banach space, Chaofeng Shi [1] gave the existence and conver- 
gence theorem for the solution of variational inclusion (1.1). 

1 This research is supported by the natural foundation of Shaanxi province of China(Grant. No.: 2006A14)and 
the natural foundation of Shaanxi educational department of China(Grant. No.:07JK421) 



Vol. 4 The system of generalized set-valued variational inclusions in q-uniformly smooth Banach spaces 51 



For a suitable choice of the mappings A, B, C , G, W, N , M, g and / £ E, a number of known 
and new variational inequalities, variational inclusions, and related optimization problems can 
be obtained from (1.1). 

On the other hand, X. P. Ding [2] introduced and studied the system of generalized vector 
quasi equilibrium problems in locally G-convex uniform spaces, and the existence results for 
the system of equilibrium problems is obtained. 

Inspired and motivated by the results in Chaofeng Shi [1], S. S. Chang [2], and X. P. 
Ding [3], the purpose of this paper is to introduce and study a class of system of set- valued 
variational inclusions. By using some new techniques, some existence theorems and approximate 
theorems for solving the set- valued variational inclusions in (/-uniformly smooth Banach spaces 
are established and suggested. The results presented in this paper generalize, improve and unify 
the corresponding results of Chaofeng Shi [1], S. S. Chang [2], Ding [3,4], Huang [5,6], Zeng [7], 
Kazmi [8], Jung and Morales [9], Agarwal [10], Liu [11,12] and Osilike [13,14]. 



§2. Preliminaries 



Definition 2.1. Let / be a finite or infinite index set. Let {Ei} ie j be a family of real 
Banach spaces, A,B,C,G : \\Ei — > CB(Ei)be set-valued mappings, Wi : D(W) C Ei — > 

iei iei 

2 Ei be an set-valued mapping, : E t — > E, be a single-valued mapping, and A/fy-, •), Mfy, ■) : 
Ei x Ei — > Ei be two nonlinear mappings, for any given /, £ Ei and A > 0, we consider the 
following problem of finding u £ £ Au,y £ Bu,z £ Cu,v £ Gu such that, for each 

ie/ 

i £ I, 



fi £ Ni(PiX,Piy) - Mi(PiZ,P z v ) + \Wi(gi(Piu)), 



(2.1) 



where Pi : ^ is a projection operator. This problem is called the generalized set- 

iei 

valued variational inclusion problem in a Banach space. 

Next we consider some special cases of problem (2.1). 

(1) If I = = f,Ni = N,Mi = M,Wi = W,gi = g, then problem (2.1) is equivalent 

to 

f £ N(x,y) - M(z,v) + Xw(g(u)), (2.2) 



This problem is introduced and studied by Chaofeng Shi [1], 



(2) Furthermore, if M = 0, W = A : D(A) — ► 2 E is an m-accretive mapping, A = T, B = F 
and C = G = 0 , then problem(2.1) is equivalent to finding q £ E,w £ Tq, v £ Fq such that 

f £N(w,v) + XA(g(q)). (2.3) 

This problem was introduced and studied by S. S. Chang [1], 



For a suitable choice for the mappings A,B. C,G, W t . N t . Mj, </j, /, and the family of spaces 
Ei , we can obtain a lot of known and new variational inequalities, variational inclusions and the 
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related optimization problems. Furthermore, they can make us be able to study mathematics, 
physics and engineering science problems in a general and unified frame [1-9]. 

Definition 2.2. [10] Let E be a real Banach space. The module of smoothness of E is 
defined by 

p E {t) = sup{^(||x + y\\ + Hz - y||) - 1; ||z|| < 1, ||y|| < t}. 

The space E is said to be uniformly smooth if lim (p E (t)/t) = 0. Moreover, E is called q- 
uniformly smooth if there exists a constant c > 0 such that p E (t) < ct q . 

Remark 2.1. [10] All Hilbert spaces, L p ( or l p ) spaces ( p > 2) and the Sobolev spaces, 
W^(p > 2) are 2-uniformly smooth, while, for 1 < p < 2, L p (or l p ) and spaces are 
p— uniformly smooth. 

Definition 2.3. Let S : E — > CB(E) be a set-valued mapping. S is said to be quasi- 
contractive, if there exists a constant r £ (0, 1) such that for any p £ Sx, q £ Sy , 

lb - 911 < r max{||z - y ||, ||z - p||, ||z - q\\,\\y - p||, \\y - g||}. 

To prove the main result, we need the following lemmas. 

Lemma 2.1. [15] Let E be a q— uniformly smooth Banach space with q > 1. Then there 
exists a constant c > 0 such that 

II tx + (1 - t)y - z\\ q < [1 - t(q - l)]|b - z\\ q + tc\\x - z\\ q 

-f(l-t‘ ? - 1 c)||x-?/|| 9 , 

for all x,y,z £ E and t £ [0, 1]. 

Lemma 2.2. [16] Suppose that {e n }, {/„}, {g n } and {■jn} are nonnegative real sequences 
such that 

^n+l ^ (l fn)&n fn9n T ^ 0; 

oo oo 

with {f n } C [0, 1], V' f n = oo, lim g n = 0 and V' 7„ < oo. Then lim e n = 0. 

A ' n — >-oo z ' n — >oo 

n— 0 n—0 

For the remainder of this paper, r and c denote the constants appearing in definition 2.3 
and lemma 2.1. We assume that is the Hausdorff metric on CB(E) defined by 

H(A, B) = max {sup d(x, B), sup d(y, A)}, A, B £ CB(E). 

xeA y£B 

Lemma 2.3. [17] Let E be a complete metric space, T : E — > CB(E) be a set-valued 
mapping. Then for any given e > 0 and any given x,y £ E,u £ Tx , there exists v £ Ty such 
that 

d(u,v ) < (1 + e)H(Tx,Ty). 

Lemma 2.4. [18] Let X and Y be two Banach spaces, T : X — > 2 5 a lower semi-continuous 
mapping with nonempty closed and convex values. Then T admits a continuous selection, i. e., 
there exist a continuous mapping h : X —> Y such that h{x) £ Tx, for each x £ X. 

Using Lemma 2.3 and Lemma 2.4, we suggest the following algorithms for the generalized 
set-valued variational inclusion (2.1). 
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Algorithm 2.1. For any given xq G n Ei, x' 0 G Ax 0 ,y' 0 G Bx 0 ,z' 0 G Cx 0 ,v' 0 G Gx 0 , 

iei 

compute the sequences {x n } and {y n } by the iterative schemes such that, for each i G I, 



(1) PiXn+l G (1 Oir^)PiX n + Otn(fi + P^n 



Ni(PiX n , Piy n ) + Mi(PiZ n , PiV n ) 



A WMPiVn))), 



(2) Piy n G (1 - Pn)PiX n + (3 n (fi + PiX n - Ni(PiX n , Piy n ) + Mi(PiZ n , Pi.v'J - \Wi(gi(PiX n ))), 

(3) x n G Ay n , ||r > jX n Pix n -\- 1 || + (1 + ■ -)Hi(PiAy n , PiAy n _j_i), 

n+1 



(4) yn G By n , \\P z y n Piy n - fi|| ^ (1 + .. )Hi(PiBy n , PiBy n +\) , 

n+1 



(2.4) 



(5) z n G Cy n , II PiZ n - PiZ„+i|| < (H —)H i (P i Cy n ,P i Cy n+ 1 ), 

n+1 

(6) v n G Gy n , || PiV n - P+ n+ i|| < (H l —)H i (P i Gy ni PiGy n+ 1 ), 

n+1 

(7) x’ n G Ax n , || PiX n - Pix' n+ 1 || < (H 1 —)H l (P i Ax n , PiAx n+ i) , 

n + 1 

(8) y' n G Bx n , ||Pij4 - Piy' n+1 \\ < (1 + —j— )H i (P i Bx n ,P i Bx n+1 ), 

n+1 

(9) 4 G Cx n , || P iZ ' n - P lZ ' n+l || < (1 + - l —)H l (P l Cx n ,P l Cx n+1 ), 

n+1 



(10) v' n G Gx n , || P iV ' n - Piv' n+1 \\ < (1 + — - - )H i (P i Gx n ,P i Gx n+1 ), 

n+1 



n = 0, 1, 2, • 



The sequence { x n } defined by (2.4), in the sequel, is called Ishikawa iterative sequence. 

In algorithm 2.1, if f3 n = 0 , for all n > 0 , then y n = x n . Take x n = x' n , y n = y' n , z n = z' n 
and v n = v' n , for all n > 0 and we obtain the following. 

Algorithm 2.2. For any given Xq G n Ei, Xo G Ax 0 ,yo G Bx 0 ,z 0 G Cx 0 ,v 0 G Gx 0 , 

iei 

compute the sequences {x n }, {x n }, { y n }, {z n } and {r> n } by the iterative schemes such that, for 
each i G I, 



PiXn -\- 1 G (1 Otn)PiX n + O n (/ ) + PiX n ^i(^PiX n , PiVn) + A4j(P^ n , PiVn) AIT) (^j (P^X^))) , 

Xn G Ax n , ||PjX n Pi^n+l|| ^ (1 + , z)Hi(PiAx n , PiAXn+x) , 

n+1 

y n G Bx n , II Piy n - PiVn+ 1 || < (H ^—)H i (P i Bx n ,P i Bx n+1 ), (2.5) 

n+1 

Zn G Cx n , || PiZ n P+n +1 I ^ (1 "t" , z)Hi(PiCx n , PiCx n -\- i), 

n + 1 



v n G Gx n , ||Pjrj n Pi+ n +i| + (1 + )Hi(PiGx n , PiGx n -\.\ ), 

n+1 



n = 0, 1, 2, ■ 



The sequence {x n } defined by (2.5), in the sequel, is called Mann iterative sequence. 
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§3. An existence theorem for the system of the generalized 
set-valued variational inclusions 

In this section, we shall establish an existence theorem for solutions of the system of the 
set- valued variational inclusions (2.1). We have the following results. 

Theorem 3.1. Let E be a (/-uniformly smooth Banach space, A,B,C,G : — > 

iei 

JJCP(p) be four set- valued mappings, for each i £ I, ip (•,•), Mj(-, •) : Ex E — > E,gi : Ei — > 
iei 

Ei,Wi : Ei — > Ei be four single- valued mappings. If I-Ni{PiA{-), PjP(-))+Mj(PjC(-), PjG(-)) — 
AWi(gi(Pi(-))) is quasi contractive, then there exists u £ JJp, ® £ -Au,?/ £ Bu,v £ Cm, 0 £ 

ie/ 

Cm, which is a solution of system of the generalized set- valued variational inclusion (2.1). 
Proof. For each i £ I , let 



S(Pu) = P« - Ni(PiA(u),PiB(u)) + Mi(PiC(u) , PiG(u)) - XW^Piu)) , 

Si(Piu) = fi + S(Piu). 

Since S is set-valued quasi contractive, S± is also set-valued quasi contractive. 

It follows from Naddler [17] that Si has a unique fixed point PiU £ Ei, that is 

Pm £ Si (Pm) =fi + P z u - Ni{PiA(u),PiB(u)) + M,(PC(m),PG(«)) - AW^Pm)). 

Thus, there exist x £ Au, y £ Bu, z £ Cm, v £ Gu such that, for each i £ I, 

fi £ Ni{PiX,Piy) - Mi{PiZ,Piv) + XWi(gi(Piu)), 

that is, m £ E is a solution of the generalized set-valued variational inclusion (2.1). 



§4. Approximate problem of the system for generalized 
set-valued variational inclusions 

In theorem 3.1, under some conditions, we have proved that there exist u £ n Ei,x £ 

ie/ 

Au, y £ Bu, z £ Cm, v £ Gu which is a solution of the system of generalized set-valued varia- 
tional inclusion (2.1). In this section, we shall study the approximate problem for the solution 
of the system of generalized set- valued variational inclusions (2.1). We have the following result. 

Theorem 4.1. Let E, be a family of (/-uniformly smooth Banach spaces, A,B,C,G : 

nP->n CB(Ei) be four set- valued continuous mappings, for each i £ I, 
i£i iei 

Ni (-,-), Mi (-, ■) : Ei x Ei Ei,gi : Ei —> Ei,Wi : Ei Ei be four single- valued continuous 
mappings, and {a n }, {/3 n } be two sequences in [0,1] satisfying the following conditions: 

(i) I - Ni(PiA(-),PiB(-)) + Mj(pC(-), PG(-)) - AWi(ff*P(-)) is quasi contractive, 

(ii) for each i £ I, A,B,C,G are M-Lipschitz continuous with the constants yi,/i 2 ,y 3 and /Z 4 , 
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(iii) 0 < h < a n , n > 0. 0 < //* < 1/2, i = 1,2, 3, 4, 

(iv) a„(g — 1 — cr q ) < 1 ,cr q < q— 1,/3 9_1 < l/c(l - cr 9 ), 

ca 9_1 + cr q (3 n (cr q — q A 1) < 1 — cr q . 

Then for any given xq £ £ Ax o,i/q £ Bxq,z' 0 £ Cxq,v' 0 £ Gx o, the sequences 

ie/ 

{x„}, {x n }, {y n }, {z n } and {f n } defined by algorithm 2.1 strongly converge to the solution 
u £ £ Au,y £ Bu,z £ Cu, v £ Gu of the system of generalized set- valued variational 

ie/ 

inclusions (2.1) which is given in Theorem 3.1, respectively. 



Proof. For each i £ I, hr (1) and (2) of (2.4), choose h n £ Wi(gi(PiX n )),k n £ WpgpPtyn)), 
such that 

Pi x n - (-1 (1 O n ) Pi<En A &n(fi A PiPn Vj(P^X n , PiPn ) T ( Pj Z rh , PiV n ) Afc n ), 

PiVn = (1 - Pn)PiXn A (3 n {fi + P i X n - N z (Pix' n , Piy'J A Mi(PiZ' n , Piv' n ) - Xh n ) . 



Let 



Pn — fi A PiPn -Vj (P^X^, PiVn) A Mi(PiZ n: Pi'V n ') \ k n , 
r n = fi A PiX n - Ni(Pix' n: Piy' n ) A Mi(Piz' n , P l v' n ) - A h n . 



Then 



Pi x n -\- 1 (1 ®-n)Pi x n A OL n Pm 

PiVn = (1 ^ Pn)PiX n A f3 n r n - (4.1) 

Since I — V i (P i A(-), P,l?(-)) + Mj(PjC(-), PjG(-)) — A Wi(giPi(-)) is quasi contractive, 

|| Pn - Piu\\ < r max{\\Piy n - PiU ||, ||p n - Piy n ||}, 

which implies that 

II Pn - Piu\\ q < r q (\\Piy n - Piu\\ q A \\p n ~ PiVn\\ q )- (4.2) 

Note that 

I bn - || < rd{PiX n ,Piy n ),n > 0. 

We consider the following cases. 

Case 1: Suppose that d(PiX n , Piy n ) = \\PiX n — Piy n \\ for some n > 0. It follows from (4.1) that 
| bn - r n\\ < r\\PiX n - Piy n || = r\\/3 n (PiX n - r„)|| 

< r(3 n \\PiX n - r„||. (4.3) 

Case 2: Suppose that d(PiX n , Piy n ) = ||Pja; n — r n || for some n > 0. Then we have 

\\p n ~ r n \\ < r\\PiX n - r n \\. (4.4) 
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Case 3: Suppose that d(PiX n , Piy n ) = || Piy n — r n \\ for some n > 0. Using (4.1), we have 

I bn - r n \\ < r\\Piy n - r n \\ 

= r || (1 - (3 n )(PiX n - r n )|| 

= r( 1 - 0 n )\\ PiX n - r„|| < r\\PiX n - r n ||. (4.5) 

Case 4: Suppose that d(PiX n , Piy n ) = \\PiX n — p n || for some n > 0. Then we have 

\\Pn-r n \\ <r\\PiX n -p n \\. (4.6) 

Case 5: Suppose that d(PiX n , Piy n ) = || Piy n ~ Pn\\ for some n > 0. It follows from (4.1) that 

I bn - r n \\ < r\\P l y n - p n || 

= r || (1 - /3 n )(PiX n - Pn) + Pn{r n Pn) || 

^ ( 1 Pn) I Pi^n Pn 1 1 4“ ^ Pn \\Pn 1 1 ? 

which implies that 

I bn - r n \\ < r ^_ ^ \\PiX n -p n || < r\\PiX n -p n \\- (4.7) 

T Pn 

It follows from (4.3)-(4.7) that 

I bn ~r„\\ q < r q \\PiX n p„ || +r q \\PiX n -r n \\ q , 

n> 0. (4.8) 

It follows from Lemma 2.1 that 

|| PiVn - PiU\\ q = ||(1 - Pn)PiXn + PnT n ~ PiU\\ q 
= ||(1 - P n )(PiX n - PiU) + Pn(r n - PiU) \\ q 
< [1 - Pn(q - l)]||-Pi* n - Piu\\ q + p n c\\r„ - PiU\\ q 

-Pn( 1 - Pn~ lc )\\ p iX n - r n \\ q ,n> 0. (4.9) 

Similarly, we have 

|| PiVn ~Pn\\ q < [1 - Pn(q ~ l)]\\PiX n - p n \\ q + P n c\\p n ~ X n || 9 

~Pn{ 1 - Pn~ lc )\\ P iXn ~ r n \\ q ,n > 0. (4.10) 

and 

\\PiX n+1 - Piu\\ q < [1 - a n {q - l)]||P l x n - Piu\\ q + a n c|bn - Piu\\ q 

-a«( 1 - < _1 c)|| PiX n -p n || 9 ,n > 0. (4.11) 



By virtue of the condition (iii) and (iv), (4.2) and (4.8)-(4.11), we have 
\\PiX n+ i - Piu\\ q < [1 - a n (q - l)]||Pia; n - Piu\\ q 
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+a n cr q [\\Piy n - Piu\\ q + \\Piy„ - p n || 9 ] 

-a n (l - a q ^ 1 c)\\PiX n — p n \\ q 
< [1 - a n (q - 1) + cr q a n (l - fi n (q - l))]||Pj£ n - Piu\\ q 
+c 2 r q a n fi n \\r n - Piu\\ q + c 2 r q a n /3 n \\p n - r n \\ q 
+[cr q a n ( 1 - fi n (q - 1)) - a n (l - a® _1 c)]||Px n - p n \\ q 
-2cr q a n fi n (l - fi q ~ x c)\\PiX n - r n \\ q 
< [1 - a n {q - 1) + cr q a n (l - fi n (q - 1)) + c 2 r 2q a n fi n \\\PiX n - Pm\\ q 
+2 cr q a n fi n (cr q + {fi q ~ l c - l))||Pj 2 ; n - r n \\ q 
+a n [c 2 r 2q l3 n + cr q ( 1 - fi n (q - 1)) - (1 - a q l ~ x c)]\\P i x n - p n \\ q 
< [1 - a n (q - 1 - cr q )( 1 + cr q /3 n )]\\PiX n - Piu\\ q 
< [1 ~ u n (q - 1 - cr q )]\\PiX n - P iU \\ q , 

Set e n = II PiX n - PiM|| 9 , fn = a n {q - 1 - cr q ), 

3n = 7n = 0,n > 0. 

It follows from the conditions (iii), (iv) and Lemma 2.2 that 



that is, 



Note that 



lim e n = 0, 

n — kx) 



lim ||PjX„ - Pi«|| = 0. 



which implies 

Thus, from (4.1), we have 



r n - P iU \\ < rmax{||PjX„ - Pju||, ||r„ - PjX n ||}, 
r n -*• P l u(n -> oo). 



PiVn (1 fin) P fin^n * Pi'U’- 

From (2.4) and conditions (ii) and (iii), we have 



\PiX n - PiX n+ 1 || < (1 H —)Hi(PiAy n , PiAy n+1 ) 

n+l 



< (1 + — + ^HWPiVn ~ PiVn+ 1||, 

n + l 

which implies that {PiX n } is a Cauchy sequence in E,. 

So, there exists x £ Ei such that PiX n — > PiX. Now we show that x £ Au. 



In fact, 



i£l 



d(P i x,P l Au) < || PiX - PiX n || + Hi{P r Ax n ,PiAu) 
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< II PiX - PiX n || + n(l H — )|| PiX n - P iU \\ -> 0. 

n + 1 



Also, we have 



Vn -> y £ Bu , z n z € Cu, v n — > t; £ Gm. 

From (2.4), we have, for each ?' £ J, 

PiU £ j) + Piit - Ni(PiX, Pi-y ) + Mi(PiZ , 



that is, 



/i £ Ni{PiX,Piy) - Mi(PiZ,PiV ) + AWj^Pj-u)). 



The required results. 



Remark 4.1. Theorem 4.1 generalizes Theorem 4.1 in Chaofeng Shi [1], the corresponding 
results in S. S. Chang [2], Agarwal [10], Liu [11, 12], Osilike [13, 14] and others. 

Remark 4.2. Since algorithm 2.2 is a special case of algorithm 2.1, from Theorem 4.1, we 
can obtain the convergence Theorem for algorithm 2.2, the details are omitted. 
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solution of matrix equation AXB = C 
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Tianshui, Gansu, 741001, P. R. of China 
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Abstract For generalized reflexive matrices R, S, i.e., R T = R, R 2 = I, S T = S, S 2 = I, 
matrix A is said to be (R,S)-symmetric matrix when RAS = A. In this paper, an iterative 
method is established for solving the matrix equation AXB = C with (R,S)-symmetric so- 
lution X. By this method, the solvability of the above matrix equation can be determined 
autonomically, and a solution or the least-norm solution of which can be obtained within 
finite iteration steps. Meanwhile, the optimal approximation solution to a given matrix Xo 
can also be derived by the least-norm solution of a new matrix equation AXB = C, where 
X = X — Xo, C = C — AXqB. Finally, given numerical examples illustrate the efficiency of 
the iterative method. 

Keywords Matrix equation, iterative method, (R,S)-symmetric solution, least-norm solution, 
optimal approximation solution. 



§1. Introduction 

Let R™ xn be the set of all m x n real matrices with rank k, SR nxn , OR nxn be the 
set of all symmetric matrices, orthogonal matrices in R nxn , respectively. Denoted by the 
superscripts T and I n the transpose and identity matrix with order n, respectively. For matrices 
A = (ai, 02 , . . . , a n ),B £ R mxn , a,i £ R rn , 38(A) and tr(A) represent its column space and trace, 
respectively; The symbol vec(-) represents the vec operator, i.e., vec(A) = (af , o^, . . . , aJ) T ; 
A 0 B stands for the Kronecker product of matrices A and B ; Moreover, ( A , B) = tr(B T A) 
is defined as the inner product of the two matrices, which generates the Frobenius norm, i.e., 
|| A|| = y/< A, A > = ^tr(A T A). 

Definition 1.1. [ 1,2 lFor given generalized reflexive matrices R G S G f?” x ”, i.e., 

R T = R, R 2 = /, S T = S, S 2 = /, we say matrix A £ R mxn is (R,S)-symmetric ((R,S)-skew 
symmetric), if RAS = A (RAS = —A). 

The set of all m x n (R,S)-symmetric ((R,S)-skew symmetric) matrices with respect to 
(R,S) is denoted by GSR mxn (GSSR mxn ). 

Definition 1.2. Assume M,N £ R pxm , where p,m are arbitrary positive integers, if 
(M, N) = 0, matrices M, N are called to be orthogonal each other. 



1 This work is supported by the Science Fundations of Tianshui Normal University (TSB0819 and TSB0721) 
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Remark 1.1. 

(a) In this paper, we let R , S be generalized reflexive matrices as in Definition 1.1. 

(b) For VF G GSR mxn , G G GSSR mxn , we can verify that (F, G) = 0. 



For above generalized reflexive matrix R, R T = R implies that the eigenvalues of R belong 
to real field, and the absolute of the eigenvalues equal to 1 since R 2 = I. Certainly, S has 
similar properties to R. Hence, we have the following assertion. 

Lemma 1.1. For given generalized reflexive matrices R , S as in Definition 1.1, there exist 
unitary matrices U\ G OR mxm , and U 2 G OR nxn , such that 



R = Lh 



/ Ir 


0 1 


1 U?, s = u 2 1 


f II 


0 


\ o 


Im—r J 


I 1 


{ 0 


T 



uL 



From Definition 1.1 and Lemma 1.1, we generate the following conclusion. 
Lemma 1.2. Suppose X G GSR mxn , then 

' x x o s 



X = U, 



U 2 



0 A'o 



where VIiG R rxl , X 2 G 

From Lemma 1.2, we can easily get a (R,S)-symmetric matrix by choosing different Ai, 
but Ui ( i = 1, 2) are fixed in Lemma 1.1. 

The problems to be discussed in this paper can be expressed as follows: 

Problem I. For given generalized reflexive matrices R G R mxm , S G R nxn , and A G 
RP xm , B G R nxq , C G R pxq , find X G GSR mxn , such that 



AXB = C. 



(1) 



Problem II. When Problem I is consistent, i.e., its solution set Se is not empty, for given 
matrix Xq G R mxn , find X G Se such that 

||X-X 0 ||= min \\X-X 0 \\. (2) 

A Eo.e 

The linear matrix equation (1) has been discussed widely with unknown dC 3 ~ 7 l , where 
X has particular structures, such as, symmetric, skew-symmetric, reflexive, anti-reflexive, Her- 
mitian etc.. In these literatures, the authors obtained the expressions of associated solutions 
and solvability conditions by general inverse or matrix decompositions : the singular value 
decomposition (SVD), the generalized SVD (GSVD)! 8 ! or the canonical correlation decompo- 
sition (CCD)! 9 !. In complex field, the inverse eigenvalue problem and procrustes problem on 
(R,S)-symmetric matrices have been studied by William [1,2]. Of course, the (R,S)-symmetric 
solution of matrix equation (1) can also be given by using GSVD or CCD. In addition, by 
establishing iterative algorithm, the symmetric solution, skew-symmetric solution, Hermitian 
least-norm solution of matrix equation (1) have been derived in [10-12], respectively. Moreover, 
the general solution and symmetric solution of matrix equations A\XBi = C\, A 2 XB 2 = C 2 
have also been investigated by Sheng [13] and Peng [14], respectively. 

Problem II is the optimal approximation problem, which occurs frequently in experimental 
design^ 13 !. j n equality (2), the matrix A'o may be obtained from experiments, but it is not 
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necessary to be needed form and the minimum residual requirement. The nearness matrix X 
is the matrix that satisfies the needed form and the minimum residual restriction. 

Motivated by the iterative methods mentioned in [10], in this paper, we will construct an 
iterative method to obtain the (R,S)-symmetric solutions of Problem I and II. This work will 
be accomplished in section 2 and section 3, respectively. In section 4, we will provide some 
numerical examples to illustrate the efficiency of the iterative method. 



§2. The iterative method for Problem I 



In this section, we will propose an iterative algorithm to solve Problem I, some lemmas 
will be given to analyze the properties of this algorithm. For any (R,S)-symmetric matrix, We 
will show that the solution of matrix equation (1) can be obtained within finite iterative steps. 
The iterative method of Problem I is stated as follows: 

Algorithm 2.1. 

Step 1: Input matrices R G R mxm , S G R nxn , and A G R pxm , B G R nxq C G R pxq . 
Choosing arbitrary matrix X\ G GSR mxn . 

Step 2: Calculate 



Ri = C — AXiB, 



Pi = ~{A T R l B T + RA t R 1 B t S ), 

k := 1. 

Step 3: If Rk = 0, stop; otherwise goto next step. 
Step 4: Calculate 

11 Rk 112 



Xfe+1 Xk ' || P k || 2 
Step 5: Calculate 



Pk- 



Rk + 1 — C — AX k+1 B — R k — | p 1, 2 AP k B. 



P k+1 = ~(A T R k+1 B T + RA T R k+l B T S) + 



l-Rfc+i 






-Pk- 



Algorithm 2.1 reveals that X i: Pi G GSR mxn (* = 1,2,...). In the sequel, we prove the 
feasibility of this algorithm. 

Lemma 2.1. The sequences Ri, Pi {i = 1,2,...) generated by Algorithm 2.1 satisfy that 



(Pj+i, Ri) — (Ri,Ri) — 



\JU 

IIP,: I 



-{Pi, Pi), 



when j > 1, 



(Ri+i,Rj) = (Ri,Rj) - 



\Pi ir \\Pi\\ 2 \\p.i-i\ 12 



( 3 ) 



Proof. When j > 1, from the Algorithm 2.1, noting that RPiS = Pi, we can obtain 

A t R,B t - RA t R 1 B t S'' 



(APiB, Rj) = ( P l ,A T R J B T ) 



p A t R j B t + RA T RjB T S 



Pi, 



2 



2 
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JML P \ + (p. AML 

’ iiif,_iip p ’- l ) + \ F " 2 



RP t S , 



RA T RiB T S 






moreover, 



( 4 ) 



(Ri+i,Rj) = {Ri- l ^(AP i B),R j 



= (Ri,Rj)~j^(AP i B,R j ) 

= (R i ,R j )-^t(P i ,A T R j B T ). 



( 5 ) 



The proof of (5) implies that (3) holds when j = 1. Furthermore, submitting (4) into (5) 
generates (3) holds. The proof is completed. 

Lemma 2.2. The sequences {Pi} generated in the iterative process are self-orthogonal, 
respectively, i.e. , 

(Ri,Rj) = 0, (Pi,Pj) = 0, i,j = 1,2, . ..,k(k > 2), i ± j. (6) 

Proof. We prove the conclusion by induction. Since (F,G) = ( G,F ), we only proof (6) 
holds when i > j. 



When k = 2, it is clear from Lemma 2.1 that 

p 2 Pi> = (R\ , Pi) - |^-jf (Pi.Pi) = 0. 

II p i II 

Applying Algorithm 2.1 yields 



P 2 Pl} = 



/ A T R 2 B t + RA t R 2 B t S \\R 2 

\ 2 + Wi 



PP 



(Pi, Pi 



Pill 

= 0 . 



P2H 2 
pill 2 

p2fPl|| 2 
Pill 2 



= (A t R 2 B t ,P 1 ) + 

= (R 2 ,AP 1 B) + 



\Ri\ 



( 7 ) 



Assume that (6) holds for k = s, that is, ( P s ,Pj ) = 0, ( P s ,Pj ) = 0, j = 1 , 2 ,..., s— 1 . Being 
similar to the proof of (7), we can verify from Lemma 2.1 and the algorithm that (R s + 1 , R s ) = 0, 



and (P S+1 ,P S ) = 0. 

Now, we can finish the proof if (R s+ i,Rj) 
noting that the assumptions and (4), then 

(R s+ i,Ri) = {R.,Ri) - I^LpPsP,#!) 



I R s 
IIP, 



-( P s ,A t R 1 B 1 



0, (P s+ i,Pj) = 0 hold. In fact, when j = 1, 
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I R a 



(Ps,Pl)= 0 



linn 2 

and connecting with (3), we have 

/p / A T R S+1 B T + RA T R S+1 B T S 

\Ps+i, P./ — ( - 



I R 



s+ll 



I Rs 



-P„P1 



/A t R s+1 B t + RA t R s+ 1 B t S „\ , ||i? s+ i|| 2 /n „ 

\ 2 ’ Pl / + -JiV^ {Ps ' Pl/ 



= ( A T R s+1 B T ,P l 
= (-R s +i) AP\B) 

IIPI 112 



IRll 



2 (Pfi+1, (Pi — P 2 )) 



I Rl I 



Pi 



r(Ps+i, P 2 ) — 0. 



(8) 



Furthermore, when 2 < j < s — 1, Lemma 2.1 and the assumptions imply that 

||PJ 2 IIPJ 2 IIP,-" 2 

ll-P, 



(Ps+i, Rj) — {Rs, Rj) || D m 2 (Ps) Pj) + 



°- 



Similar to the proof of (8), we can obtain (P s +:i , Pj) = 0. Hence, the conclusions hold. 

Remark 2.1. From Lemma 2.2, we know that, Ri (i = 1,2,... ,pq) can be regarded as 
an orthogonal basis of matrix space R pxq for their orthogonality each other, so a solution of 
Problem I can be obtained by the Algorithm 2.1 at most pq + 1 iteration steps. 

Lemma 2.3. When matrix equation (1) is consistent, provided that X is an arbitrary 
solution, then Rj-, Pk generated by Algorithm 2.1 satisfy that 

(X-X k ,P k )= || R k || 2 , fc = 1,2, ... . 

Proof. When k=l, from the Algorithm 2.1 and Lemma 2.2, we get 

a-x I ,p l) = (x-x 1 ,AP P±P^hPs' 



A t R 1 B t + RA t R 1 B t S \ / A t R 1 B t + RA t R 1 B t S 

= ( X, ) — ( X\ , 



= ( X 1 A t R 1 B t ) - {X l ,A T R 1 B T ) 

= {A{X-X 1 )B,R 1 ) 

= l|Pi|| 2 - 

Assume that the conclusion holds for k = s, i.e., (X — X S ,P S ) = || R s || 2 , then 



(X-X e+1 ,P S ) = (X-X„- 



I Rs 



Ps 



;Ps,Ps) 



= ((X-X S ),P S )-^X(P S ,P S \ 



= IIpJI 2 -KC||p s II 2 
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= 0, 



and 



A t R s+ 1 B t + RA t R s+ 1 B t S lli?. 



<X-X S+ 1 ,P S+1 > = {X-X a+U - 

= \{x- x s+1 ,A T R s+1 B T + RA t R s+ 1 B t S) + 

= (X-X s+ 1 ,A t R s+ 1 B t ) 

= (A(X-X S+1 )B,R S+1 ) 



ls+l| 



\R 



■S + l| 



\R« 



R, 



-(X-X a+1 ,P S ) 



= IIP 



■S+1| 



Therefore, we complete the proof by the principle of induction. 

Remark 2.2. Lemma 2.3 implies that, if the linear matrix equation (1) is consistent, then, 
Ri = 0 if and only if Pi = 0. In other words, if there exist a positive number t such that Rt ^ 0 
but P t = 0, then matrix equation (1) is inconsistent. Therefore, the solvability of Problem I 
can be determined automatically by Algorithm 2.1 in the absence of roundoff errors. 



Based on the previous analysis, we have the following conclusion whose proof is omitted. 

Theorem 2.1. When matrix equation (1) is consistent, for arbitrary initial matrix X- t g 
GSR mxn , a solution of which can be obtained within finite iterative steps in the absence of 
roundoff errors. 

To facilitate the statement of our main results, we cite the following conclusion. 

Lemma 2.4. (see [10]) Suppose that the consistent linear equations My = b has a solution 
yo £ R(M T ), then yo is the least-norm solution of which. 

Theorem 2.2. When matrix equation (1) is consistent. If initial iteration matrix X\ = 
A T HB t + RA t HB t S, where arbitrary H £ R pxq , or especially, let Xi = 0 £ R pxq , then the 
solution generated by the iteration method is the unique least-norm solution of Problem I. 

Proof. Algorithm 2.1 and Theorem 2.1 imply that, if let X\ = A T HB T + RA T H B T S , 
where H is an arbitrary matrix in R pxq , a solution X* of Problem I is of the form X* = 
A t YB t + RA t YB t S. So it is enough to show that X* is the least-norm solution of matrix 
equation (1). 



Considering the following linear matrix equations with X £ GSR mxn 

I AXB = C 
] ARXSB = C 



(9) 



Obviously, the solvability of (9) is equivalent to that of matrix equation (1). Denote 
vec(X*) = x*, vec(X) = x, vec(Y) = y , vec(C) = c, then the above matrix equations can be 
transformed into 

/ B T ® A 
^ (B T S) ® ( AR ) 




In addition, the solution X* can be rewritten as 
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x . = ( bT ® a ) T (A<=a(( BT ® A ) T ) 

\(B T S)®(AR) ) \y J y\(B T S)®(AR) ) j 

which implies, from Lemma 2.4, that X* is the least-norm solution of the linear matrix equations 
(9). We complete the proof. 



§3. The solution of Problem II 



Suppose that matrix equation (1) is consistent, i.e., Se is not empty. It is easy to verify 
that Se is a closed convex set in matrix space GSR mxn under Fronbenius norm, so the optimal 
approximation solution is unique. Without loss of generality, we can assume the given matrix 



X 0 € GSR ” 



in Problem II because of the orthogonality between (R,S)-symmetric and (R,S)- 



skew symmetric matrix. In fact, for X £ Se, we have 



LY-X 0 f = 



X- 



X 0 + RX 0 S X 0 - RX 0 S 



X - 



Xn + RX 0 S 



X a - RXnS 



Let X = X — Xq, C = G — AXqB, then Problem II is equivalent to find the least-norm 
solution X* £ GSR mxn of the following matrix equation 



AXB = C. (10) 

From Theorem 2.2, if let initial iteration matrix X\ = A T HB T + PA T HB T P, where 
arbitrary H £ R pxq , or especially, let X\ = 0 £ R mxn , we can obtain the unique least-norm 
solution X* of matrix equation (10) by applying Algorithm 2.1. Furthermore, the unique 
optimal approximation solution X to Xq can be obtained by X = X* + Xq. 



§4. Numerical examples 

Example 1 . Input matrices A, B , C, and generalized reflexive matrices R, S as follows: 



l 5 


-3 


0 


3 


0 


2 


8 1 




< 1 


0 


0 


0 


0 


0 


0 


0 


-4 


-6 


4 


-6 


0 


-14 




0 


0 


-1 


0 


0 


0 


0 


-6 


0 


7 


0 


7 


3 


1 




0 


-1 


0 


0 


0 


0 


0 


0 


5 


-3 


-5 


-3 


0 


3 


, R = 


0 


0 


0 


0 


0 


0 


1 


4 


-7 


0 


7 


0 


-8 


-3 




0 


0 


0 


0 


1 


0 


0 


-1 


0 


-6 


0 


-5 


9 


0 




0 


0 


0 


0 


0 


-1 


0 



\0-3 0 3 -7 0 1/ \ 0 0 01000/ 
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B = 



-3 

0 

6 

-2 

-1 

0 



5 

4 

-1 

4 



-5 -2 

9 9 



7 

0 

-2 



5 

4 

-1 

4 

-6 



— 12 \ 
-6 
3 
5 
2 



( 0 1 

1 0 



S = 



0 0-1 



0 0 
0 0 



0 

0 

0 

0 

-1 



c = 



-9 1 


1-9 2 


7 


O 

O 


o 

0 

1 


1 0 


254.7246 


- 420.4859 


- 4.5129 


- 444.0827 


- 420.4859 


95.2199 


- 32.1326 


- 896.3380 


- 642.6184 


- 414.3636 


- 896.3380 


78.8874 


285.7966 


- 281.1569 


169.3670 


- 306.3903 


- 281.1569 


472.0991 


- 219.4309 


907.1669 


255.3442 


601.2957 


907.1669 


- 332.3348 


- 63.6234 


- 937.6031 


- 462.9456 


- 449.7040 


- 937.6031 


149.2335 


352.1195 


- 511.6855 


- 185.5221 


- 550.9943 


- 511.6855 


224.8234 


300.0260 


- 987.7629 


- 247.5968 


- 678.6078 


- 987.7629 


862.4447 



Next, we will find the associated solutions of Problems I and II by Algorithm 2.1 and 
Matlab software. For the influence of roundoff errors, Ri will usually unequal to zero in the 
Fronbenius norm. Therefore, for any chosen positive number e, however small enough, e.g., 
e = l.Oe — 010, whenever ||i? fc || < s, stop the iteration, and X k is regarded as a required 
solution. 

(I) The solutions of Problem I: 

Choose an arbitrary initial matrix 



A , = 



-^36 = 



12.0207 


12.0207 


0 


14.0000 


6.3639 


- 6.3639 


- 6.9999 


21.9996 


- 14.1420 


10.6065 


21.9996 


1.0000 


- 21.9996 


6.9999 


- 14.1420 


- 10.6065 1.0000 


21.9996 


- 34.9993 


- 62.9988 


44.5473 


- 19.0917 - 15.9997 3.9999 


9.8994 


9.8994 


0 


16.0000 


12.7278 


- 12.7278 


18.3846 


- 18.3846 


- 75.0000 


0 


- 24.0414 - 24.0414 


- 62.9988 


- 34.9993 


- 44.5473 


- 19.0917 - 3.9999 


15.9997 


2.1 and 36 iteration 


steps, we 


obtain a solution of Problem I, 


6.3639 


6.3639 


0 


4.0000 


2.1213 - 


2.1213 \ 




0.0221 


7.9080 


- 4.0470 


2.7378 


0.8945 - 


8.1973 




- 7.9080 


- 0.0221 


- 4.0470 


- 2.7378 


- 8.1973 0.8945 




- 0.9779 


- 6.0917 


2.3169 


- 5.7474 


1.8945 3.8025 


5 


4.2426 


4.2426 


0 


- 0.0000 


2.8284 - 


2.8284 




3.5355 


- 3.5355 


8.0000 


0 


1.4142 1.4142 




- 6.0917 


- 0.9779 


- 2.3169 


- 5.7474 


- 3.8025 - 


1.8945 / 





with ||i? 36 || = 3.8765e - Oil < e, and ||X 36 || = 27.3649. 



If let initial matrix Xi = A T HB T + PA T HB T P , where arbitrary H £ R 7x 6 , then the 
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solution Xk generated by Algorithm 2.1 is the unique least-norm solution of matrix equation 



(1) by Theorem 2.2. Choosing arbitrary matrix 










/ -4 


1 


-2 


0 


1 


0 > 




6 


0 


-1 


1 


0 


-1 




—2 


1 


0 


-1 


1 


0 


H = 


1 


0 


0 


2 


— 2 


1 




0 


-1 


0 


-1 


0 


0 




2 


0 


1 


0 


-1 


1 




V o 


-1 


0 


1 


1 


0 / 


then the least-norm solution A'* of Problem I 


is 







/ 6.3639 


6.3639 


0 


4.0000 


2.1213 


-2.1213 


-0.3925 


8.3779 


-4.2023 


3.3713 


-0.6099 


-6.4524 


-8.3779 


0.3925 


-4.2023 


-3.3713 


-6.4524 


-0.6099 


-1.3925 


-5.6218 


2.1616 


-5.1139 


0.3901 


5.5474 


4.2426 


4.2426 


0 


0.0000 


2.8284 


-2.8284 


3.5355 


-3.5355 


8.0000 


0 


1.4142 


1.4142 


\ -5.6218 


-1.3925 


-2.1616 


-5.1139 


-5.5474 


-0.3901 



with \\R 37 \\ = 8.4542e - Oil < e, and ||X*|| = 26.9134. 
(II) The solution of Problem II 

Suppose that the given matrix Xq G GSR 7x6 is 





< -1.4142 


-1.4142 


0 


- 1.0000 


-2.8284 


2.8284 \ 






2.5000 


2.5000 


-3.5355 


-4.2426 


-2.5000 


-6.4999 






-2.5000 


-2.5000 


-3.5355 


4.2426 


-6.4999 


-2.5000 




*o = 


-0.5000 


-7.4999 


6.3639 


-3.5355 


- 1.0000 


6.9999 






0 


0 


0 


9.0000 


-1.4142 


1.4142 






4.2426 


-4.2426 


7.0000 


0 


2.1213 


2.1213 






V -7.4999 


-0.5000 


-6.3639 


-3.5355 


-6.9999 


1.0000 / 




Compute Cq = 


AX 0 B, let X = X 


- X 0 , C 


= C-C 0 


then we 


can obtain the optimal 


approximation solution X of Problem II by finding 


the least-norm solution X* of equation (10), 


that is, 


( 7.7781 


7.7781 


0 


5.0000 


4.9497 


-4.9497 ^ 






-2.8110 


5.7857 


-0.6363 


7.4894 


2.1856 


-0.2953 






-5.7857 


2.8110 


-0.6363 


-7.4894 


-0.2953 


2.1856 




X* = 


-0.8110 


1.7858 


-4.1718 


-1.7029 


1.6856 


-1.7953 


7 




4.2426 


4.2426 


0 


-9.0000 


4.2426 


-4.2426 






-0.7071 


0.7071 


1.0000 


0 


-0.7071 


-0.7071 






\ 1.7858 


-0.8110 


4.1718 


-1.7029 


1.7953 


-1.6856 ) 





Hence the solution of Problem II is 
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/ 6.3639 


6.3639 


0 


4.0000 


2.1213 


-2.1213 \ 


-0.3110 


8.2857 


-4.1718 


3.2468 


-0.3144 


-6.7952 


-8.2857 


0.3110 


-4.1718 


-3.2468 


-6.7952 


-0.3144 


-1.3110 


-5.7141 


2.1921 


-5.2384 


0.6856 


5.2046 


4.2426 


4.2426 


0 


0 


2.8284 


-2.8284 


3.5355 


-3.5355 


8.0000 


0 


1.4142 


1.4142 


\ -5.7141 


-1.3110 


-2.1921 


-5.2384 


-5.2046 


-0.6856 / 



Example 2. Given matrices A, B , C being the same as example 2 in [8], 



A = 



and 



1-10 3 

-1 -3 -4 4 

3 14 2 



, B = 



( 2 -10 1 \ 
-3 0 1-1 

0-241 

V 1 -2 1 1 / 



, c = 



20 3 -22 

24 24 -72 

16 -18 28 





( l 


0 


0 


0 > 




( 0 


1 


0 


0 > 




0 


0 


0 


-1 




1 


0 


0 


0 


R = 


0 


0 


1 


0 


, 5 = 


0 


0 


-1 


0 



\ 0 -1 0 0 / 



\ 0 0 0 \ ) 



2 ^ 
6 

- 2 



Choosing positive number e = l.Oe — 005, we will stop the iteration when ||Pfc|| < £, or 
||i?fc|| > e, but ||P/c|| < £. Let initial iterative matrix X\ = 0 £ P 4x4 , by the iterative method, 
we obtain that ||i?6|| = 2.1149e + 003 > £, and ||P 6 || = 7.4838e — 009 < e. Then, Theorem 2.1 
implies that Problem I has no solution in GSR 4xA . 
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On an equation of the Smarandache function 

Lixiang Cai 

Department of Mathematics, Northwest University, Xi’an, Shaanxi, P.R. China 

Abstract For any positive integer n, let 5(n) denotes the Smarandache function. The main 
purpose of this paper is using the elementary method to study the solvability of the equation 
S(a) 2 + S(b)~ 2 = S(c)~ 2 , and obtain its all positive integer solutions. 

Keywords The Smarandache function, Pythagorean triple, equation, positive integer solut- 
ions. 



§1. Introduction and result 

For any positive integer n, the famous Smarandache function 5(n) is defined as the smallest 
positive integer m such that in! is divisible by n. That is, 

S(n) = min{TO : n | in!, m £ N}, 



where N denotes the set of all positive integers. For example, the first few values of S(n) are: 
5(1) = 1, 5(2) = 2, 5(3) = 3, 5(4) = 4, 5(5) = 5, 5(6) = 3, 5(7) = 7, 5(8) = 4, 5(9) = 6, 
5(10) = 15, 5(11) = 11, 5(12) = 4, 5(13) = 13, 5(14) = 7, 5(15) = 5, 5(16) = 6, 5(17) = 17, 
5(18) = 6, 5(19) = 19, 5(20) = 5, •••. From the definition of 5(n) we can easily get the 
following conclusions: if n = p a , S(n) = ap, a < p; 5(n) < ap, a > p. Let n = p^p^ 2 • • •p“ r 
denotes the factorization of n into prime powers, then S(n) = max {5(p“’)|. 

l<i<r 

About other properties of 5(n), many scholars had studied it, and obtained some interesting 
results, see [1], [2], [3], [4] and [5]. For example, in reference [5], Xu Zhefeng proved that for 
any real number x > 1, 



£(5(") - PH) 2 

n<x 



2C(|)* § 

3 In a; 



+ o 




where P{n) is the greatest prime divisor of n, £(s) is the Riemann zeta-function. 

In reference [6], Rongi Chen and Maohua Le proved that the equation 5 2 (n) + S(n) = kn 
has infinitely many positive integer solutions for every positive integer k. 

On the other hand, Bencze [7] proposed the following problem: Find all positive integer 
solutions of the equation 



111 



5(a) 2 S{b) 2 5(c) 2 ' 



(1) 
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About this problem, it seems that none had studied it yet, at least we have not seen any 

related papers before. In this paper, we use the elementary method to study this problem, and 

solved it completely. That is, we shall prove the following: 

Theorem. The equation (1) has infinite positive integer solutions. And each (a, 6, c) 

satisfy the equation (1) if and only if S 2 (a ) + S 2 (b) = z 2 , which means 5(a), 5(6), z are 

Pythagorean triples, and GCD(S(a ), 5(6)) = d > 1, z | d 2 , 5(c) = where GCD(S(a), 5(6)) 

\ z \ 

denotes the greatest common divisor of 5(a) and 5(6). 

From this Theorem we can obtain many positive integer solutions of the equation (1). For 
example, if we take = 3, =4, ^ = 5, d = 5, 5(c) = 12, then we can get a = 5 3 , 6 = 5 4 , 

c = 3 5 . If d = 10, then a = 5 7 , 6 = 5 9 , c = 3 10 . 



§2. Proof of the theorem 



In this section, we shall complete the proof of our theorem. First we need the following: 
Lemma. For every positive integer to, the equation 5(n) = to has positive integer solution. 
Proof. From the definition of 5(n) and reference [8] we can easily get this conclusion. 
Now we use this Lemma to prove our theorem. For convenience, we denote 5(a) by x, 5(6) 
by y, 5(c) by w. Then the equation (1) become 



1 1 _ 1 
x 2 y 2 w 2 



(2) 



It is clear that if one of 5(a), s(b) and 5(c) is 1, then (a, 6, ) does not satisfy the equation 
(1). So without loss of generality we can assume that x > 2, y > 2, w > 2. Obviously the 
solvability of the equation (1) is equivalent to the solvability of the equation (2). Now we 
consider the positive integer solutions of the equation (2). From equation (2) we can obtain 



2 2 / 2 I 2 \ 2 

x y = (x + y )w 



( 3 ) 



and 



x 2 y 2 
c 2 + y 2 



( 4 ) 



Suppose that the equation (1) has positive integer solutions, then there exist x, y, w 
satisfied the equation (3). Noting that x > 2, y > 2, w > 2, form the equation (3) we can get 
that x 2 + y 2 is a complete square, in other words, there exists an positive integer z such that 
x 2 + y 2 = z 2 , and so x, y, z are Pythagorean triples. Hence we can get 



(xyy 2 

ItJ =w - 



( 5 ) 



Suppose GCD(x, y) = d, we can get GCD(x, y) = GCD(x , z) = GCD(y , z) = d in 
equation (5). Suppose x = dx' , y = dy ' , z = dz\ then we have 

(dx’ij 



(6) 
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from equation (6), w is integer, then zJ \ dx'y'. Note that GCD{x ', z') = GCD{y ' , z') = 1 and 

GCD(x'y', z') = 1, then we can get z' \ d and d > z' > 2. Because z = dz' , z \ d 2 , we can get 
xy 

w = — from equation (5). 

z 

At last, we get the positive integer solutions of equation (2), and the solutions satisfying 

x 2 + y 2 is a complete square, in other words, there exists a positive integer z such that x 2 + y 2 = 

xy 

z 2 , and GCD(x, y) = d > 2, z \ d 2 , ui = — . From Lemma we can get the positive integer 

z 

solutions of equation (2). Since there are infinite Pythagorean triples x, y, z satisfied the 
above situation, so the equation (2) has infinite positive integer solutions. And the solvability 
of equation (1) is equivalent to the solvability equation (2), so the equation (1) also has infinite 
positive integer solutions. This completes the proof of Theorem. 
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Abstract In this paper, we used the elementary method to study the divisibility properties 

of the integer, and improved the Smarandache divisibility theorem. 
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§1. Introduction and result 



Divisibility theory is the fundamental part in mathematics. The study of the divisibility 
theory is very important and interesting contents in number theory, many authors pay much at- 
tention to this topic, and some deeply results have been obtained. For example, F. Smarandache 
[1] proved that for any positive integer to, we have 

to | (a m — a) • (to — 1)!. 



Maohua Le [3] improved the above result, and obtained the following conclusion: 



to | (a m - a) ■ 



TO 

L 2"! 



Because (to, to — 1) = 1, it is obvious that to | ( a m — a) ■ (to — 2)!. In this paper, we using the 
elementary method to study this problem, and give another form of the Smarandache divisibility 
theorem. That is, we shall prove the following conclusion: 

Theorem. If a and to are two integers with to > 1, then we have 



( a m — a) • (to — 2)! 





'TO — 1' 




4 


4 





where [x] denotes the largest integer less than or equal to x. 



§2. Some preliminary lemmas 



In this section, we shall give several simple Lemmas which are necessary in the proof of 
our theorem. They are stated as follows: 

Lemma 1 . If a and to are two integers with m > 0, then ( a m — a)(m — 2)! can be divided 
by to. 

Proof. (See reference [2]). 

Lemma 2. Let a and b are two positive integers, S = {1, 2, 3, ...., a}. Then the number 
of all positive integers in S which can be divided by b is — ^ 
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Proof. (See reference [1]). 

Lemma 3. Little Fermat’s Theorem: For any integer a and any prime p , we have the 
congruence 

a p = a mod p. 

Proof. (See reference [3]). 



§3. Proof of the theorems 



In this section, we shall complete the proof of our theorem directly. It is clear that we just 
need to prove that 



m — 1 



( a m — a) ■ (to — 2)! 



to 



(1) 



Now we separate to into three cases: 

(a). If to = 1, 2, 3, 4, then from Lemma 1 we know that (1) is true. 











' TO — 1 


5, then 


TO — 1 

4 


= i. 


4 


4 


(a m - a) 


• (to — 


2)! 


CO 

"cT 

1 

m 

3 




m 






5 



We also know that there exists at least an even number among a, (a + 1) and (a — 1). So 
we must have 4 | (a 5 — a) • (5 — 2)!, note that (4, 5) = 1, we may immediately deduce that (1) 
is correct. 



If to = 8, then 



to — 1 



= 1. 4 



TO — 1 



= 4. 



( a m — a) ■ (to — 2)! (a 8 — a) • 6! a ■ ( a 7 — 1) • 6! 



It is clear that there exists an even number between a and ( a 1 — 1), so we have 2 | (a 8 — a), 
4 2 | 6!. That is to say, our Theorem is true if to = 8. 

(c). If to > 5, and to yf 8, then we can discuss it in two cases: 

If to is an odd number, we can write m = Ak + f, where t = 1 or 3, k G N. Then 



4L 



'TO — 1' 




Ak + t — 1 


4 


= 4 


4 



= A k . 



If to = Ak + 1, then we have the estimate 

OO 

E 



to — 2 




TO — 2 


+ 


TO — 2 


2? 




2 


4 



= (2k — 1) + (fc — 1) + • • • > (3 k - 2) > 2k. 



(2) 
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If m = 4k + 3, then we also have 



OO 



E 



m — 2 




to — 2 


+ 


TO — 2 


23 




2 


4 



= 2k + k +■•■> 3k > 2k. 



( 3 ) 



From (2) and (3) we may immediately deduce that 4 k is a divisor of (4 k — 2)1. Let (to — 2)! = 
4 k • D, where D is a positive integer. Note that (4k + t, 4 k ) = 1, so (4 k + t ) | ( a m — a) ■ D. 
Therefore, if m > 5 is an odd number, then (1) is correct. 

If to is an even number, we can write m = 2 r ■ q, where q is an odd number. It is clear that 
if r = 1, then from the above we know that the (1) is true. So we can assume that r >2. This 
time, in order to prove our theorem, we just need to prove that 



Similarly, we also have 



2 2 ^.q- 2 +r | (a 2 '’ q -a)(2 r g-2)\ 

q 



oo 



E 



TO — 2 




TO — 2 


+ 


TO — 2 


+ • 


• + 


m — 2 


2i 




2 


4 


2 r 









= 2 r ~ 1 q- l + 2 r ~ 2 q-l + ■■■ + q- 1 + ■ ■ • 
> 2 r ~ 1 q — 1 + 2 r ~ 2 q — 1 + • • • + g — 1 
= q(2 r ~ 1 + 2 r ~ 2 + ■ ■ ■ + 1) - r 

= q2 r - q - r. 



( 4 ) 



Note that r > 2 and (2, q) = 1, so we have (q-2 r — q— r) — (2 ( ' r ~ 1 ^ -q—2+r) = 2 r ~ 1 q—q+2—2r > 0. 
Thus, 2 2 9 - 2 +r j g an divisor of (to — 2)!. Since q < to — 2, we also have q | (to — 2)!. Therefore, 

from (2), (3) and (4) we know that (1) is correct, if to > 5 ( m ^ 8 ) is an even number. 
Combining (a), (b) and (c) we may immediately complete the proof of our Theorem. 
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§1. Introduction and result 

For any positive integer n, the F.Smarandache permutation sequence {P(n)} is defined 
as P(n) = 135 • • • (2 n — l)(2n) • • • 42. For example, the first few value of the sequence {P(n)} 

are: P(l) = 12, P( 2) = 1342, P( 3) = 135642, P(4) = 13578642, This sequence 

was introduced by professor F.Smarandache in reference [1], where he asked us to study its 
elementary properties. About this problem, many people had studied it, and obtained a series 
valuable results, see references [2], [3] and [4], In reference [5], F.Smarandache proposed the 
following problem: Is there any perfect power among the permutation sequences? That is, 
whether there exist positive integers n, m and k with k > 2 such that P(n) = m k . 

The main purpose of the paper is using the elementary method to study this problem, 
and solved it completely. At the same time, we also obtained some other properties of the 
permutation sequence {P(n)}. That is, we shall prove the following: 

Theorem 1. There is no any perfect power among the permutation sequence. 

Theorem 2. Among the permutation sequence {P(n)}, there do not exist the number 
which have the factor 2 k , where k > 2; There exist infinite positive integers a £ {P(n)} such 
that 3 2 divide a. 



§2. Proof of the theorems 

In this section, we shall use the elementary method to prove our Theorems directly. First 
we prove Theorem 1. For any positive integer n > 2, note that P(n) be an even number and 

P(n) = 135 • • • (2 n - 1)(2 n) • • • 42 

= 10 ai + 3 x 10“ 2 + • • • + (2n - 1) x 10“" + (2 n) x 10“" +1 + • • • + 4 x 10 + 2. (1) 

So from (1) we know that 

P(n) = 10 ai + 3 x 10“ 2 + • • • + (2 n - 1) x 10“" + (2n) x 10“" +1 + -- -+ 4xl0 + 2 = 0 mod 2 
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and 



P(n) = 10 Q1 +3 x 10“ 2 + • •• + (2n- 1) x 10“" + (2n) x 10“" +1 + • •• +4 x 10 + 2 

= 2 mod 4. (2) 

That is to say, 2 | P(n) and 4f P(n). So P(n) is not a perfect power. Otherwise, we can write 
P(n) = m k , where k and m > 2. Since P(n) be an even number, so to must be an even number. 
Therefore, P(n) = m k = 0 mod 4. Contradiction with 4 f P(n). This proves Theorem 1. 

Now we prove Theorem 2. For any positive integer n > 2, from (2) we may immediately 
deduce that 4 r f P(n), where r > 2 is an integer. 

Now we write P(n) as: 

P(n) = 10 ai + 3 x 10“ 2 + • • • + (2n - 1) x 10“" + (2 n) x 10“" +1 + • • • + 4 x 10 + 2. (3) 

Note that 10 = 1 mod 9, from the properties of the congruence we may get 10* = 1 mod 9 for 
all integer t > 1. Therefore, 

P(n) = 10 Q1 + 3 x 10“ 2 + ■ • • + (2n - 1) x 10“" + (2n) x 10“" +1 + • • • + 4 x 10 + 2 

= 1 + 2 + 3 + -- - + -- - + 2 n = n • (2 7Z -t- 1) mod 9. (4) 

From (4) we know that P(n) = 0 mod 9, if n = 0 or 4 mod 9. That is to say, for all integers 
n > 2, P(9n) and P(9n + 4) can be divided by 9 = 3 2 . This completes the proof of Theorems. 
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Abstract For any positive integer n > 1, the Smarandache quotients Q(n) is defined as the 
smallest positive integer k such that n ■ k is a factorial number. That is, Q(n) = min{fc : 
n ■ k = m!}. The main purpose of this paper is using the elementary method to study 
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§1. Introduction and results 

In his book “Only problems, not solutions”, professor F. Smarandache introduced many 
functions, sequences and unsolved problems, many authors had studied it, see references [1], 
[2] and [3]. One of the unsolved problems is the Smarandache quotients sequence {Q(n)}, it 
is defined as the smallest positive integer k such that n ■ k is a factorial number. That is, 
Q(n) = min {A; : n • k = ml}, where m is a positive integer. For example, from the definition of 
Q(n) we can find that the first few values of Q(n) are Q{ 1) = 1, Q{ 2) = 1, Q( 3) = 2, Q{ 4) = 6, 
Q(5) = 24, Q(6) = 1, Q(7) = 720, Q( 8) = 3, Q(9) = 80, Q(10) = 12, Q(ll) = 3628800, 
Q(12) = 2, Q(13) = 479001600, Q(14) = 360, Q(15) = 8, Q(16) =45, 

In reference [4], professor F. Smarandache asked us to study the properties of the sequence 
{<3(n)}. About this problem, some authors had studied it, and obtained several simple results. 
For example, Kenichiro Kashihara [5] proved that for any prime p, we have Q(p) = (p — 1)!. 

In this paper, we use the elementary method to study the properties of an infinite series 
involving the Smarandache quotients sequence, and give an interesting identity. That is, we 
shall prove the following conclusion: 

Theorem 1. Let d(n ) denotes the Dirichlet divisor function, then we have the identity 

+00 +00 J/ |\ 

E l 171 ' d( m -) 

Q(n ) • n ^ (m + 1)! 

n— 1 v 7 TTi—1 v 7 

Let <3i(2n — 1) denotes the smallest positive odd number such that Q\(2n — 1) • (2 n — 1) 
is a two factorial number. That is, Qi(2n — 1) • (2 n — 1) = (2m — 1)!!, where (2m — 1)!! = 
1 x 3 x 5 x • • • x (2m. — 1). Then for the sequence {Qi(2n — 1)}, we can also get the following: 
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Theorem 2. Let d(n) denote the Dirichlet divisor function, then we also have the identity 

+oo 1 +oo 



E 



i 



Qi(2n — 1) • (2n — 1) 






to • d((2m — 1)!!) 
( 2 to + 1 )!! 



§2. Proof of the theorems 

In this section, we shall use the elementary method to complete the proof of our theorems 
directly. For any positive integer n, let Q(n) = k, from the definition of Q(n) we know that 
k is the smallest positive integer such that n ■ k = to!. So k must be a divisor of to!. This 
implies that k | to!, and the number of all k (such that k | to!) is d{m\). On the other hand, if 
k | (to — 1)!, then Q(n) ^ k , and the number of all k (such that k \ (to — 1)!) is d((m — 1)!). So 
the number of all k (such that k \ ml and k f (to — 1)!) is d(m\) — d((m — 1)!). That is means, 
the number of all k (such that Q{n) = k and Q[n ) • n = to! ) is d(m\) — d((m — 1)!). From this 
we may immediately get 



+oo 1 

V 


+oo 

V 


+oo 

V 


l 


+oo 

v 


#{Q(n) • n 


= to!} 




h ■ n 


m= 1 


n= 1 


to! 


Z-^ 

m—1 




to! 










Q{n)-n=m\ 
















= i + 


d{m\ 

m — 2 


)-<*((' 

ml 


m — 


1)!) 


+oo 

= E 

m—1 


d(m\) 

ml 


+oo 

E 

m—1 




+oo 

V 


(to + 1) • 


d(m\) 


~d( 


to!) 


+oo 

v 


to • d(ml) 




m= 1 


0 


m + 1) 






m—1 


(to- + 1)! 


5 



where d(m) is the Dirichlet divisor function, and #{Q(n) • n = to!} denotes the number of all 
solutions of the equation Q(n ) • n = ml. This proves Theorem 1. 

Similarly, we can also deduce Theorem 2. This completes the proof of Theorems. 
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§1. Introduction 

The Smarandache bisymmetric determinant natural sequence (SBDNS), introduced by 
Murthy [1], is defined as follows. 

Definition 1.1. The Smarandache bisymmetric determinant natural sequence, {SBDNS(n)}, 
is 



1 2 
2 1 



12 3 

2 3 2 

3 2 1 



12 3 4 

2 3 4 3 

3 4 3 2 

4 3 2 1 



1 






A first few terms of the sequence are 1, -3, -8, 20, 48, -112, -256, 576, 

The following result is due to Majumdar [2]. 

Theorem 1.1. Let a n be the n-th term of the Smarandache bisymmetric determinant 
natural sequence. Then, 



1 
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n — 2 


n — 1 


n 
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n — 1 
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n — 1 


rn -1 
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n — 1 


n — 2 


= (— 1) L 2 J (n+ 1) 2 n ~ 2 
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Let {S n } be the sequence of n — th partial sums of the sequence {a n }, so that 



n 




This paper gives explicit expressions for the sequence {SVi}- This is given in Theorem 3.1 in 
Section 3. In Section 2, we give some preliminary results that would be necessary for the proof 
of the theorem. We conclude this paper with some remarks in the final section, Section 4. 

§2. Some preliminary results 

In this section, we derive some preliminary results that would be necessary in deriving the 
expressions of S n in the next section. These are given in the following two lemmas. 

Lemma 2.1. For any integer m > 1, 





fc=l,3,...,(2m-l) 



Proof. Since 





fc=l,3,...,(2m— 1) 



is a geometric series with common ratio 2 4 , the result follows. 
Lemma 2.2. For any integer m > 1, 




fc=l,3,...,(2m— 1) 



Proof. Denoting by s the series on the left above, we see that 



s = 1 + 3 • 2 4 + 5 • 2 8 + • • • + (2m - 1) • 2 4(m ~ 1) , 
so that, multiplying throughout by 2 4 , we get 



2 4 s = 1 • 2 4 + 3 • 2 8 + • • • + (2m - 3) • 2 4(m " 1) + (2m - 1) • 2 4m • 



Now, subtracting (**) from (*), we have 



(1 - 2 4 )s = 1 + 2 • 2 4 [l + 2 4 + • • • + 2 4(m " 2) l - (2m - 1) • 2 4m 




— (2 4m+1 - 17) - (2m - 1) • 2 4m 
15 



which now gives the desired result. 
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§3. Main results 

In this section, we derive the explicit expressions of the n-th partial sums, S n , of the 
Smarandache bisymmetric determinant natural sequence. 

From Theorem 1.1, we see that, for any integer k > 1, 

+ fl2fe+i = ( — l) fe (6fc + 5)2 2 *- A; 1 \ 



0’2k+2 + fl2fe+3 — (— l) fc+1 (6fc + 11)2 



-\2k 



so that 



«2 k 



+ d2fc+l + fl2fc+2 + «2/c+3 — 3(— l) fc+1 (6fc + 13)2 2 ^ fe 



(1) 



Letting 



we can prove the following result. 



S n — ®1 + 02 + ■ • • + cin 



Theorem 3.1. For any integer m > 0, 

1) . S 4m+ i =\-m- 2 2 ( 2m+1 ) + § • 2 4m - ^ ^ {(60m + 31)2 4m " 1 - 3}; 

2) . S im+ 2 = 2 4m+3 - H . 2 2 l 2m+1 ) - A = - A {(10m + ll)2 4m+1 + 3}; 

5 25 Zb Zb 

3) . S Am+ 3 = • m • 2 4 ^ +1 ) - H • 2 2 ( 2m+1 ) - ^ ^ A {(60m + 61)2 4m+1 + 3}; 

4) . S 4m+4 = l ■ m • 2 4m+5 + 1 • 2 4 ( 2m + 2 ) - A = 1 | (5m + 8)2 4 ( m +D - 3}. 

Proof. To prove the theorem, we make use of Lemma 2.1 and Lemma 2.2, as well as 
Theorem 1.1. 

1). Since can be written as 

<54771+ 1 = Ol + 02 + • • • + 0 4m+1 = Oi + (a2fc + 02fc+l + 02fc+2 + O 2 /C+ 3 ); 

fe=l,3,...,(2m— 1) 



by virtue of (1), 



5 4 m+ i=oi + 3 Y (-l) fc+1 (6fc + 13)2 2 ( fe - 1 ) 
fc=l,3,...,(2m-l) 

= Oi + 3 | 6 Y fc2 2(fe_1) + 13 Y 2 2(fe_1) 

I fe=l,3,...,(2m— 1) fc=l,3,...,(2m-l) 



Now, appealing to Lemma 2.1 and Lemma 2.2, we get 

S 4m+1 = 1 + [^7 (2m - l)2 4m - 4( 24m+1 - 17 ) 
5 25 



+ ^(2 4 m -l) 
5 
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which now gives the desired result after some algebraic manipulations. 

2). Since 



S 4m+2 — \ I 4“ &4m+2j 



from part (1) above, together with Theorem 1.1, we get 



’4m+2 — 



- • m 
5 



22(2m+l) | 9 4m 

25 ' 25 



- (4m + 3) • 2 4m , 



which gives the desired expression for Si m+ 2 after algebraic simplifications. 

3). Since 



*54m+3 — ^4m+2 4” d4 m _|_3 — 



_i 2 4m + 3 - 2 2 ( 2m+1 ) 

5 25 ’ 25 



- (4?n + 4) • 2 4m+1 , 



we get the desired expression for .S' 4 m+ 3 after simplifications. 

4). Since 



*$4771+4 — *$4t 77+3 4” 4 4m+4 — 



-- m o 4m+3 - — 2 2 ( 2m+1 ^ — 

5 25 25 



- (4?n 4- 5) • 2 4m+2 , 



the result follows after some algebraic simplifications. 
The case when m = 0 can easily be verified. 
Hence, the proof is complete. 



§4. Remarks 

Theorem 3.1 in the previous section gives the n-th term of the sequence of partial sums, 
{S'„}, in all the possible four cases. The following lemmas prove that, in each case, S n is indeed 
an integer. 

Lemma 4.1. For any integer m > 0, 2 4m_1 (60m 4- 31) — 3 is divisible by 25. 

Proof. The result is true for rn = 0, 1. So, we assume its validity for some positive integer 
m. Now, since 

[2 4m+3 {60(m + 1) + 31} - 3] - {2 4m - 1 (60m + 31) - 3} = 25(36 to + 57)2 4m " 1 , 

it follows, by virtue of the induction hypothesis, that 2 4m+3 {60(m4- 1) 4- 31} — 3 is also divisible 
by 25. Thus, the result is true for m + 1 as well, completing induction. 

Lemma 4.2. For any integer m > 0, 2 4m+1 (10?n 4- 11) 4- 3 is divisible by 25. 

Proof. is by induction on m. The result is clearly true for m = 0,1. Now, assuming its 
validity for some positive integer m, since 

[2 4m+5 {10(m + 1) + 11} + 3] - {2 4m+1 (10m + 11) + 3} = 25(6m + 13)2 4m+1 , 

is divisible by 25, it follows that the result is true for m + 1 as well. This completes the proof. 
Lemma 4.3. For any integer to > 0, 2 4m+1 (60m 4- 61) 4- 3 is divisible by 25. 
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Proof. is by induction on to. The result is clearly true for to = 0, 1. We now assume that 
the result is true for some positive integer to. Then, since 

[2 4m+5 {60(?n + 1) + 61} + 3] - {2 4m+1 (60m, + 61) + 3} = 25(36 to + 75)2 4m+1 , 



this, together with the induction hypothesis, shows that the result is true for m + 1 as well. 
This completes the proof by induction. 

Lemma 4.4. For any integer to > 0, 2 4 l m+1 )(5TO + 32) — 3 is divisible by 25. 

Proof. We first assume that the result is true for some positive integer to. Now, since 

[ 2 4( m +2){ 5 ( m + x) + 13 | _ 3 ] _ {2 4 ( m+1 )(5TO + 13) - 3} = 25(3m + 8)2 4(m+1) , 



this together with the induction hypothesis, shows that 2 4 ( m+2) {5(TO + 1) + 13} — 3 is also 
divisible by 25. This, in turn, shows that the result is true for to + 1 as well. To complete, we 
have to prove the validity of the result for to = 0, 1, which can easily be checked. 

The Smarandache bisymmetric arithmetic determinant sequence, introduced by Murthy 
[!], is 





a 


a + d 


a, 








a + d 


a 



a 


a + d 


cl + 2d 


a + d 


a T 2d 


a + d 


cl + 2 d 


a + d 


a 



The n-th term of the above sequence has been found by Majumdar [2] to be 



(- 1 ) 




(2d) n ~ 1 . 



Note that the Smarandache bisymmetric determinant natural sequence is a particular case of 
the Smarandache bisymmetric arithmetic determinant sequence when a = 1 and d = 1. 

Open Problem: To find a formula for the sum of the first n terms of the Smarandache 
bisymmetric arithmetic determinant sequence. 
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§1. Introduction and results 



For any fixed natural number fc > 2 and any positive integer n, we call b k (n) as a fc-power 
complements function of n, if b k {n) denotes the smallest positive integer such that n • bk(n) is 
a perfect fc-power. Especially, we call b 2 (n), 63 (n), b±[n) as the square complements function, 
cubic complements function and quartic complements function respectively. In reference [1], 
Professor F.Smarandache asked us to study the properties of the fc-power complements function. 
About this problem, there are many authors had studied it, and obtained a series results. For 
example, in reference [2], Wenpeng Zhang calculated the value of the series 



+ OO 

E 



71=1 



1 

(n • b k (n)) s ’ 



where s is a complex number with Re(a) > 1, fc = 2, 3, 4, and prove that 



+ OO 

E 



71=1 



1 

(n • b 2 (n)) s 



c 2 ( 23 ) 
C(4s) ’ 



where £(s) is the Riemann-zeta function. 

Ruiqin Fu [3] discussed the asymptotic properties of In 62 f ci! ) , and proved that 

( f — A In 5 77, 

n • exp r 

\ (In In n) S 

where n > 2 is an positive integer, A > 0 is a constant. 

But for the properties of the fc-power complements function of n, we still know very little 
at present. In this paper, we use the elementary method to study the asymptotic properties 
of &fc(n!), and give an interesting asymptotic formula for ln 6 fc(n!). That is, we shall prove the 
following general conclusion: 
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Theorem. For any fixed natural number k > 2 and any positive integer n, we have the 
asymptotic formula 



ln&fc(n!) = n I k — 



1 



d{kd + 1) 



+ O n ■ exp 



—Ain 5 n 
(In In n) s 



where A > 0 is a constant, and exp(y) = e v . 



§2. Two simple lemmas 



To complete the proof of the theorem, we need the following two simple lemmas: 

Lemma 1. Let n\ = p^p^ 2 ■ • -p“ s denotes the factorization of n\ into prime powers. Then 
we have the calculate formula 



b k (n\) = 



where the orcl function is defined as: 



ord(pi) = < 



h{pT)bk{P2 2 )---h{p^) 

ord(p i) ord(p 2 ) ord(p a ) 

Pi 'Pi '"Pa 

k — 1, ifctj = km + 1, 

k — 2, ifctj = km + 2, 

0, ifcti = k{m + 1). 



where m = 0, 1, 2, • • • . 

Proof. See reference [ 3 ]. 

Lemma 2. For any real number x > 2 , we have the asymptotic formula 



9{x) = E ln P = 

p<x 



1 + O \ x ■ exp 



—Tin 5 x 
(In In x) i 



where exp(y) = e v , and A > 0 is a constant. 
Proof. See reference [ 4 ] and [ 5 ]. 



§3. Proof of the theorem 

In this section, we shall complete the proof of Theorem. First from Lemma 1 we have 

In b k (n\) = ln(6 fc (p“ 1 )6fe(p2 2 )-"^fc(p“ a )) 

= E] ord{p) In p 

p<n 

k]ord(p) 

= (k — 1) E lnp + (k — 2) E In p + • • • + E lnp 

f<p<n f<P<f 

+ (fe-i) E In p-\ h ^2 hO(l). 

fcT2 <P— fcTl 2k — 1 
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Let n be a positive large enough, if a prime factor p of n\ in the interval ( — , n] , then the power 
of p is 1 in the standard factorization of n! . Similarly, if a prime factor p of n! in the interval 

Tl Tl 

( — , — ], then the power of p is 2 in the standard factorization of n!; If a prime factor p of n\ in 

~ Tl Tl 

the interval ( — , — ], then the power of p is 3 in the standard factorization of n\, ■ ■ ■ . 

Then from the above formula and Lemma 2 we have 



In bk{n\) = (k - 1) (o(n) - 9 (^)) + (k - 2) (o ( 
+ (*- 1 ) 



+ 



k - 



+ o 



lim n 

M—too 



lim n 

M — >-oo 



= n 




S d ( u+i ) 



where A > 0 is a constant, and exp(y) = e v . 

This completes the proof of Theorem. Absolutely, this theorem can be used to study 
k- power residues numbers. 
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Abstract This paper gives an alternative approach to the solution of the LP problem with 
equality constraints. In this scheme, the constraint system of linear equations is first solved 
by elementary row operations (if such a solution exists). In the second stage, the solution 
of the system of linear equations is directly plugged in the objective function and then it is 
minimized. This scheme thus avoids the use of the simplex method. 

Keywords Linear programming, equality constraint, elementary row operations. 



§1. Introduction 



In a recent paper, Ru, Shen and Xue [1] considered the problem of finding an initial basic 
feasible solution (bfs) of the LP problem of the form 



min z = C l X s.t. AX = B, X > 0, 



(1) 



where 



C = 





an a 12 ■ • • a i n 






u \ 




021 022 • • • 02 n 




U \ 




( xi \ 


C2 








b 2 




X2 




,A = 




,B = 




,x = 








Oil 0^2 ‘ * * Ojin 










V c " / 








^ bm j 




\x n J 



y Oml ®m2 ‘ ‘ ‘ O mn J 



In (1), without loss of generality, we may assume that B > 0; for otherwise, if some < 0 
(1 < * < to), then multiplying the i-th constraint throughout by (-1), we get one with (—hi) > 0. 

To find an initial bfs of the LP problem (1), the procedure suggested by Ru et al. [1] is as 
follows : 

Step 1 : Starting with the augmented matrix ( A \ B) : reduce it to a form with an identity 
matrix of desired order, by elementary row operations only. 

Note that, the elementary row operations applicable in this case are 
(a) multiplying any row of ( A \ B) by a positive constant, 



90 



A.A.K. Majumdar 



No. 3 



(b) multiplying any row of ( A | B) by a constant and adding it to a second row. 

Throughout this paper, the first row operation would be denoted by Ri — > kRi (which 
means that the i-th row Ri is multiplied by the constant k , where k > 0), and the second row 
operation would be denoted by Ri — > Ri + kRj (which means that the j- th row Rj is multiplied 
by the constant k and the resulting row is added, term-by-term, to the i-tli row R,, to get the 
new i-th row). Note that the latter row operation is allowed only if bi + kbj > 0. 

If rank(A) = rank (A \ B) = l{< to), then using the above elementary row operations, it 
is possible to find a unit matrix of order l , which may be exploited to find an initial bfs of the 
LP problem (1). It may be mentioned here that, if l < to, then only l of the to system of linear 
equations AX = B are independent, and the remaining l — m equations can be ignored. It may 
also be mentioned here that, if rank(A) ^ rank{A \ B), then the system of linear equations 
AX = B is inconsistent, and as such, the LP problem (1) is infeasible. 

Step 2 : With the initial bfs found in Step 1 above, form the initial simplex tableau in the 
usual manner. 

Step 3 : Proceed in the usual way of the simplex method to find an optimal solution. 

The method outlined above is nevertheless not a new one. In fact, the same procedure has 
been followed by Papadimitriou and Steiglitz [2] . 

In the next section, we propose an alternative approach. 



§2. An alternative approach 

In this section, we suggest an alternative approach to the solution of the LP problem (1). 
The procedure is quite simple, and can be done in the following two steps. 

Step 1 : Reduce the system of linear equations AX = B to a form containing the unit 
matrix of requisite order, using elementary row operations only. 

Step 2 : With the solution found in Step 1, consider the problem of minimizing z = C f X . 

Thus, the method works with the solution of the system of linear equations AX = B and 
the objective function only, avoiding the simplex method completely. 

We illustrate our method with the help of the following examples, due to Ru et al. [1]. 

Example 2.1. Consider the following LP problem : 

min z = 4xi + 3x2 s.t. 

1 2 

+X 2 +xX3 --X 4 =2 

f , 

2, X3 =3 

—6x2 + 4 X 4 = 0 

Xi,X 2 ,X 3 ,X4 > o. 

To solve the above LP problem, we first solve the constraint system of linear equations, 
using elementary row operations only. We thus start with the augmented matrix (A \ B ) and 



:Xi 



2 Xl 

3xi 
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proceed as follows, using the indicated row operations: 







( ® 


1 




1 

2 


2 
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2 N 
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2 




1 




1 
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2 
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\ 




(■ A\B) = 


3 

2 


0 
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r 3 —+ 


R% + 6ii!i 
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0 ) 








6 
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-> +3 
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— * R\ + R 2 
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-> /?, - |i? 2 






5 
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5 




R‘2 


— > R-2 + « R 2 


A 


















z 


















\ 6 


0 


3 


0 


12 , 


r 3 - 


-> R:i - f i?2 




0 


0 


1 


0 


0 


) 



The final tableau gives the following solution : 



Xi = 2, x 3 = 0, X 2 = 1 + -X 4 = 1 + -t; (t > 0 is any real number ), 
o o 

2 2 

and the objective function is z = 4xi + 3^2 = 8 + 3(1 + -t) = 11 + -t. 

o _ < o 

Since the objective is to minimize z, it is clear that 2 is minimized when t = 0, and the optimal 
solution of the given LP problem is thus x* = 2, x\ = 1, x 3 = 0, x% = 0; min 2 * = 11. 

Example 2.2. To solve the LP problem 



min z 


= — 3xi + x 2 + x 3 


s.t. 


Xl 


— 2X2 +X 3 +X4 


= 11 


— 4 :Xi 


+X 2 + 2 X 3 


X 5 3 


to 

H 


+X 3 


= 1 



XuX2,X3,Xl,X5 > 0 , 



we first solve the constraint system of linear equations by elementary row operations. This is 
done below. 



(A\B) 



/ 1 


-2 


1 


1 


0 






( -7 


0 


-4 


1 


2 


0 


-1 


3 


Ri — * Ri + 2R2 


-4 


1 


V -2 


0 


1 


0 


0 


1 ) 




l -2 


0 



5 

2 

1 



( 

Ri — > R\ + 5i?3 

V 



3 

0 

-2 



0 

1 

0 



0 

0 

1 



1 


-2 


17 > 


0 


-1 


3 


0 


0 


1 


1 


-2 


12 > 


0 


-1 


1 


0 


0 


1 / 



The final tableau is equivalent to the following system of linear equations 



3 xi 


+X4 —2X5 


= 12 


X2 


-X5 


= 2 


—2x\ 


+X 3 


= 1 



with the solution 



X\ = t\,x$ = t 2 ,X 2 = l+f2,£3 = l+2t\,Xi — 12— 3ti+2t2; (t± > 0,t2 > 0 are any real numbers) 
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and the objective function is z = 2 — t\ + t2- 

Since the objective is to minimize z, it is evident that t2 should be as small as possible and t\ 
should be as large as possible. Now, since X 4 = 12 — 3fi + 2 f 2 > 0, we see that, we must have 
t 2 = 0, ti = 4. Hence, the desired optimal solution is 



x * = 4 , 2:3 = 1 , 2:3 = 9, x\ — 0, x§ = 0; min z* = -3. 
Example 2.3. To solve the LP problem 



min z 


= 4xi + 32:3 




s.t. 


1 

I^ 1 


1 

+X2 


2 

- 3 X 4 


= 2 


2* 1 


+r 3 




= 3 


32q 


- 62:2 


+ 2:4 


= 0 






> 0 , 





we start with the augmented matrix (A\B) and proceed as follows, using the indicated row 
operations : 



fill 



(A\B) = 



° 4 



3 

0 3 



\ 



\ 3 -6 0 4 0 / 



f 



1 1 1 



R3 — > R3 — AR\ 



\ 



3 

„ _ „ 0 3 

2 4 

\ 0 0 0 0 0 / 



- 0 - 



r 3 

Ri 

R‘2 



R 3 + 6f?i 



R\ — -R'2 
2 ^ 
3 * 2 



Z 1 1 1 



2 ° 4 

\ 6 0 3 

( 0 1 1 



3 

0 3 



\ 



1 ° 2 



0 12 / 
3 

0 2 



\ 



\ 0 0 1 

The final tableau is equivalent to the following system of linear equations 



0 0 / 



X2 



x 3 



-2:4 



= 1 



Xi + -X3 



= 2 



with the solution 



1 12 

X3 = ti,X\ = 2 — -t\,X4 = t2, X2 = 1 — -t\ + -t2\ (ti > 0 , <2 > 0 are any real numbers ), 
and the objective function is z = 8 + t\. 

Since the objective is to minimize z, we must have t\ = 0. Hence, the desired solution is 



x\ = 2 , X2 = 1 + —t, x 3 = 0 , 2:4 = t(t > 0 is any real number); min z* = 8 . 
o 



The above solution shows that the given LP problem has infinite number of solutions. 
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§3. Some observations and remarks 



In the previous section, we have illustrated our method with the help of three simple 
examples. However, in large-scale problems, the success as well as the computational efficiency 
of the method remains to be checked out. 

Clearly, the method is useful when the correct basis has been found. This point is illustrated 
with the help of the following example, due to Papadimitriou and Steiglitz [2]. 

Example 3.1. Consider the LP problem below : 

min z = X\ + X2 + X3 + X4 + x$ s.t. 



32q + 2 iC 2 +213 — 1 

5 a;i +x 2 +£3 +£4 = 3 

2 a;i + 5 ic 2 +x 3 +2:5 = 4 

2;i,X2,X3,ai4,2;5 > 0. 

We first solve the constraint system of linear equations by elementary row operations, using the 
indicated row operations. 





/32100l\ 


Ri — * \R\ 


f 1 1 1 0 0 n 


( A\B) = 


5 1110 3 


R2 — * R2 — Ri 


1 0 i 1 0 § 




\251014y 


R3 — * R3 — 51 ?i 


l -¥ o -1 0 1 | J 



The final tableau is equivalent to the following system of linear equations 



t:Xi +x 2 

7 

f 1 

- x 1 



11 



1 

+ ^3 

I 

+ -X3 +£4 

r^3 



1 



+2:5 = 



with the solution 



X\ =t\,xz = tiiti > 0, > 0 are any real numbers ), 



131 571 3 11 3 

X2 “ 2 “ 2 h ~ 2 h ' X4 -2~2 h ~ 2 h ' X5 ~2 + Y h + 2 25 

9 3 

and the objective function is z = — + — (ti + f 2 ). 

Since the objective is to minimize z, it is evident that we must have ti = 0, f 2 = 0, giving the 
desired optimal solution 

* / * i~. * ^ ^ ^ 

Xi =0,x 2 = -,x 3 = 0,x 4 = -,x 5 = - ; mm z = 

On the other hand, the solution below 
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Ri - 


-> Ri - Ri 
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— > 1?3 — R\ 
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2 
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1 


3 
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(A\B) = 
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leads to the system of equations 

3aq +2x 2 +£3 = 1 

2 x\ —X2 +X4 = 2 

— X\ -1-3x2 TX 5 — 3 

whose solution is 



X\ =t\,X 2 = t 2 {t\ > 0,^2 > 0 are any real numbers ), 



ar 3 = 1 — 3ii — 2 1 2 , X4 = 2 - 2ti + f 2 , x 5 = 3 + ti - 3 1 2 ; 



and the objective function is z = 3 [2 — (£1 + £ 2 )]. 

Thus, the problem of finding the solution of the given LP problem reduces to the problem of 
solving the LP problem 

max t\ + t 2 s.t. 

3fi +2£ 2 < 1 

2 1\ — 1 2 + 2 
— £1 T3£2 ^ 3 

£i, t 2 > 0. 

The optimal solution of the above LP problem is ti = 0, t 2 = 

Hence, the optimal solution of the original LP problem is 



1 ^3 

x l = 0, x 2 = 7,, x 3 = 0 > x 4 = 0^5 = o 1 111111 ^ 



9 

2 ' 



The above example shows that, in the alternative approach proposed in this paper, as well 
as the method suggested by Ru et al. [1], the choice of the unit matrix (of desired order) during 
the solution of the system of linear equations AX = B (by elementary row operations) may 
save the number of iterations considerably. However, no particular criterion is specified in this 
respect. It thus still remains to check how far these methods can compete with the existing 
method. 
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Abstract Two triangles T(a, b, c) and T(a ,b,c ) are said to be pseudo Smarandache related 
if Z(a) = Z(a ), Z(b) = Z(b ), Z(c) = Z(c ), where Z(.) is the pseudo Smarandache function, 
and T(a,b,c) denotes the triangle with sides a,b and c. This paper proves the existence of 
an infinite family of pairs of dissimilar Pythagorean triangles that are pseudo Smarandache 
related. 

Keywords Smarandache function, pseudo Smarandache function, Smarandache related 
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§1. Introduction and result 

The Smarandache function, denoted by S(n), is defined as follows. 

Definition 1.1. For any integer n > 1, ( Z + being the set of all positive integers), 

S(n) = inin{m : m £ Z + , n|m!}. 

The following definition is due to Sastry [6] . 

Definition 1.2. Two triangles T(a 1 b,c ) (with sides of length a,b and c) and T(a ,b ,c ) 
(with sides of length a , b and c ) , are said to be Smarandache related if 

S(a) = S(a'),S(b) = S(b'),S(c ) = S(c). 



Definition 1.3. A triangle T(a,b,c) is said to be Pythagorean if and only if one of its 
angles is 90°. 

Thus, a triangle T(a, b, c) (with sides of length a, b and c) is Pythagorean if and only if 
a 2 + b 2 = c 2 . 

Sastry [6] raised the following question : Are there two distinct dissimilar Pythagorean 
triangles that are Smarandache related? Recall that two triangles T(a , b, c) and T(a ,b , c ) are 
similar if and only if the corresponding three sides are proportional, that is, if and only if 

a b c 

a b' c 



Otherwise, the two triangles are dissimilar. 
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The following result, due to Ashbacher [1], answers the question in the affirmative. 

Theorem 1.1. There are an infinite family of pairs of dissimilar Pythagorean triangles 
that are Smarandache related. 

The proof of Theorem 1.1 is rather simple : For any prime p > 17, the two families 
of dissimilar Pythagorean triangles T(3p, Ap, bp) and T(5p, 12p, 13p) are Smarandache related, 
since S(3p) = S(Ap) = S(5p) = S(l2p) = S^lSp) = p. 

Ashbacher [1] introduced the concept of pseudo Smarandache related triangles, defined as 
follows. 

Definition 1.4. Two triangles T(a, 6, c) and T(a ,6 , c ) are said to be pseudo Smarandache 
related if Z(a) = Z{a ), Z(b) = Z{b ), Z(c) = Z{c ). 

In Definition 1.4 above, Z(.) denotes the pseudo Smarandache function. Recall that the 
pseudo Smarandache function is defined as follows. 

Definition 1.5. For any integer n > 1, Z(n) is the smallest positive integer to such that 

TYiy 777, -)- \ J 

1 + 2 H + to = is divisible by n. Thus, 



Ashbacher [1] used a computer program to search for dissimilar pairs of Pythagorean 
triangles that are pseudo Smarandache related. He reports some of them, and conjectures 
that there are an infinite number of pairs of dissimilar Pythagorean triangles that are pseudo 
Smarandache related. 

In this paper, we prove the conjecture of Ashbacher [1] in the affirmative. This is done 
in Theorem 3.1 in Section 3. We proceed on the same line as that followed by Ashbacher. 
However, in the case of the pseudo Smarandache related triangles, the proof is a little bit more 
complicated. The intermediate results, needed for the proof of Theorem 3.1, are given in the 
next Section 2. 

§2. Some preliminary results 

The following result, giving the explicit expressions of Z(3p), Z(Ap) and Z(5p), are given 
in Majumdar [5]. 

Lemma 2.1. If p > 5 is a prime, then 



2 





if 3|(p - 1), 
if 3|(p + 1). 



Lemma 2.2. If p > 5 is a prime, then 



P~ 1, ^ 8|(p- 1), 



2(4,)= P ’ if8 K»>+l). 



3p-l, if 8|(3p- 1), 

3 p, if 8|(3p+l). 
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Lemma 2.3. If p > 7 is a prime, then 



Z(5p) 



p- 1, 


if 10|(p-l), 


p , 


if 10|(p+l), 


2p — 1 , 


if 5 (2p — 1), 


2 p, 


if 5 (2p + 1). 



The explicit expressions of Z{Y2p) and Z(13p) are given in the following two lemmas. 



Lemma 2.4. If p > 13 is a prime, then 



P- 1, 


if 24|(p-l) 


P , 


if 24|(p+l), 


3p- 1, 


if 8 (3p — 1), 


3 P, 


if 8 (3p + 1), 


7p — 1, 


if 24 (7p — 1), 


7 p, 


if 24|(7p+l). 



Proof. By definition, 

• r io ,m(m+ 1) . m(m+l) 

Z(12p) = mm {to : 12p\ } = mm {to : p\ — }. (1) 

If p|to(to + 1), then p must divide either to or to + 1, but not both, and then 24 must divide 
either to + 1 or to respectively. In the particular case when 24 divides p — 1 or p + 1, the 
minimum m in (1) may be taken as p— 1 or p respectively. We now consider the following eight 
cases that may arise : 

Case 1 : p is of the form p = 24 a + 1 for some integer a > 1. 

In this case, 24|(p — 1). Therefore, Z(l2p) = p — 1. 

Case 2 : p is of the form p = 24 a + 23 for some integer a > 1. 

Here, 24|(p + 1), and so, Z(l2p) = p. 

Case 3 : p is of the form p = 24a + 5 for some integer a > 1, so that 8|(3p + 1). 

In this case, the minimum m in (1) may be taken as 3 p. That is, Z(12p) = 3 p. 

Case 4 : p is of the form p = 24a + 19 for some integer a > 1. 

Here, 8|(3p — 1), and hence, Z(12p) = 3p — 1. 

Case 5 : p is of the form p = 24 a + 7 for some integer a > 1. 

In this case, 24|(7p — 1), and hence, Z(12p) = 7 p. 

Case 6 : p is of the form p = 24 a + 17 for some integer a > 1. 

Here, 24|(7p+ 1), and hence, Z(l2p ) = 7p + 1. 

Case 7 : p is of the form p = 24 a + 11 for some integer a > 1. 

In this case, 8|(3p — 1), and hence, Z(12p) = 3p - 1. 

Case 8 : p is of the form p = 24 a + 13 for some integer a > 1. 

Here, 8|(3p+ 1), and hence, Z{12p) = 3 p. 
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Lemma 2.5. For any prime p > 17, 



Z{l3p) = < 



P~ 1, 
P, 

2p-l, 
2 P, 

3p- 1, 
3p, 

4p — 1, 
4p, 

5p — 1, 
5p, 

6p- 1, 
6 P, 



Proof. By definition, 

Z(13p) = : 13p| } 



if 13|0-1), 
if 13|0 + 1), 

if 13|(2p — 1), 
if 13|(2p + 1), 
if 13|(3p — 1), 
if 13|(3p+l), 
if 13|(4p — 1), 
if 13|(4p+l), 
if 13|(5p — 1), 
if 13|(5p+l), 
if 13|(6p — 1), 
if 13|(6p+l). 



min{m : p\ 



m{m + 1) 
26 



}• 



We have to consider the twelve possible cases that may arise : 
Case 1 : p is of the form p = 13a + 1 for some integer a > 1. 
In this case, 13|0 — 1), and so, Z(13p) = p — 1. 

Case 2 : p is of the form p = 13a + 12 for some integer a > 1. 
Here, 13|0+ 1), and hence, Z{13p) = p. 

Case 3 : p is of the form p = 13a + 2 for some integer a > 1. 
In this case, 13|(6p + 1), and hence, Z{13p) = 6 p. 

Case 4 : p is of the form p = 13a + 11 for some integer a > 1. 
Here, 13|(6p — 1), and hence, Z(l3p) = 6p — 1. 

Case 5 : p is of the form p = 13 a + 3 for some integer a > 1. 
In this case, 13|(4p + 1), and hence, Z(l3p) = 4 p. 

Case 6 : p is of the form p = 13a + 10 for some integer a > 1. 
Here, 13|(4p — 1), and hence, Z(13p) = 4p — 1. 

Case 7 : p is of the form p = 13a + 4 for some integer a > 1. 
In this case, 13|(3p + 1), and hence, Z(13p) = 3 p. 

Case 8 : p is of the form p = 13a + 9 for some integer a > 1. 
Here, 13|(3p — 1), and hence, Z(13p) = 3p — 1. 

Case 9 : p is of the form p = 13 a + 5 for some integer a > 1. 
In this case, 13|(5p + 1), and hence, Z(13p) = 5 p. 

Case 10 : p is of the form p = 13a + 8 for some integer a > 1. 
Here, 13|(5p — 1), and hence, Z(13p) = 5p — 1. 

Case 11 : p is of the form p = 13a + 6 for some integer a > 1. 
In this case, 13|(2p + 1), and hence, Z(13p) = 2 p. 



(2) 
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Case 12 : p is of the form p = 13a + 7 for some integer a > 1. 
Here, 13|(2p — 1), and hence, Z(13p) = 2p — 1. 



§3. Main result 



We are now state and prove the main result of this paper in the following theorem. 

Theorem 3.1. There are an infinite number of pairs of dissimilar Pythagorean triangles 
that are pseudo Smarandache related. 

Proof. We consider the pair of dissimilar Pythagorean triangles 

T{3p 1 4p, 5p) and T(5p, 12p, 13p), (3) 

where p is a prime of the form 

p = (2 3 • 3 • 5 • 13)n + 1 = 1560n + 1, n G Z+ . (4) 

By Lemma 2.1 - Lemma 2.5, 

Z(3p) = Z(4p) = Z(5p) = Z(l2p) = Z(13p) = p - 1, 

so that the triangles T{3p, 4p, bp) and T(5p, 12p, 13p) are pseudo Smarandache related. Now, 
since there are an infinite number of primes of the form (4) (by Dirichlet’s Theorem, see, for 
example, Hardy and Wright [3], Theorem 15, pp. 13), we get the desired infinite number of 
pairs of dissimilar Pythagorean triangles that are pseudo Smarandache related. 

It may be mentioned here that, the family of pairs of triangles (3), where p is a prime of 
the form 



p = 1560n — 1, n G Z + , (5) 

also forms (dissimilar) pseudo Smarandache related Pythagorean triangles. 



§4. Some remarks 



The pseudo Smarandache function Z(ri) is clearly not bijective. However, we can define 
the inverse Z~ 1 (m) as follows : 



Z 1 (m) = {n £ Z + : Z(n ) = to} for any integer m > 3, 



( 6 ) 



with 



Z~\ 1) = 1 ,Z“ 1 



(2) = 3. 



(7) 
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As has been pointed out by Majumdar [5], for any m £ Z + , the set Z 1 (jn) is non-empty 

r ffl( Tfl H - 1 ) 

and bounded with as its largest element. Clearly, n £ if and only if the 

following two conditions are satisfied : 

/ 1 \ ., m(m+ 1) 

(1) n divides , 

£(£ 1) 

(2) n does not divide — -for any t with 1 < l < to — 1. 

We can look at (6) from a different point of view : On the set Z + , we define the relation 
SR as follows : 



For any 711,712 £ Z + R 712 if and only if Z{n\) = Z(ri 2 )- (8) 

It is then straightforward to verify that R is an equivalence relation on Z + . It is well- 
known that an equivalence relation induces a partition (on the set Z + ) (see, for example, Gioia 
[2], Theorem 11.2, pp. 32). The sets Z~ 1 (m),'in £ Z + , are, in fact, the equivalence classes 
induced by the equivalence relation R on Z + , and possesses the following two properties : 

OO 

(1) E Z~\m) = Z+ 

m—1 

(2 ) Z -1 (toi) n Z~ 1 (m 2 ) = 0, if mi ^ m 2 . 

Thus, for any n £ Z + , there is one and only one m £ Z + such that n £ Z~ 1 (m). 

A different way of relating two triangles has been considered by Aslibacher [1], which is 
given in the following definition. 

Definition 4.1. Given two triangles, T(a, b 7 c) and T(a , b , c ), where 

a + b + c = 180 = a + b + c , (9) 



they are said to be pseudo Smarandache related if 

Z(a) = Z(a), Z(b) = Z(b'),Z(c) = Z(c). 



The difference in Definition 1.4 and Definition 4.1 is that, in the former, the sides of the 
pair of triangles are pseudo Smarandache related, while their angles, measured in degrees, are 
pseudo Smarandache related in the latter case. Note that the condition (9) merely states the 
fact the sum of the three angles of a triangle is 180 degrees. 

Using a computer program, Ashbacher searched for pseudo Smarandache related pairs of 
triangles (in the sense of Definition 4.1). He reports three such pairs. 

However, in this case, the equivalence classes Z -1 (?n),TO £ Z + , might be of some help. 
The condition (9) can be dealt with by considering the restricted sets Z~ 1 (m\'K) : 

Z -1 (m|7r) = {n £ Z+ : Z(n) = m, 1 < n < 178}. (10) 

Table 4.1 gives such restricted sets related to our problem. 

Clearly, two pseudo Smarandache related triangles T(a, b, c) and T(a ,b , c ) must satisfy 
the following condition : 

a, a £ Z~ 1 (mi\n);b 7 b £ Z -1 (to2|7t), c, c S ^^ 1 (to 3 |7r) for some 7711,7712,777.3 £ Z + . 
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Thus, for example, choosing 



a, a' G {2,6} = Z~\Z\ir)\ b, b G {44,48,88,132,176} = Z~\32\tt), 



we can construct the pseudo Smarandache related triangles T(2,48,130) and T(6,44, 130). 
Again, choosing a, a', 6, b' G {25, 50, 75, 100, 150} = Z _1 (24|7t), we can form the pair of pseudo 
Smarandache related triangles T(25, 150,5) and T(75,100, 5) with the characteristic that 



Z(a) = Z(a') = 24= Z(b) = Z(b') (and Z(c) = Z(c') = 4). 



Choosing a, a ,b,b G {8, 20, 24, 30, 40, 60, 120} = Z 1 (15| 7r) , we see that the equilateral triangle 
T(60, 60, 60) is pseudo Smarandache related to the triangle T(20, 40, 120)! 



Ashbacher [1] cites the triangles T(4, 16, 160) and T(14, 62, 104) as an example of a pseudo 
Smarandache related pair where all the six angles are different. We have found three more, 
given below : 

(1) T( 8, 12, 160) and T(40, 36, 104), 

(with Z( 8) = Z(40) = 15,Z(12) = Z( 36) = 8,Z(160) = Z( 104) = 64), 

(2) T(16, 20, 144) and T(124, 24, 32), 

(with Z( 16) = Z( 124) = 31, Z( 20) = Z( 24) = 15, Z( 144) = Z(32) = 63), 

(3) T(37, 50, 93) and T(74,75,31), 

(with Z(37) = Z(74) = 36, Z(50) = Z(75) = 24, Z(93) = Z(31) = 30). 



Table 4.1. Values of Z 1 (m|7r) = {nG Z + : Z(n) = to, 1 < n < 180} 
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m 


Z 1 (TO7r) 


m 


Z 1 {vn 1 7r ) 


m 


Z 1 (m 7r) 


i 


{1} 


35 


{90} 


80 


{81,108,162,180} 


2 


{3} 


36 


{37,74,111} 


82 


{83} 


3 


{2,6} 


39 


{52,130,156} 


83 


{166} 


4 


{5,10} 


40 


{41,82,164} 


84 


{ 170} 


5 


{15} 


41 


{123} 


87 


{116,174} 


6 


{7,21} 


42 


{43,129} 


88 


{89, 178} 


7 


{4, 14,28} 


43 


{86} 


95 


{152} 


8 


{9,12,18,36} 


44 


{99,110,165} 


96 


{97} 


9 


{45} 


45 


{115} 


100 


{101} 


10 


{11,55} 


46 


{47} 


102 


{103} 


11 


{22,33,66} 


47 


{94,141} 


106 


{107} 


12 


{13,26,39,78} 


48 


{49,56,84,98,147,168} 


108 


{109} 


13 


{91} 


49 


{175} 


111 


{148} 


14 


{35, 105} 


51 


{102} 


112 


{113} 


15 


{8,20,24,30,40,60,120} 


52 


{53,106} 


120 


{121} 


16 


{17,34,68,136} 


53 


{159} 


124 


{125} 


17 


{51,153} 


54 


{135} 


126 


{127} 


18 


{19,57,171} 


55 


{140,154} 


127 


{64} 


19 


{38,95} 


56 


{76,114,133} 


128 


{ 172} 


20 


{42, 70} 


58 


{59} 


130 


{131} 


21 


{77} 


59 


{118,177} 


136 


{137} 


22 


{23} 


60 


{61,122} 


138 


{139} 


23 


{46,69,92,138} 


63 


{32,72,96,112,144} 


148 


{149} 


24 


{25,50,75,100,150} 


64 


{80,104,160} 


150 


{151} 


25 


{65} 


65 


{143} 


156 


{157} 


26 


{27,117} 


66 


{67} 


162 


{163} 


27 


{54,63,126} 


67 


{134} 


166 


{167} 


28 


{29,58} 


69 


{161} 


168 


{169} 


29 


{87, 145} 


70 


{71} 


172 


{173} 


30 


{31,93,155} 


71 


{142} 


178 


{179} 


31 


{16,62,124} 


72 


{73,146} 


255 


{128} 


32 


{44,48,88,132,176} 


78 


{79} 






34 


{85,119} 


79 


{158} 







Ashbacher [1] reports that, an exhaustive computer search for pairs of dissimilar search for 
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all pseudo Smarandache related triangles T(a , b, c) and T(a ,b ,c) ( a+b+c = 180 = a +b +c ) 
with values of a in the range 1 < a < 178, revealed that a cannot take the following values : 



1, 15, 23, 35, 41, 45, 51, 59, 65, 67, 71, 73, 77, 79, 82, 83, 86, 87, 89, 
90, 91, 97, 101, 102, 105, 107, 109, 113, 115, 116, 118, 121, 123, 125, 126, 127, 
131, 134, 135, 137, 139, 141, 142, 143, 148, 149, 151, 152, 153, 157, 158, 159, 
161, 163, 164, 166, 167, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178. 



( 11 ) 



Can the table of the sets Z _1 (m 1 7r) , Table 4.1, be utilized in explaining this observation? Of 
course, if a is too large, then there is a possibility that no two dissimilar triangles exist: Very 
large value of a very often forces the two triangles to be similar. For example, if a = 174 £ 
Z -1 (87|7t), then the three possible pairs of values of ( b , c) are (1, 5), (2, 4) and (3, 3). Note that 
the two sets Z~ 1 ( l|7r) = {1} and Z _1 (2|7t) = {3} are singleton. Thus, if a = 174, then the two 
triangles T(a, b, c) and T(a ,b , c ) must be similar. On the other hand, if a = 116 £ Z~ 1 (87\n), 
then we cannot find b , c with b + c =64. Now, we consider the case when a = 164. This case 
is summarized in the tabular form below : 



a=164 


a'=164 


a =82 


a =41 


Remark 




(1,15) 


Similar Triangle 


Not Possible 


Not Possible 


Both belong to singleton sets 


(b,c) 


(2,14) 


Similar Triangle 


Not Possible 


Not Possible 


b' + c < 34 




(3,13) 


Similar Triangle 


Not Possible 


Not Possible 


3 belongs to a singleton set 




(4,12) 


Similar Triangle 


Not Possible 


Not Possible 


b + c <64 




(5,11) 


Similar Triangle 


Not Possible 


Not Possible 


b + c <65 




(6,10) 


Similar Triangle 


Not Possible 


Not Possible 


b + c <16 




(7,9) 


Similar Triangle 


Not Possible 


Not Possible 


b' + c < 57 




(8,8) 


Similar Triangle 


Not Possible 


Not Possible 







However, if a = 165, we can get two dissimilar pseudo Smarandache related triangles, namely, 
T(165, 7, 8) and T(81, 21, 60). 

A closer look at the values listed in (11) and those in Table 4.1 reveals the following facts: 

(1) The numbers not appearing in any 180 degrees triplets are all belong to singleton sets, with 
the exception of 3 (though we conjecture that 3 cannot appear in any triplet), 47, 64, and 103. 
For the last three cases, we can form the following examples : 

(a) T(47, 76, 57) and T(47, 114, 19) (with Z{ 76) = 56 = Z(114),Z(57) = 18 = Z( 19)), 

(b) T(64, 12, 104) and T(64, 36, 80) (with Z(12) = 8 = Z(36), Z(104) = 64 = Z(80)), 

(c) T(103, 11, 66) and T(103, 55, 22) (with Z( 11) = 10 = Z(55),Z(66) = 11 = Z( 22)). 

(2) In most of the cases, if a value does not appear in the triplet, all other values of the corre- 
sponding Z -1 (m 1 7r) also do not appear in other triplet, with the exception of 145 £ Z~ 1 (29\tt) 
(87 £ Z~ 1 (29\n) does not appear in any triplet), and 146 £ Z _1 (72|7r)(73 £ Z~ 1 (72\n) does 
not appear in any triplet). In this connection, we may mention the following pairs of triplets : 
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(a) T(145, 14,21) and T(145,28,7) (with Z( 14) = 7 = Z(28), Z(21) = 6 = Z( 7)), 

(b) T(146, 4, 30) and T(146, 14,20) (with Z{ 4) = 7 = Z(14),Z(30) = 15 = Z(20)). 

(c) Both 54,63 £ Z _1 (72|7r) can appear in one or the other triplet, but 126 € Z _1 (72|7r) 
cannot appear in any triplet. 

It is also an interesting problem to look for 60 degrees and 120 degrees pseudo Smarandache 
related pairs of triangles, in the sense of Definition 4.1, where a 60 (120) degrees triangle is one 
whose one angle is 60 (120) degrees. We got the following two pairs of pseudo Smarandache 
related 60 degrees triangles : 

(1) T(8, 60, 112) and T(24,60,96), (with Z( 8) = Z(24) = 15,Z(112) = Z(96) = 63), 

(2) T(32, 60, 88) and T(72,60,48), (with Z(32) = Z(72) = 63,Z(88) = Z(48) = 32), 
while our search for a pair of pseudo Smarandache related 120 degrees triangles went in vain. 

We conclude the paper with the following open problems and conjectures. The first two 
have already been posed by Ashbacher [1], 

Problem 1. Are there an infinite number of pairs of dissimilar 60 degrees triangles that 
are pseudo Smarandache related, in the sense of Definition 1.4? 

Problem 2. Is there an infinite family of pairs of dissimilar 120 degrees triangles that are 
pseudo Smarandache related, in the sense of Definition 1.4? 

Ashbacher [1] reports that a limited search on a computer showed only four pairs of dis- 
similar 60 degrees psedo Smarandache related triangles, while the number is only one in the 
case of 120 degrees triangles. 

In the formulation of Definition 4.1, the number of pairs of dissimilar Smarandache related 
or pseudo Smarandache related triangles is obviously finite. In fact, since the number of integer 
solutions to 

| ct T b T c = 180 
1 a > 1,6 ^ l,c> 1 

is C( 179,2) (see, for example, Johnsonbaugh [4], Theorem 4.5, pp. 238), the number of such 
triangles cannot exceed C(179, 2). But this is a very crude estimate. The next problems of 
interest are 

Problem 3. Is it possible to find a tight or better upper limit to the number of pairs of 
dissimilar triangles that are pseudo Smarandache related, in the sense of Definition 4.1? 

Problem 4. Is it possible to find a good upper limit to the number of pairs of dissimilar 
60 degrees triangles that are pseudo Smarandache related, in the sense of Definition 4.1? 

Problem 5. Is it possible to find another pair of 60 degrees pseudo Smarandache related 
triangles, in the sense of Definition 4.1, where all the six angles are acute? 

By a limited search, we found only two pairs of dissimilar 60 degrees pseudo Smarandache 
related triangles, the ones already mentioned above. 

Conjecture 1. There is no pair of dissimilar 120 degrees triangles that are Pseudo Smaran- 
dache related, in the sense of Definition 4.1. 

Conjecture 2. There is no pair of dissimilar 90 degrees triangles that are Pseudo Smaran- 
dache related, in the sense of Definition 4.1. 
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Abstract For any positive integer n, the famous F. Smarandache function S(n) is defined as 
the smallest positive integer m such that n \ m\. That is, 5(n) = min{m : n \ ml, n £ N}. 
The main purpose of this paper is using the elementary method to study the number of the 

solutions of the congruent equation l s ( n_1 ) + 2 s C-i) -| -f (n — + 1 = 0 ( mod n), 

and give its all prime number solutions. 

Keywords F. Smarandache function, divisibility, primitive root. 



§1. Introduction and result 



For any positive integer n, the famous F. Smarandache function S(n) is defined as the 
smallest positive integer m such that n | ml. That is, 5(n) = min {to : n | m!, n £ N}. For 
example, the first few values of S(n) are 5(1) = 1, 5(2) = 2, 5(3) = 3, 5(4) = 4, 5(5) = 5, 
5(6) = 3, 5(7) = 7, 5(8) = 4, 5(9) = 6, 5(10) = 5, 5(11) = 11, 5(12) = 4, • • • . About the 
elementary properties of 5(n), many authors had studied it, and obtained a series interesting 
results, see references [1], [2], [3] and [4]. For example, Xu Zhefeng [2] studied the value 
distribution problem of 5(n), and proved the following conclusion: 

Let P(n) denotes the largest prime factor of n, then for any real number x > 1, we have 
the asymptotic formula 



E - p w ) 2 

n<.x 



2C(iH 

3 In a; 



O 




where £(s) is the Riemann zeta-function. 

Lu Yarning [3] studied the solutions of an equation involving the F. Smarandache function 
5(n), and proved that for any positive integer k > 2, the equation 



5(toi + m 2 + • • • + m k ) = 5(toi) + 5(to 2 ) + • • • + S(m k ) 



has infinite group positive integer solutions (mi, m 2 , • • • , rri k ). 

Rongji Chen [4] proved that for any fixed r £ N with r > 3, the positive integer n is a 
solution of the equation 



S(n) r + 5(n) r_1 + • • • + 5(n) = n 
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if and only if 

n = p{p r - 1 +p r ~ 2 + ■•• + !), 
where p is an odd prime such that 

p r ~ 1 +p r - 2 + --- + 1 | (p-1)!. 

On the other hand, in reference [5], C. Dumitrescu and V. Seleacu asked us to study the 
solvability of the congruent equation 



1 S(n-l) + 2 S(n-l) + l)S(n-l) + 1 = Q mod 



n. 



(i) 



About this problem, it seems that none had studied it yet, at least we have not seen any 
related papers before. In this paper, we use the elementary method and the properties of the 
primitive roots to study the solvability of the congruent equation ( 1 ), and obtain its all prime 
solutions. That is, we shall prove the following conclusion: 

Theorem. Let n be a prime, then n satisfy the congruent equation 

1 S(n-i) + 2 S(n- 1 ) + i)S(n-i) + 1 = 0 mod n 



if and only if n = 2, 3 and 5. 

It is clear that our Theorem obtained all prime solutions of the congruent equation (1). 
About the general solutions of the congruent equation ( 1 ) is still an unsolved problem. 

2. Proof of the theorem 

In this section, we shall complete the proof of our Theorem directly. We only discuss the 
prime solutions of ( 1 ). 

(1) For n = 2, since 2 | l 5 ! 1 ) + 1 = 2 , so n = 2 is a prime solution of ( 1 ). 

(2) For n = 3, since 3 | l s ( 2 ) + 2 S ( 2 ) + 1 = 6 , so n = 3 is a prime solution of (1). 

(3) For n = 5, since 5 | l s ( 4 ) + 2 S ( 4 ) + 3 S ( 4 ) + 4 5 ^ 4 ) + 1 = 355, so n = 5 is also a prime 
solution of ( 1 ). 

(4) For prime n = p > 7, it is clear that p has at least a primitive root. Let g be a primitive 

root of p, that is to say, (g l — 1 , p) =1 for all 1 < i < p — 2 , the congruences 

p p_1 = 1 mod p and g rn ^ p ~ 1 '> = 1 mod p ( 2 ) 

hold for any positive integer m. 

Then from the properties of the primitive root mod p we know that g°, g 1 , • • • , g p ~ 2 is a 
reduced residue class. So we have the congruent equation 

1 S(n-l) + 2 S(n-l) + !)S(n-l) ( 3 ) 

= gO-S(p-l) + gl-S(p-l) ^ h (p-2).s(p-l) 

g(p-l)-S(p~l) _ l 

~ mod p. 
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It is clear that for any prime p > 7, we have S(p— 1) < p— 2, so (g s ( p ^ —l,p)= 1. Therefore, 
from (2) and (3) we have 



(p— 1)-S(p— 1) _ -i 

1 S'(n— 1) + 2 S(n- 1) + l)S(n-l) = 9 ^ = Q mod p. 

v ’ gS(p-l) _ 1 

So from this congruence we may immediately get 

1 S(n-l) + 2 S(n- 1) + . . . + ( n _ l)S(n-l) + 1 = 1 mod p _ 



Thus, if prime p > 7, then it is not a solution of (1). This completes the proof of Theorem. 
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Abstract We introduce the notion of <r-finite Ponomarev-system (f,M,X,{V n }) to give a 
consistent method to construct an mssc-mapping / with covering-properties onto a space 
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Keywords cr-finite Ponomarev-system, relatively compact metric, sn-network, so-network, 

cs-network, ^-network, 2-sequence-covering, 1-sequence-covering, sequence-covering, compact- 
covering, mssc-mapping. 



§1. Introduction 

An investigation of relations between spaces with countable networks and images of sep- 
arable metric spaces is one of interesting questions on generalized metric spaces. Related to 
characterizing a space X having a certain countable network V by an image of a separable 
metric space M under some covering-mapping /, many results have been obtained. In the past, 
E. Michael [18] proved the following results. 

Theorem 1.1. ([18], Proposition 10.2) The following are equivalent for a space X. 

(1) A is a cosmic space. 

(2) X is an image of a separable metric space. 

Theorem 1.2. ([18], Theorem 11.4) The following are equivalent for a space X. 

(1) X is an H 0 -space. 

(2) A is a compact-covering image of a separable metric space. 

Also, he posed the following question, on page 999. 

Question 1.3. Find a better explanation for images of separable metric spaces? 

In [26], Y. Tanaka and Z. Li proved the following result. 

Theorem 1.4. ([26], Corollary 8 (1)) The following are equivalent for a space A, where 
“subsequence-covering” can be replaced by “pseudo-sequence-covering”. 

(1) A has a countable cs*-network (resp., cs-network, sn-network). 
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(2) A is a subsequence-covering (resp., sequence-covering, 1-sequence-covering) image of a sep- 
arable metric space. 

Recently, Y. Ge [8] sharpened Theorem 1.2 as follows. 

Theorem 1.5. ([8], Theorem 12) The following are equivalent for a space X. 

(1) X is an H 0 -space. 

(2) X is a sequence-covering, compact-covering image of a separable metric space. 

(3) X is a sequentially-quotient image of a separable metric space. 

Taking these results into account, the following question naturally arises. 

Question 1.6. When V is a countable network (cs-network, sn-network, so-network) for 
X; how nice can the mapping / and the metric domain M be taken to? 

Around this question, the author of [2] has shown that, when V is a countable cs-network, 
/ and M in Theorem 1.5 can be an mssc-mapping and a relatively compact metric space, 
respectively. More precisely, the following has been proved. 

Theorem 1.7. ([2], Theorem 2.5] The following are equivalent for a space X. 

(1) X is an H 0 -space. 

(2) A is a sequence-covering, compact-covering mssc- image of a relatively compact metric 
space. 

(3) X is a sequentially-quotient image of a separable metric space. 

In the above result, / and M can not be any compact mapping and any compact metric 
space, respectively; see [2, Example 2.11 & Example 2.12]; and the key of the proof is to 
construct a sequence-covering, compact-covering mssc-mapping / from a relatively compact 
metric space M onto X. 

By the above, we are interested in the following question. 

Question 1.8. Find a consistent method to construct an mssc-mapping with covering- 
properties onto a space from some relatively compact metric space? 

In this paper, we answer Question 1.8 by introducing the notion of cr-finite Ponomarev- 
system (/, M, X, {P n }) to give a consistent method to construct an mssc-mapping / with 
covering-properties onto a space X from some relatively compact metric space M . As applica- 
tions, we systematically obtain characterizations on mssc-images of relatively compact metric 
spaces under certain covering-mappings, which gives an answer for Question 1.6, particularly, 
for Question 1.3. 

Throughout this paper, all spaces are regular and Ti, N denotes the set of all natural 
numbers, ui = N U {0}, and a convergent sequence includes its limit point. Let V be a family 
of subsets of X. Then ]JP , and f]V denote the union ]J{P : P 6 P}, and the intersection 
f]{P : P & V}, respectively. A sequence {x n : n £ u>} converging to x 0 is eventually in Ac X, 
if {x n : n > no}U{xo} C A for some no £ N, and it is frequently in A if {x nk : k £ N}U{a;o} C A 
for some subsequence {x nk : k £ N} of {x n : n £ N}. 

Definition 1.9. Let V be a family of subsets of a space X. 

(1) For each x £ A', V is a network at x in A, if a; £ f) V, and if x £ U with U open in A, then 
there exists P £ V such that x £ P C U. 

(2) V is a cs-network for A [10] if, for every convergent sequence S converging to x £ U with 
U open in A, there exists P £ V such that S is eventually in P C U. 
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(3) V is a cs*-network for X [5] if, for every convergent sequence S converging to x £ U with 
U open in X , there exists P £ V such that S is frequently in P C U. 

(4) V is a fc-network of X [19] if, for every compact subset K C U with U open in A', there 
exists a finite T C V such that K c [J T C U. 

Definition 1.10. Let V = \_}{P X '■ x £ X} be a family of subsets of a space X satisfying 
that, for each x £ X, V x is a network at x in A, and if U, V £ V x , then W C U D V for some 
W£V X . 

(1) V is a weak base for X [21], if G C A' such that for each x £ G, there exists P £ ’Pa- 
satisfying P C G, then G is open in X. 

(2) P is an sn-network for X [14], if each member of V x is a sequential neighborhood of x in 
A. 

(3) P is an so-network for X [12], if each member of P x is sequentially open in X. 

(4) The above V x is respectively a weak base, an sn-network , and an so-network at x in X [12]. 

Definition 1.11. ([4]) Let P be a subset of a space X. 

(1) P is a sequential neighborhood of x if, for every convergent sequence S converging to x in 
A, S is eventually in P. 

(2) P is a sequentially open subset of X if, for every x £ P, P is a sequential neighborhood of 

x. 

Definition 1.12. Let X be a space. 

(1) X is relatively compact, if A is compact. 

(2) X is an H 0 -space [18] (resp., sn-second countable space [7], so-second countable space, g- 
second countable space [21], second countable space [3]), if X has a countable cs-network (resp., 
countable sn-network, countable so-network, countable weak base, countable base). 

(3) A is a sequential space [4], if every sequentially open subset of A is open. 

Remark 1.13. ([15]) (1) For a space, weak base =>■ sn-network => cs-network. 

(2) An sn-network for a sequential space is a weak base. 

Remark 1.14. (1) It is easy to see that “compact metric =>■ relatively compact metric => 
separable metric”, and these implications can not be reversed from Example 2.14 and Example 
2.15. 

(2) It is well-known that a space X is an H 0 -space if and only if X has a countable fc-network 
(cs*-network), see [23, Proposition C], for example. 

Remark 1.15. Let / : X — > Y be a mapping. 

(1) / is an mssc-mapping[13], if X is a subspace of the product space X n of a family 

nEN 

{X n : n £ N} of metric spaces, and for every y £ Y, there exists a sequence { V v , n : n £ N} of 
open neighborhoods of y in Y such that each p n (f~ 1 {V Vtn )) is a compact subset of X n , where 
Pn ■ n Xi — > X n is the projection. 

i£N 

(2) / is an 1-sequence-covering mapping [14] if, for every y £ Y, there exists x v £ f^ 1 (y) 
such that whenever {y n '■ n £ N} is a sequence converging to y in Y there exists a sequence 
{x n : n £ N} converging to x y in X with each x n £ f~ 1 (y n ). 

(3) / is a 2-sequence-covering mapping [14] if, for every y £ Y, x £ / -1 (y), and sequence 
{y n :neN] converging to y in Y, there exists a sequence {x n : n £ N} converging to x in X 
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with each x n £ f 1 (y n )- 

(4) / is an 1-sequentially quotient mapping [17] if, for every y £ Y, there exists x y £ f~ l {y) 
such that whenever {y n : n £ N} is a sequence converging to y in Y there exists a sequence 
{xk '■ k £ N} converging to x y in X with each Xk £ f~ 1 (yn k )- 

(5) / is a sequence-covering mapping [20] if, for every convergent sequence S in Y, there exists 
a convergent sequence L in X such that f(L) = S. 

(6) / is a compact-covering mapping [18] if, for every compact subset K of Y, there exists a 
compact subset L of X such that f(L) = K. 

(7) / is a pseudo-sequence-covering mapping [11], if for every convergent sequence S of Y, there 
exists a compact subset K of X such that f(K) = S. 

(8) / is a sequentially-quotient mapping [1] if, for every convergent sequence S in Y , there 
exists a convergent sequence L in X such that f(L) is a subsequence of S. 

For terms are not defined here, please refer to [3, 24]. 



§2. Results 

Definition 2.1. Let V = {P n : n £ N} be a countable network for a space X. Because 
X is T\ and regular, we may assume that every member of V is closed. For each n £ N, let 
V n = {X} U {Pi : i < n} = {P a : a £ A n }, where A n is a finite set, and let every A n be 
endowed with the discrete topology. Setting 

M = ja = (a*) £ : {P ai : i £ N} forms a network at some point x a £ xj. 

ieN 

Then M, which is a subspace of the product space Ai, is a metric space. Since X is T\ and 

»e N 

regular, x a is unique for each a £ M. We define / : M — > X by /(a) = x a for each a £ M. 
Then / is an ?nssc-mapping from a relatively compact metric space onto X by the following 
Theorem 2.3. The system (f,M,X,{P n }) is a cr-finite Ponomarev-system. 

Remark 2.2. For more details on Ponomarev-systems and images of metric spaces, see 
[9,16,25], for example. 

Theorem 2.3. Let (/, M, X, {V n }) be the system in Definition 2.1. Then the following 
hold. 

(1) M is a relatively compact metric space. 

(2) / is an mssc-mapping. 

Proof. (1) For every i £ N, since A t is finite, Aj is compact metric. Then At is compact 

metric, and so M C Ai is. It implies that M is relatively compact metric. 

ieN 

(2) It suffices to prove the following facts (a), (b), and (c). 

(a) / is onto. 

Let x £ X. Since V is a countable network for X, (' P) x = {P £ V : x £ P} is a countable 
network at x in X. We may assume that {V) x = {Px,j '■ j € N}, where P x j £ Vi( :] ) for some 
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i(j) £ N satisfying i(j ) < i(j + 1). For each i £ N, taking P ai = if i = i(j) for some j £ N, 
and otherwise, P ai = A. Then (a,) £ and {P ai : i £ N} = (V) x forms a network at x 

ie N 

in X. Setting a = (a,), we get a £ M and /(a) = x. 

(b) / is continuous. 

Let x = /(a) £ U with U open in X and a £ M. Setting a = (a*) £ A i: where {P ai : i £ N} 

ieN 

forms a network at x in A'. Then there exists n £ N such that x £ P an C U. Setting 
M a = {b = (fifi £ M : fi n = a n }, then M a is an open neighborhood of a in M. For each 
b £ M a , we get f(b) £ Pg n = P arl C U. It implies that f(M a ) C U. 

(c) / is an mssc-mapping. 

Let x £ X. For each * £ N, taking 1 4,i = X. Then {V Xt i : i £ N} is a sequence of open 
neighborhoods of x in X. Since Ai is finite, Ai is compact. Then Pi(f~ 1 (V X} i)) = Pi(f~ 1 ( A)) C 
Ai is compact. It implies that / is an mssc-mapping. 

Remark 2.4. In the view of (a) in the proof (1) => (2) of [2, Theorem 2.5], / is onto by 
choosing a = (afi satisfying /(a) = x , where x £ P ai £ Vi. Unfortunately, this argument is 
not true. In fact, we may pick x £ P ai = X £ V z . Then {P ai : i £ N} = {A} does not form 
a network at x in X. It is a contradiction. The wrong argument is corrected as in the proof 

(2). (a) of Theorem 2.3. 

Theorem 2.5. Let (/, M, X, {V n }) be a er-finite Ponomarev-system. Then the following 
hold. 

(1) If V is a cs-network for A, then / is sequence-covering. 

(2) If V is a fc-network for A, then / is a compact-covering. 

(3) If V is a cs*-network for A, then / is pseudo-sequence-covering. 

(4) If V is an sn-network for A, then / is 1-sequence-covering. 

(5) If V is an so-network for A, then / is 2-sequence-covering. 

Proof. (1) Let S = {x m : m £ ui} be a convergent sequence converging to Xq in A. 
Suppose that U is an open neighborhood of S' in A. A family A of subsets of A has property 
cs(5, U) if: 

(i) A is finite. 

(ii) for each P£ A, tt^PnScPcU. 

(iii) for each x m £ S, there exists unique P Xm £ A such that x m £ P Xm . 

(iv) if x 0 £ P £ A, then S' — P is finite. 

For each i £ N, since A = {A} C Vi has property cs(S, A) and Vi is finite, we can assume 

that 

{A C Vi : A has property cs(S, A)} = { A : j = n^i + 1, . . . , nj, 

where no = 0. By this notation, for each j £ N, there is unique i(j) £ N such that A^j) has 
property cs(S, A). Then for each j £ N, we can put Ai(j) = {P a : a £ Ej }, where Ej is a finite 
subset of Aj . 

For each j £ N, m £ u> and x m £ S, it follows from (iii) that there is unique o:y m £ Ej 
such that x m £ P ajtm € A^j). Let a m = (<^j,m) £ n <= n Aj. We shall prove that 

ieN jeN 

{Puj m '■ j £ N} is a network at x m in A. In fact, let x m £ U with U open in A, we consider 
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two following cases. 

(a) If m = 0, then S is eventually in P XQ C U for some P Xo £ V. For each x £ S — P Xo , 

let x £ P x C X — (S — {cc}) for some P x £ V. Then Q = {P Xo } U {P x : x £ S — P Xo } C V 
has property cs(S,X). Since Q is finite, Q C Vi for some iel. It implies that Q = A(j) 
for some i(j) £ N with j £ {m - i + 1 , . . . , rZj}. Since x m = x 0 £ P aj 0 , P ajt0 = Px 0 - Hence 
Xm — XO £ Pctj,a P • 

(b) If m ^ 0, then S — {a; m } is eventually in P Xo C X — {x m } for some P Xo £ V. For each 

x £ (S — {x m }) — P Xo , let x £ P x C X — (S — {a;}) for some P x £ V, and let x m £ P Xm C 
[/n(X-(S-{a; m })) for some P Xm £ V. Then U = {H a;o }U{P Xm }U{P x : x £ (5-{a: m })-P l0 } 
has property cs(S,X). Since H is finite, H C Vi for some i £ N. It implies that H = Ayy 
for some i(j ) £ N with j £ {nj_i + 1 Since x m £ Pa j<m , Pa jtm = P Xm ■ Hence 

Xm £ Potj^rn P ' 

By the above, for each m £ u> we get a m = (ay m ) £ M satisfying /(a m ) = x m . For each 
j £ N, since families H and Q are finite, there exists rrij £ N such that ay m = ay o if m > rrij. 
Hence the sequence {ay m ’■ m £ N} converges to ay o in Aj. Thus, the sequence {a m : m £ N} 
converges to ao in M. Setting L = {a m : m £ w}, then L is a convergent sequence in M and 
f(L) = S. This shows that / is sequence-covering. 

(2)Let K be a compact subset of X. Suppose that V is an open neighborhood of K in X. A 
family B of subsets of X has property k(K. V) if: 

(i) B is finite. 

(ii) P fl K ^ 0 for each P £ B. 

(hi) Kc\JBcV. 

For each i £ N, since B = {X} C Vi has property k(K , X) and Vi is finite, we can assume 

that 

{B C Vi : B has property k(K , X)} = { B i(j ) : j = rii-i + 1, . . . , nj, 

where no = 0. By this notation, for each j £ N, there is unique i(j) £ N such that B, tjj has 
property k(K,X). Then for each j £ N, we can put Hyy = {P a : a £ Fj }, where Fj is a finite 
subset of Aj . 

Setting L = {a = (aj) £ ]^[ Fj : P (K fl P a J ^ 0}, we shall prove that L is a compact 

ieN -ieN 

subset of M satisfying /(F) = K, hence / is compact-covering, by the following facts (a), (b), 
and (c). 

(a) L is compact. 

Since L C Fj and Fj is compact, we only need to prove that L is closed in Fj. Let 
iSN ieN ieN 

a = (aj) £ Fj — L. Then P (KnP ai ) = 0. Since KnP ai is compact, there exists zo £ N such 

ieN i£N 

that P| (A'nP a J = 0. Setting W = {b = (/3j) £ Fj : /3j = aj if z < z 0 }. Then W is an open 

neighborhood of a in Fj and W fl L = 0. If not, there exists b = (/3j) £ W fl L. Since b £ F, 

zeN 

P (K fl Pp J 7^ 0, hence P (K fl Pp J ^ 0. Since b £ W, P (K n P ai ) = p (ifnP ft )/0. 

iGN i<io i<i 0 i<i 0 
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This is a contradiction of the fact that P (K ft P ai ) = 0. 

i<i 0 

(b) L C M and f(L) C K. 

Let a = ( at ) £ L , then a £ n F.j and P(/\ (~1 P Qi ) / 0- Taking x £ P (K fl P Qi ). If 

ieN i£N ieN 

{P ai : f € N} is a network at x in X , then a £ M and f(a) = x, hence L C M and f(L) C K. 

So we only need to prove that {P ai : i £ N} is a network at x in X. Let V be an open 
neighborhood of x in X. There exist an open subset W of K such that x £ W, and compact 
subsets W and K — W such that W C V and K — W C X — {a?}, where W is the closure 
of W in K. Since V is a fc-network for X, there exist finite families T\ C V and C V such 
that W K clJPcy and K—W C [JP 2 C X — {x}. We may assume that P fl K ^ 0 for each 
P £ T\ U P 2 - Setting T = U P 2 , then T has property k(I \ , X). It implies that T = B^j) for 
some i(j ) £ N with j £ {n,_i + 1, . . . , rij}. Since x £ P aj £ B i ^' ) ,P aj £ T\. Then x £ P a . C V. 
This prove that {P a . : j £ N} is a network at x in X. 

(c) I< c f(L). 3 

Let x £ K. For each i £ N, there exists £ F t such that x £ P ai . Setting a = (a*), then 
a £ L. Furthermore, /(a) = x as in the proof of (b). So K C f{L). 

(3) Let S' be a convergent sequence of X. By using [26, Lemma 3] and as in the proof of (2), 
here S plays the role of the compact subset K , there exists a compact subset L of M such that 
f(L) = S. Then, / is pseudo-sequence-covering. 

(4) Let V = [j{V x : x £ X}, where each V x is an sn-network at x in X. Let x £ X. We may 
assume that V x = {P x ,j ■ j £ N}, where P x j e ^*(j) with some i(j ) £ N satisfying i(j) < i(j+ 1). 
For each * £ N, take P ai = P-ijj) if i = i(j) for some j £ N, and otherwise, P ai = X. Then 
(c*i) G n Ai, and {P ai : i £ N} = V x forms a network at x in X. Setting a x = (ai), then 

ieN 

a £ M and a x £ / - 1 (ie). For each n £ N, setting B n = {b = (Pi) £ M : Pi = a* if i < n}. 
Then {B n : n £ N} is a decreasing neighborhood base at a x in M. For each n £ N, we get 
f(B n ) = P| P a .. In fact, for each b = (Pi) £ B n , f(b) = p Pp. £ Pp. = P a . if i < n. Then 

i<~n iGN 

m £ p P a . . It implies that f(B n ) C P P Qi . Conversely, if x £ P P ai , then x = f(b) for 

i<n z<n i<n 

some b = (Pi) £ M. Setting c = ("fp, where 7 * = ai if i < n, and 7 j = Pi- n if i > n. Then 
{P 7i : i £ N} is a network at x in X. For each i > n, since Pi- n £ Ai_ n C Aj, 7 j £ A, : . Then 
c £ B n and /(c) = x. It implies that P P ai C f(B n ). 

i<n 

Suppose that the sequence {xj : j £ N} converges to x in X. For each isN, since P ai is 
a sequential neighborhood of x in X , there exists j) £ N such that x :j £ P a . for each j > ji. 
For each n £ N, setting j(n) = max{j,; : i < n}, then Xj £ P P Qi = f(B n ) for each j > j(n). 

i<n 

It implies that f~ 1 (xp nB„/0 for each j > j(n). We may assume that 1 < j(n) < j(n + 1). 
For each j £ N, taking 



a j e 



/ 1 &j) if 

/ _1 (^)nP„ if j(n) < j < j(n+ 1 ). 



Then the sequence {aj : j £N} converges to a x in M. In fact, let U be a neighborhood of 
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a x , then there exists n £ N such that a x £ B n C U. Obviously, aj £ B n C U for each j > j(n), 
so {a,j : j £ N} converges to a x . Note that f(uj) = Xj for every j £ N. It implies that / is 
1-sequence-covering . 

(5) Let x £ X and a x £ f~ 1 (x). Setting a x = (a*), then on £ Aj for each i £ N, and {P ai : i £ 
N} forms a network at x in X. For each n £ N, setting B n = {b = ( fli ) £ M ■. (3i = on \i i < n} . 
Then {B n : n £ N} is a decreasing neighborhood base at a x in M, and for each n £ N, 
f{B n ) = a<„ P ai as in the proof of (4). 

Suppose that sequence {xj : j £ N} converges to x in X. For each i £ N, since P ai is 
sequentially open, there exists ji £ N such that Xj £ P ai for each j > ji. As in the proof (4), 
we get a sequence {aj : j £ N} converges to a x in M such that f{a,j) = Xj for each j £ N. It 
implies that / is 2-sequence-covering. 

By using Theorem 2.3 and Theorem 2.5 above, we obtain characterizations of ?nssc-images 
of relatively compact metric spaces under certain covering-mappings, which sharpen results in 
[18,26] and more. Firstly, we sharpen Theorem 1.1 as following. 

Corollary 2 . 6 . The following are equivalent for a space X. 

(1) A is a cosmic space. 

(2) X is an mssc-image of a relatively compact metric space. 

(3) X is an image of a separable metric space. 

Proof. We only need to prove (1) => (2), other implications are routine. Since X is a cosmic 
space, X has a countable network V . Then the cr-finite Ponomarev-system (/, Af, A', {V n }) 
exists. Therefore, / is an mssc-mapping and M is relatively compact metric by Theorem 2.3. 
It implies that A is an mssc-image of a relatively compact metric space. 

Remark 2 . 7 . Corollary 2.6 is an answer of Question 1.3. 

Next, we extend partly Theorem 1.4 as follows. 

Corollary 2 . 8 . The following are equivalent for a space A. 

(1) A is an s?r-second countable space. 

(2) A is an 1-sequence-covering, compact-covering mssc-image of a relatively compact metric 
space. 

(3) A is an 1-sequentially-quotient image of a separable metric space. 

Proof. (1) => (2). Since A is an sn-second countable space, A has a countable closed 
sn-network Vi and a countable closed fc-network P 2 - Set V = 'P l UV 2 , then V is a countable sn- 
network and fc-network for A. It follows from Definition 2.1 that the cr-finite Ponomarev-system 
(/, M, A, {Vn}) exists. Therefore, / is an 1-sequence-covering, compact-covering mssc-mapping 
from a relatively compact metric space M onto A by Theorem 2.3 and Theorem 2.5. 

(2) => (3). It is obvious. 

(3) (1). Let / : M — > A be an 1-sequentially-quotient mapping from a separable 
metric space M onto A. Since M is separable metric, M has a countable base B. We may 
assume that B is closed under finite intersections. For each x £ A, there exists a x £ f~ 1 {x) 
such that whenever {x n : n £ N} is a sequence converging to x in A there exists a sequence 
{afc : k £N} converging to a x in M with each a*, £ f~ 1 (x nk ). Let B x = {B £ B : a x £ B} and 
V x = f{B x ). We shall prove that V = [_ }{V x ■ x £ A} is a countable sn-network for A. Since 
B is countable, V is countable. It suffices to prove the following facts (a), (b), and (c) for every 
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x£X. 

(a) V x is a network at x in X. 

It is straightforward because B x is a neighborhood base at a x in M. 

(b) If Pi, P 2 £ V x , then P C P\ fl P 2 for some P £V X . 

Let Pi = f(Bi),P 2 = f(B 2 ) with Bi,B 2 £ B x . Then Bi D B 2 £ B x . Setting P = 
f{Bi D B 2 ), we get P £V X and P C Pi fl P 2 . 

(c) Each P £ V x is a sequential neighborhood of x. 

Let P = f(B ) with B £ B x , and let {x n : n £ N} be a sequence converging to x in 
X. For each subsequence {x ni : i £ N} of {x n : n £ N}, since / is 1-sequentially-quotient, 
there exists a sequence {a,k : k £ N} converging to a x in M with each a,k £ f~ 1 {x 1li k ). 
Since {a fc : k £ N} U {a x } is eventually in B , {x ni :ieN}U {x} is frequently in P. Then 
{x n :neN}U {a;} is frequently in P. It follows from [6, Remark 1.4] that P is a sequential 
neighborhood of x. 

Corollary 2 . 9 . The following are equivalent for a space X. 

(1) X is a g-second countable space. 

(2) X is an 1-sequence-covering, compact-covering, quotient mssc-image of a relatively compact 
metric space. 

(3) X is an 1-sequentially-quotient, quotient image of a separable metric space. 

Proof. (1) => (2). It follows from Corollary 2.8 that X is an 1-sequence-covering mssc- 
image of a relatively compact metric space under some mapping /. Since / is an 1-sequence- 
covering mapping onto a g-second countable space, / is a quotient mapping by [15, Lemma 
3.5]. It implies that X is an 1-sequence-covering, quotient mssc-image of a relatively compact 
metric space. 

(2) =>■ (3). It is obvious. 

(3) => (1). Since X is a quotient image of a separable metric space, X is a sequential space 
by [4, Proposition 1.2]. Then A' is a sequential space having a countable sn-network V by 
Corollary 2.8. It follows from Lemma 1.13 that V is a weak base for X. Then X is a g-second 
countable space. 

In the following, we prove a result on spaces having countable so-networks. 

Corollary 2 . 10 . The following are equivalent for a space X. 

(1) X is an so-second countable space. 

(2) A is a 2-sequence-covering, compact-covering mssc-image of a relatively compact metric 
space. 

(3) X is a 2-sequence-covering image of a separable metric space. 

Proof. (1) => (2). Since X is an so-second countable space, X has a countable closed 
so-network V 1 and a countable closed fc-network V 2 . Set V = Vi UV 2 , then V is a countable so- 
network and fc-network for X. It follows from Definition 2.1 that the cr-finite Ponomarev-system 
(/, M, X, {Vn}) exists. Therefore, / is a 2-sequence-covering, compact-covering mssc-mapping 
from a relatively compact metric space M onto X by Theorem 2.3 and Theorem 2.5. 

(2) => (3). It is obvious. 

(3) => (1). Let / : M — > X be a 2-sequence-covering mapping from a separable metric 
space M onto X. Since M is separable metric, M has a countable base B. For each x £ X, 
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let B x = {B £ B : f~ 1 (x) fl B / 0}, and V x be the family of finite intersections of members 
of f(B x ). We shall prove that V = IJ {V x : x £ A'} is a countable so-network for X. Since B 
is countable, V is countable. It suffices to prove the following facts (a), (b), and (c) for every 
x £ X. 

(a) V x is a network at x in A'. 

It is straightforward because B x is a base of f _1 (x) in M. 

(b) If Pi, P 2 £ V x , then P C Pi (~1 P 2 for some P GV X . 

It is straightforward by choosing P = Pi n P 2 . 

(c) Each P £ V x is sequentially open. 

For each y £ f(B) with B £ B x , we get f~ 1 {y ) D B/0. Let {y n : n £ N} be a sequence 
converging to y in A. Then there exists a sequence {a n : n £ N} converging to a £ f~ 1 {y) D B 
such that f(a n ) = y n for each n £ N. Since B is open, {a n : n £ N} U {a} is eventually in B. It 
implies that {y n : n £ N}U{y} is eventually in /(P). Hence f(B) is a sequential neighborhood 
of y in X. Therefore, f(B) is sequentially open. Since each P £ V x is some intersection of 
finitely many members of f(B x ), P is sequentially open. 

Corollary 2.11. The following are equivalent for a space X. 

(1) A is a second countable space. 

(2) A is a 2-sequence-covering, compact-covering, quotient mssc - image of a relatively compact 
metric space. 

(3) A is a 2-sequence-covering, quotient image of a separable metric space. 

Proof. (1) =>• (2). It follows from Corollary 2.10 that A is a 2-sequence-covering mssc- 
image of a relatively compact metric space under some mapping /. Since / is a 2-sequence- 
covering mapping onto a second countable space, / is a quotient mapping by [15, Lemma 3.5]. 
It implies that A is a 2-sequence-covering, quotient mssc - image of a relatively compact metric 
space. 

(2) => (3). It is obvious. 

(3) =>• (1). Since X is a quotient image of a separable metric space, A is a sequential 
space by [4, Proposition 1.2]. Then X is a sequential space having a countable so-network V 
by Corollary 2.10. For each P £ V, since P is sequentially open and A is sequential, P is open. 
Then V is a countable base for X, i.e. , A is a second countable space. 

Remark 2.12. By using Theorem 2.3 and Theorem 2.5, we also get the proof of Theorem 
1.7 again, which was presented in [2]. 

From the above results, we obtain preservations of certain spaces under covering-mappings 
as follows. 

Corollary 2.13. Let / : X — > Y be a mapping. Then the following hold. 

(1) If X is an sn-second countable space and / is 1-sequentially-quotient, then Y is an sn-second 
countable space. 

(2) If A is a g-second countable space and / is 1-sequentially-quotient, quotient, then Y is a 
^-second countable space. 

(3) If X is an so-second countable space and / is 2-sequence-covering, then Y is an so-second 
countable space. 

(4) If X is a second countable space and / is 2-sequence-covering, quotient, then Y is a second 
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countable space. 

Proof. We only need to prove (1), by similar arguments, we get the other. Since X is an 
sn-second countable space, A' is an 1-sequentially-quotient image of a separable metric space 
under some mapping g by Corollary 2.8. It is easy to see that / o g is also an 1-sequentially- 
quotient mapping. Then Y is an 1-sequentially-quotient image of a separable metric space 
under the mapping fog. Therefore, Y is an sn-second countable space by Corollary 2.8. 

Finally, we give examples to illustrate the above results. 

Example 2.14. ([2], Example 2.8) A relatively compact metric space is not compact. 

Example 2.15. ([2], Example 2.9) A separable metric space is not relatively compact. 

Example 2.16. A 2-sequence-covering, compact-covering mapping from a separable met- 
ric space is not an mssc-mapping. 

Proof. Let / be the mapping in [2, Example 2.10]. Then / is also 2-sequence-covering. 
This complete the proof. 

Example 2.17. A second countable space is not any image of a compact metric space. It 
implies that “relatively compact metric” in the above results can not be replaced by “compact 
metric” . 

Proof. Let R be the set of all real numbers endowed with the usual topology. Then R is a 
second countable space. Since R is not compact, R is not any image of a compact metric space. 

Example 2.18. A g-second countable space is not any 1-sequence-covering, compact- 
covering compact image of a separable metric space. It implies that “mssc-image” in Corollary 
2.8 can not be replaced by “compact image”. 

Proof. See [22, Example 2.14(3)]. 
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